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AVANT-PROPOS 


Il est assez présomptueux, sans doute, de publier un nouveau traité d'Électromagnétisme, alors que tant : 
d'excellents ouvrages sont déjà disponibles : celui-ci aura l’excuse de sa présentation sous forme de fascicule, 
plus modeste qu’un livre véritable. 

Le public que nous cherchons à atteindre est, avant tout, celui des étudiants de Premier Cycle (MP-PC) 
des Facultés, ou d’établissements de niveau équivalent : nous souhaiterions leur donner ici quelques indica- 
tions sur l’organisation générale de ce cours, et quelques conseils d'utilisation. 

L'’Electromagnétisme dérivé des équations de Maxwell suppose de fréquentes manipulations d’analyse 
vectorielle : nos étudiants, malheureusement, y sont en général mal préparés. Nous avons tenté de combler 
cette lacune par deux chapitres de « calcul » (les deux premiers) les exercices donneront, d'autre part, de 
multiples occasions de s'entraîner. 

Une fois introduites les équations de Maxwell, intérêt se concentre dans deux directions principales. 

L'une, plus « théorique », vise à montrer l’incohérence de l’Électromagnétisme et de la Mécanique clas- 
sique, et la nécessité « d'inventer » la Relativité. 

L'autre développe cette importante conséquence des équations de Maxwell que constituent l’existence 
etla propagation des ondes électromagnétiques. Elle est l’occasion d’une confrontation suivie avec le problème 
général de l’utilisation, en Physique, d’un formalisme un peu élaboré. L'expérience générale a montré que, bien 
que les programmes ne fassent mention que de la propagation dans le vide, il est pratiquement inévitable qu’on 
remplisse un peu ce vide avant d'y propager des ondes, faute de quoi les exercices dont on peut 
disposer manquent par trop de variété. Nous avons cru sage de tirer les conséquences de cet état de fait, et de 
consacrer deux importants chapitres (chap. 11 et 12 du deuxième fascicule) à la propagation dans des milieux 
linéaires. (Nous excluons toutefois toute justification générale des équations constitutives à partir des 
relations microscopiques.) 

Il a semblé bon, pour la cohérence de l’ensemble, et pour fournir aux étudiants parfois pointilleux que 
forment nos collègues mathématiciens les justifications qu'ils affectionnent, d'introduire quelques développe- 
ments complémentaires, tant en ce qui concerne les outils mathématiques (nous disons quelques mots des 
distributions) et les calculs fondamentaux de la théorie (potentiels retardés), que certaines de ses conséquences 
où la nécessité de faire des approximations et de les discuter est l’occasion de retours en arrière qui devraient 
être fructueux (dipôle oscillant ; potentiels de Liénard et Wiechert). Ces compléments sont présentés en petits 
caractères, et pourront être omis en première lecture. 

Ce fascicule dérive d’un « polycopié » utilisé, pendant trois ans, à la Faculté des Sciences de Paris, et qui 
s'était au départ beaucoup inspiré du cours de Physique Fondamentale du Professeur P. NOZIÈRES. Il s’est 
enrichi de remarques et exercices dus en particulier à J. L. BERDOT, M. CHERKI, A. DESNOYERS, B. PERRIN, 
J. F. PÉTROFF et J. TUCHENDLER ; qu’ils trouvent ici mes sincères remerciements. 


INTRODUCTION 


H. 


BUT GÉNÉRAL DE L'ÉTUDE QUI SUIT. 


Notre propos sera ici de donner les bases de la description classique (*) de l'ensemble des phénomènes électro- 
magnétiques. 

Cet exposé a été préparé par celui de l’Électrostatique (2) (qui ne concerne que des charges au repos et leurs effets), 
et de la Magnétostatique (°) (qui ne concerne que les courants stationnaires dont l'intensité et la répartition dans l’espace 
ne dépendent pas du temps.) Nous chercherons désormais à décrire une situation plus générale où les charges électriques 
qui interviennent peuvent se déplacer, et de manière absolument quelconque. Notons que cet électromagnétisme 
des régimes variables a fait une première apparition quand nous nous sommes intéressés au phénomène d’induction 
électromagnétique (Ÿ). 

La théorie classique de l'Électromagnétisme repose sur les célèbres équations de Maxwell, proposées par leur auteur 
voici maintenant une centaine d'années : leur importance mérite amplement qu'on en fasse le titre de ce fascicule. 

Notons dès maintenant que, comme c’est toujours le cas pour une théorie physique, les équations de Maxwell 
ne sont pas « démontrables », de manière logique, ou axiomatique; au moment de leur publication, elles étaient seulement 
cohérentes avec les résultats connus, en particulier d’Électrostatique, et de Magnétostatique, mais « extrapolaient » 
largement. Au reste, leur caractère « esthétique » (qui s’affirmera dans leur formulation relativiste (*)) ne saurait suffire 
à convaincre de leur validité. Elles furent d’ailleurs accueillies avec réticence par beaucoup. Ce type de conflit ne peut être 
résolu que par un verdict expérimental. Ce fut Hertz qui, finalement, força l'adoption de la théorie de Maxwell en 
mettant en évidence l'existence et la propagation d’ondes électromagnétiques, à une vitesse égale ou voisine de la vitesse 
de la «lumière» (les radiations lumineuses n'étant d’ailleurs qu’un cas particulier des ondes électromagnétiques), 
phénomènes qui avaient été prévus par le physicien anglais. (En même temps qu'elles imposaient un progrès théorique 
décisif, les expériences de Hertz fondaient la radioélectricité, dont il est inutile de souligner l’importance pratique.) 

Rapidement cependant (en une quinzaine d'années), le succès même de la théorie de Maxwell plongeait la physique 
moderne dans la première de ses grandes crises. Il apparut en effet que la Mécanique newtonienne classique était 
incompatible avec l'Électromagnétisme. Cette incompatibilité se manifestait en particulier au niveau de l’existence 
d’une vitesse « absolue » (gardant la même valeur dans tous les repères) : la vitesse c de la lumière. La constance de c, 
conséquence des équations de Maxwell expérimentalement confirmée, contredit les formules newtoniennes de composition 
des vitesses. Toutes ces difficultés furent finalement résolues grâce à Einstein et à sa théorie de la Relativité Restreinte. 
Une nouvelle mécanique se fondait, compatible avec l’électromagnétisme de Maxwell, et dont la mécanique newtonienne 
représente une approximation quand les mobiles considérés gardent des vitesses faibles devant c. 

Il serait très souhaitable que l'étudiant qui s'initie aux fondements de l’Électromagnétisme et de la Relativité, 
et qui en trouvera exposés dans des ouvrages tels que celui-ci les aspects en quelque sorte techniques, prenne le temps 
de s’informer parallèlement sur cette phase de l’histoire des Sciences, passionnante, et très révélatrice du fonctionnement 
de « l'esprit scientifique » ($). 


LE PROBLÈME DES CHANGEMENTS DE RÉFÉRENTIELS, 


Une théorie physique se propose de corréler les valeurs d’un certain nombre de grandeurs qui décrivent l’état et 
l’évolution des systèmes physiques. 

D’une manière générale, ces grandeurs sont définies par rapport à un «référentiel» (on dira aussi « repère », 
«système d’axes », « système de coordonnées »), et leur valeur dépend du référentiel choisi. Le mouvement d’un mobile, 
par exemple, peut être examiné par plusieurs observateurs liés à différents référentiels; si ces référentiels sont mobiles 
les uns par rapport aux autres, la vitesse du mobile changera avec le référentiel : la mécanique classique nous a familiarisés 
avec cette possibilité de travailler dans différents repères, et avec cette dépendance des grandeursphysiques vis-à-vis du repère 
où l’on se place. 

Cependant, le repère n’est qu’un intermédiaire : le système physique et la succession de phénomènes dont il est le siège 
ont un caractère « intrinsèque », c’est-à-dire indépendant du repère utilisé pour les décrire. Pour s’assurer que les 
prédictions de la théorie sont effectivement indépendantes du repère utilisé pour la mettre en œuvre, le mieux est de 
travailler simultanément dans plusieurs référentiels, puis de comparer les résultats obtenus et de montrer qu'ils sont bien 
les représentations d’un même phénomène. 


(1) Classique est ici synonyme de non quantique. 

(2) Cf. par exemple, dans la même collection, le fascicule Électricité-2; Électrostatique. Électrocinétique. 

(3) Cf. p. ex. le fascicule Électricité-3, Électromagnétisme. 

(4) Cf. p ex. dans la même collection, le fascicule Relativité Restreinte par Y. SIMON. 

(5) Des éléments très intéressants se trouveront, en particulier, dans les ouvrages de G. BACHELARD tels que L'activité 
rationaliste de la physique contemporaine (P.U.F., Paris, 2° édit., 1965). 
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Outre cet intérêt théorique, l’utilisation en parallèle de plusieurs systèmes de coordonnées peut répondre à une nécessité 
pratique. Si un centre spatial terrestre cherche à téléguider une rencontre de deux vaisseaux cosmiques, il calculera leurs 
trajectoires dans un repère lié à la Terre, ou à l’ensemble Terre-Lune, etc.; mais il devra donner aux équipages, pour 
qu'ils manœuvrent correctement, des indications aussi simples et immédiatement utilisables que possible, et, en premier 
lieu, relatives à un repère lié à la cabine où ils sont enfermés. 

Pour toutes ces raisons, il est indispensable de savoir procéder à un changement de référentiel, c'est-à-dire de pouvoir 
calculer les valeurs de grandeurs physiques dans un repère (5%) quand on les connaît dans un repère (S). 

Nous garderons cette exigence présente à l'esprit dans les développements qui suivent : en même temps que nous 
apprendrons à calculer les champs électrique ou magnétique dans un repère donné en fonction des charges et des courants 
dans ce repère, nous chercherons à définir le mode de transformation de ces champs dans un changement du référentiel. 
Cependant, ce programme ne pourra être rempli qu'après que nous aurons abordé l’étude de la Relativité Restreinte, tant celle-ci 
est intimement liée à l’Électromagnétisme. 


INVARIANCE DES GRANDEURS PHYSIQUES. 


Une théorie physique s'exprime par la définition de certaines grandeurs physiques et des êtres mathématiques qui 
les représentent, et par l'affirmation de «lois», c'est-à-dire de relations entre ces grandeurs, relations formulées 
mathématiquement. 

Ainsi peut-on dire que la Mécanique classique fait intervenir, pour un mobile ponctuel dans un repère galiléen (S) 


le rayon-vecteur (P), le temps (^, la vitesse (8), l'accélération (y), la masse (m), la force subie (È, etc... avec les « lois » : 


— decinématique : 
D = dř/dt 7 = dë/dt (Électromagn. int. 1) 
— de dynamique 
F = mÿ (Électromagn. int. 2) 


Nous savons changer de repère (au moins de repère galiléen). Si de (S) nous passons à (S'), en translation à la 
vitesse ï par rapport à (S), le temps reste £ = t; le rayon vecteur devient : 


F — üt 3 
e | (Électromagn. int. 3) 


i + 


la vitesse : v=0—- ü 
Paccélération reste ÿ = Y, la force F = F, la masse m' = m et l’on a encore : 
D = dr/dt 7 = dÿ/d” P = m7. 


Autrement dit, les relations fondamentales de la Mécanique classique, (Int. 1) et (Int. 2) sont invariantes dans un 
changement de repère galiléen. 

Ce faisant, la Mécanique classique va au-delà de l'exigence formulée plus haut de se prêter aux changements de 
référentiels. Elle satisfait à un principe plus général (dont la nécessité s’est d’ailleurs imposée dans la deuxième moitié du 
xx" siècle, à l’époque même des travaux de Maxwell), qu’on peut appeler principe d'invariance, ou principe de relativité, 
et qu'on peut énoncer ainsi : 

à une théorie physique donnée, on peut associer une classe de référentiels (C) telle que les lois exprimant cette théorie 
aient même forme mathématique dans tous les référentiels de (C). 

Dans le cas de la Mécanique classique, (C) est constituée par l’ensemble des repères galiléens. Il n'est pas possible de 
distinguer deux repères galiléens par une expérience de Mécanique classique, car les lois de ladite mécanique y sont les 
mêmes. (On sait par contre que les effets de forces d’inertie permettent de reconnaître un repère non galiléen d’un repère 
galiléen.) Autrement dit, il n'existe pas de repère galiléen privilégié, « absolu »; tous les repères galiléens sont physiquement 
équivalents (ce qui ne veut pas dire également adaptés à la description d’un système donné), et l’on ne peut mettre en 
évidence que des mouvements relatifs de ces repères. D'où le nom de principe de relativité (1). 

La dénomination « principe d’invariance » se justifie d'elle-même. Comme les formules de changements de repère à 
l'intérieur de (C) forment un groupe, soit G, on dira que G est le groupe d'invariance de la théorie. Ainsi la Mécanique 
classique admet pour groupe d’invariance le groupe des transformations galiléennes (Int. 3). 


(1) Quand (C) comprend une partie seulement des repères possibles, par exemple les repères galiléens, on dit que la théorie 
satisfait à un principe de relativité restreinte. Sinon, il y a « Relativité générale ». 
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L'Électromagnétisme, fondé sur les équations de Maxwell, doit lui aussi satisfaire à un principe d’invariance : par un 
procédé d’ailleurs détourné, nous verrons ultérieurement que ses équations sont invariantes dans le groupe des « trans- 
formations de Lorentz », qui est aussi le groupe d’invariance de la Relativité Restreinte. (Ici encore se manifeste l’étroite 
liaison entre électromagnétisme et Relativité.) Par suite, les groupes d’invariance de la Mécanique classique et de 
l’Électromagnétisme de Maxwell ne coïncident pas : Cest wne nouvelle expression, plus abstraite, de leur incompatibilité 


déjà signalée. 


FORMULATION COVARIANTE DUNE THÉORIE PHYSIQUE. 


Une fois admise la nécessité, en fondant une théorie physique, de préciser son groupe d’invariance, il devient avantageux 
de donner à cette théorie une formulation mathématique qui rende évidente cette invariance. 

Nous rappellerons, au chapitre 1, que la notion de vecteur est directement liée au comportement des composantes 
dans une rotation du système d’axes. Soit A et B deux vecteurs, qui dans un premier repère (S), ont pour composantes Ay, Ay, 4, 
et Bo B,, B, L'égalité : 


Â=8B (Électromagn. int. 4) 


nous indique que l’on a, dans (S) : | 
(Électromagn. int. 5 


Si de (S) nous passons à un autre repère (S'), les composantes de A et B deviennent À,, Ap À, et Ba Bp Bz les formules 
de passage des A; aux A; coïncidant avec celles qui permettent le calcul des B; à partir des B;. De (Int. 5) on déduit donc : 


A&A =B, ; 4 =B, ; A =B, (Électromagn. int. 6) 


Autrement dit, les formules (Int. 5) sont invariantes dans le groupe des rotations à trois dimensions. Mais légalité 


vectorielle (Int. 4) suffit en fait à assurer cette invariance, qui est pour ainsi dire contenue dans la définition de À et B 
comme vecteurs. 
Il en serait de même d’une égalité 


A=B 


entre deux quantités dont on aurait vérifié qu’elles sont des scalaires (+). 

Il apparaît (et c’est le cas pour l’Électromagnétisme et la Relativité Restreinte) que certains changements de repère 
peuvent s'identifier à des rotations de systèmes d’axes (au besoin dans des espaces à plus de trois dimensions : ainsi 
« l'espace-temps », à quatre dimensions, de la Relativité Restreinte), et que les grandeurs physiques que décrit la théorie 
peuvent s’interpréter comme l'équivalent, dans ces espaces, soit de scalaires, soit de composantes de vecteurs (ou de tenseurs). 
Si les lois de la théorie s’expriment par des égalités entre ces « scalaires », ces « vecteurs » (ou ces « tenseurs »), on est alors 
assuré de leur invariance dans les changements de repère en question. On dit que la théorie est mise sous forme covariante. 

Nous indiquerons (fascicule « Relativité Restreinte ») comment des équations équivalentes aux équations de Maxwell (2) 
peuvent être mises sous forme covariante : on appréciera alors le caractère synthétique de la formulation obtenue. 


(1) Plus généralement l'invariance serait automatiquement assurée pour une égalité entre « tenseurs ». Cette notion sera 
mentionnée au chap. 1. 
(2) C'est d’ailleurs aussi le cas des équations de Maxwell elles-mêmes. 


RAPPELS MATHÉMATIQUES - 1 


Les notions rappelées dans ce chapitre ne doivent pas être entièrement nouvelles ; l'étude de l’électrosta- 
tique et de la magnétostatique a fourni de multiples occasions de les mettre en œuvre. Ce chapitre n’est donc 
qu'un résumé, et non une introduction ({). 


I. REPÉRAGE D'UN POINT DE L'ESPACE A TROIS DIMENSIONS. 


Un point M de R? peut être repéré par ses trois coordonnées cartésiennes, x, y, z, dans un système d’axes 


qu’on supposera orthonormé (sauf mention contraire) : (Ox, Oy, O2). 
Pour tenir compte de symétries éventuelles, cylindrique, sphérique, etc..., il peut être avantageux d’adop- 
ter des coordonnées curvilignes, en particulier : 


1° cylindriques : 


On remplace les deux coordonnées cartésiennes x et y par les coordon- 
nées polaires, p et 0, dans le plan (Ox, Oy). On a donc, m étant la projection 
de OM sur (Ox, Oy) 


0 = (Ox, Om) = Arctg(y/x) nu. 


Élect ciei 
p = Om spp y=p'sin0 (Électromagn. 1-1) 


Les trois coordonnées cylindriques de M sont p, 8 et z. (L’axe Oz est dit 


axe des coordonnées cylindriques considérées.) at 
IG. 1. 


2° sphériques : 


avec les notations adoptées pour les coordonnées cylindriques, on repère M 
par : 
r= OM = (x? + y? + 22)" 


6 = (Oz, OM) en prenant ici 0<0<7r;0 = Arc tg£ 


p = (Ox, Om) = Arc tg (y/x) 


(Électromagn. 1-2) 


Fc. 1.2 


On a donc : 


z=r"cos 0 p=r'sinð >: 
(Electromagn. 1-3) 


x = r sin 0 cos @ y = rsin sin @ 


On peut enfin fixer la position de M par la donnée de son rayon-vecteur, à partir d’une origine que nous 
noterons O (et qui sera par exemple confondue avec l’origine d’un système d’axes). Cerayon-vecteur sera noté : 
F, AM), Fy, OM, M... 

On remarquera que le module Fi — OM de ce rayon-vecteur coïncide avec la coordonnée sphérique que 
nous avons notée r. On peut poser : 

F=r (Électromagn. 1-4) 


Supposons que, de M, nous passions au point infiniment voisin M’; F subit un accroissement que nous 
pouvons noter dr : 


d? = MM' 


(1) On trouvera plus de détails dans notre fascicule Rappels et compléments de mathématiques, dans la même collection. 
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Mais r = OM devient, dans le même temps, r = OM’ =r + dr. 
dr = OM' — OM 


est en général distinct de MM' = |MM| = |df. 

(En fait, si K est l'intersection avec OM du cercle de centre O passant 
par M',ona :dr = MK = MH où H est le pied de la perpendiculaire abaissée 
de M' sur OM, le calcul étant limité au premier ordre en dr ou |dř].) Diffé- 
rentions d’ailleurs : 


il vient : 


dr = = ü dF (Électromagn. 1-5) 


(où ü est le vecteur unitaire de OM); le produit scalaire à - dF représente précisément MH. 
Revenons en coordonnées cartésiennes; on a : 
r = x? + y? + 7? 
A y et z fixés, x variant de dx, il apparaît une variation dr de r avec 


r'dr = x:dx. 


On en déduit ELA == 
ôx r 
et de même CRE (Électromagn. 1-6) 
ðy r 
Ca 
ôz r 


VECTEURS ET CHAMPS DE VECTEURS. 


Il est inutile d’insister sur la définition géométrique du vecteur, donné par son origine et son extrémité. 
Partant de cette définition, on peut : 


— repérer le vecteur par ses composantes, par exemple par rapport à un système d’axes cartésiens. 


— donner les « formules de changements d’axes » qui permettent de passer des composantes par rapport 
à un premier système d’axes auxcomposantes par rapport à un second système déduit du premier par une rota- 
tion, une symétrie-plan, une inversion, etc... 


Inversement, si un procédé de calcul, un algorithme, permet de calculer trois grandeurs x, y, z dans un 
premier système d’axes (S), et trois autres grandeurs x’, y’, z' dans un deuxième système ($°), et si l’on passe 
de (x, y, Z) à (x’, y’, zZ’) quels que soient (S) et (S”) par les formules de changement d’axes valables pour les 
composantes d’un vecteur défini géométriquement, on dira que l’algorithme définit un vecteur, de compo- 
santes (x, y, Z) ou (x’, y’, z’) suivant le système d’axes considéré, (et alors même qu'on a perdu, pour la déter- 
mination de ce vecteur, toute référence géométrique). 

Notons : 


— que cette définition du vecteur dans R? sera généralisée par la suite (cf. Fascicule « Relativité 
Restreinte ») : on travaillera alors dans un espace à quatre dimensions, sur des « quadrivecteurs ». On perd 
la référence géométrique; l’analyse du comportement dans un changement de « repère » devient essen- 
tielle. 


RAPPELS MATHÉMATIQUES - 1 (SUITE) 


— que, dans R? même, le vecteur n’est qu’un cas particulier d’une catégorie mathématique plus vaste, 
le « tenseur ». (Cette notion sera elle-même introduite dans « Relativité Restreinte ».) 


(— qu'un autre type de tenseur, plus simple encore que le vecteur, est le scalaire, autrement dit le nombre, 
réel, complexe au besoin, qu’on calcule par référence à un certain système d’axes, mais dont la valeur est fina- 
lement indépendante du système. Un exemple : la longueur d’un vecteur défini par ses composantes dans un 
système d’axes : Í = (x? + y? + 22)1/2)) 

Outre les « vrais » vecteurs considérés jusqu'ici, on peut, dans R°, introduire des « pseudo-vecteurs » (ou 
vecteurs axiaux) (le champ magnétique est un vecteur axial; le produit vectoriel de deux -« vrais » vecteurs est 
un vecteur axial). (En fait, ces pseudo-vecteurs sont des tenseurs, mais qui n’ont que trois composantes dis- 
tinctes, qu’onidentifie à celles d’un vecteur.) Géométriquement, les pseudo-vecteurs se distingueront comme suit 
des vecteurs au sens propre : 

— dans une rotation propre (rotation d’un certain angle autour d’un certain axe, l’angle et l’axe étant 
d’ailleurs quelconques, opération qui transforme un trièdre « droit » en un autre trièdre droit, et une main 
droite en une « autre » main droite) vecteurs et pseudo-vecteurs se transforment de la même manière. 


— mais dans une rotation impropre (produit d’une rotation propre par l’inver- y” 
sion, ou symétrie par rapport à un point; dans cette catégorie figurent en particulier 
Pinversion ellemême, et les symétries par rapport à des plans, opérations qui trans- 
forment un trièdre « droit » en un trièdre « gauche », et une main droite en une main 
gauche) un vecteur se transforme simplement en son opposé (pour l'inversion), son 
symétrique par rapport au plan de symétrie (pour une symétrie-plan), tandis qu’un y 
pseudo-vecteur se transforme en l’opposé de ce transformé géométrique : il est, par Fic. 1.4 
exemple, invariant dans l'inversion; dans la symétrie par rapport au plan P, si le 
vecteur Ÿ est un « vrai » vecteur, il donne V; s’il est un pseudo-vecteur, il donne P" = — v. 


On convient parfois de noter Ÿ (surmonté d’une flèche droite) un vecteur et Y (flèche courbe) un pseudo- 
vecteur. On omettra au besoin cette distinction, l’essentiel étant de garder présente à l’esprit la nature exacte 
des êtres qu’on manipule, et d'éviter, dans toute la mesure du possible, de soulever des problèmes qui ne 
peuvent être résolus que par le recours à un formalisme plus poussé. Dans l’espace à trois dimensions, on 
bénéficie de la possibilité de traiter certains tenseurs comme des vecteurs tant que n’interviennent que des 
rotations propres, des trièdres qui sont tous de même sens (tous droits par exemple); on s’efforcera de s’en 
tenir là. 

(Par contre, il est essentiel de n’écrire que des équations qui respectent ce qu’on peut appeler l’« homogé- 
néité tensorielle » : un scalaire ne peut être égal qu’à un autre scalaire, un vecteur à un autre vecteur. C’est 


une faute très grossière d'écrire, s étant un scalaire, et un vecteur :s = Ñ. Si l’on distingue vecteurs et pseudo- 
vecteurs, il ne peut de même y avoir égalité qu'entre deux vecteurs ou deux pseudo-vecteurs, mais pas entre un 
vecteur et un pseudo-vecteur.) 

— Nous n’avons mentionné jusqu'ici que les composantes cartésiennes des vecteurs par rapport à un 
repère (système d’axes) fixé une fois pour toutes, (pour un problème donné, s’entend). Or on a fréquemment à 
considérer des champs de vecteurs, c’est-à-dire des vecteurs fonctions de point (ou des fonctions de point à 
valeurs vectorielles). Soit M) un tel champ de vecteurs : aux points M d’un certain domaine de l’espace, on 
fait correspondre les vecteurs Į (eux-mêmes, en général, dans un espace à trois dimensions). 

Il peut être utile, comme on l’a vu, de repérer M par des coordonnées curvilignes, cylindriques ou sphé- 
riques par exemple. Dans ce cas, on a en général intérêt à donner les composantes du vecteur V(M) associé 
à ce point non pas dans un repère indépendant de M, mais dans un système d’axes orthonormé, S(M), qui est 
lui-même fonction de M. 

Ce repère est défini de la manière suivante. Soit x,,x,, X3 les trois coordonnées curvilignes de M, par 
exemple : 


pre 
V 


X,=P ; X% =0 , x}=Z (coordonnées cylindriques) 


x” =r , X1=0 , x; =@ (“en sphériques”) 


Les conditions 
x, = x (M) = constante 


xa = X{M) 
x3 = X3(M) 
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définissent trois surfaces qui se coupent en M : 


— en coordonnées cylindriques x, = p = constante correspond au cylindre ayant pour axe l’axe Oz des 
coordonnées, qui passe par M; x, = 0 = constante correspond au plan (M, Oz}; x, = z = constante, au 
plan perpendiculaire à Oz passant par M. (Fig. 1.5.) 


— en coordonnées sphériques x, = r(M) est la sphère de centre O passant par M; x, = 0 = OM), le 
cône d’axe Oz et demi-angle au sommet 4; x, = @ = œ(M) le plan « méridien » (Oz, M). (Fig. 1.6.) 


Fic. 1.6 


Les conditions x; = x;(M); x; = x (M) définissent une courbe C, passant par M, qui est l'intersection de 
deux des surfaces précédentes, et le long de laquelle seule la troisième coordonnée, x,, varie. Soit t, le vecteur 
unitaire de la tangente à C, en M. 

L'ensemble des trois vecteurs i4, ,, U, définit le système d’axes S(M) que nous avions annoncé. On véri- 
fiera immédiatement que, pour les deux types de coordonnées curvilignes que nous avons définis, S(M) est 
orthonormé (1). , 

C’est par ses composantes dans S(M) qu’on pourra définir le vecteur V(M). On les notera K, Vp, Vy 
ou V, K etc... de manière à rendre plus explicite l'orientation de l’axe correspondant : 

— en coordonnées cylindriques, C, = C, est laxe HM intersection du plan méridien de M et du plan 
perpendiculaire à Oz passant par M, autrement dit l'axe perpendiculaire à Oz mené par M. La composante 
correspondante V; est dite radiale. C, = C, est le cercle d’axe Oz passant par M; la composante V, sera dite 
tangentielle. (C, = C, est la parallèle à Oz menée par M.) 

— encoordonnées sphériques, C, = C, est laxe OM ; la composante correspondante F; sera dite radiale. 
C, = Cest le cercle « méridien » de centre O, passant par M, situé dans le plan (M, Oz). C, = C, est le cercle 
« parallèle » d’axe Oz passant par M. 

(Dans le même contexte, on peut remarquer que la variation élémentaire dx, de x, à partir de x (M), 
x; et x; restant fixes, fait passer de M en M,, point infiniment voisin sur C, Au premier ordre 
en dx,, On à : 

MM, = h, 1: dx, (Électromagn. 1-1) 


où h, est un paramètre (en général fonction de x, x; x,) qui permet de passer de la variation dx, de x, à la dis- 
tance parcourue sur C,. Ainsi : 


— «en cylindriques» : A, = h, = 1; h, =h=p;h;=h,=1 


zZ 


— «en sphériques » : h, = h, = 1; h, = hy = r; h, = h,= r'sinð. 


Dans (1-7), MM, est la variation du rayon vecteur F = F{x,, X2, x.) de M qui accompagne la variation 
de x, à x; et x; fixés. Si les trois coordonnées curvilignes subissent de « petits » accroissements, 7 subira l’ac- 
. -b Eggs = . . a . . 
croissement dř somme des MM, ainsi calculés (toujours au premier ordre par rapport aux dx,) : 


3 
d= X hy't’ dx {Électromagn. 1-8) 
k=1 


(1) Un choix convenable des orientations permet à S(M) d’être direct. 


RAPPELS MATHÉMATIQUES - 1 (SUITE) 


HT. 


Rappelons enfin que les champs de vecteurs interviennent fréquemment dans le calcul de deux scalaires : 
— la circulation sur une courbe C (joignant deux points À et B, ou courbe fermée) : 


g = f V(M): dM 

+ (Électromagn. 1-9) 
€ = p (M) - dM | 

c 


(le point M, auquel est associé le vecteur V(M) du champ de vecteurs considéré, décrit la courbe C :dM est 
l’accroissement de son rayon-vecteur pour un déplacement infiniment petitsur €. Ilest essentiel de remarquer 


que F et dM interviennent par leur produit scalaire; il faut donc connaître non seulement les modules de ces 
vecteurs, mais aussi leur angle, c’est-à-dire l'angle de la tangente à C en M avec }(M).) 
— le flux à travers une surface S (limitée par une certaine courbe C, ou fermée) : 


ð= [ V(M)-dS (Électromagn. 1-10) 
(5) 


(le point M de S, auquel est associé Ÿ(M ), est le point moyen d’un élément de surface d’aire dS; dS est un vec- 
teur normal à S en M et de longueur égale à dS. Ici encore, les questions d'orientation sont primordiales (!)). 
(Dans le cas d’un champ de vecteurs « à flux conservatif», dont la « divergence » est nulle (vide infra), ® ne 
dépend pas de S mais seulement de la courbe C qui limite cette surface, si elle n’est pas fermée. Dans ce cas, et 
dans ce cas seulement, on pourra parler du flux à travers la courbe C — et non à travers la surface S.) 


GRADIENT DUNE FONCTION SCALAIRE. 


Soit (M) une fonction de point à valeurs scalaires (on pourrait dire un « champ scalaire », par analogie 
avec un champ de vecteurs) définie sur un certain domaine de R°. 

On peut repérer M par trois coordonnées cartésiennes, x, y, Z, relatives à un repère orthonormé; {(M) 
devient alors une fonction de trois variables : 


JM) = f(x, y z). 


Ses dérivées partielles constituent les composantes, dans le repère orthonormé utilisé, dun champ de 
vecteurs : le « gradient de f », noté grad f : 


_ ôf _ ôf _ ôf ; 
grad, f É à grad, f = y grad, f = a7 (Electromagn. 1-11) 


(Nous avons ici, notons-le au passage, un bon exemple d’une définition non géométrique d’un vecteur. 
Les trois dérivées partielles de M), dans un autre repère orthonormé où les coordonnées de M 
sont x’, y', z', deviennent 


eA e UE 
ôx' °? ôy 87 
Le à | , of of 
On peut définir un vecteur gradient parce qu’on passe des 3x AUX = par les formules de changement 


d’axes caractéristiques des composantes d’un vecteur; la démonstration de cette possibilité est ici laissée au 
lecteur (!).) 


(1) Cf. Rappels et Compléments de Mathématiques, chap. 10. 
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Si les coordonnées de M subissent des accroissements dx, dy, dz (qu’on peut considérer comme compo- 
santes de la différentielle dM du rayon vecteur de M), f subit une variation 


-2 ax p Lay p 2L. 
DÉS PRO: az dz 


(d’après la définition des dérivées partielles), ce qu’on peut récrire 


| df = grad f: dM | (Électromagn. 1-12) 


Cette « formule du gradient » peut d’ailleurs être considérée comme définissant ce dernier, et de manière 
«intrinsèque », indépendamment de la référence à un système d’axes. On en déduira très commodément les 
composantes du gradient dans le système cartésien local rattaché à un ensemble de coordonnées curvilignes 
(système d’axes S(M) du $ II° de ce chapitre). Si grad, f est la composante suivant l’axe de vecteur unitaire i, 
nous tirons de (1-8) et (1-12) : 


3 
df = Ÿ h,:grad, f dx, 
k=1 


En identifiant à 


3 of 
df = ), z>; dx 
f 2, ÔXy g 
Nous obtenons : 
1 ôf z 
grad, f = =: 5 (Electromagn. 1-13) 
k OXk 
d’où : 
— en coordonnées cylindriques 
of 
grad, f = P 
„3S Lo 
grad, f = — 36 grad, f = En 
— en coordonnées sphériques : 
à 
grad, f = a 


1 ô 1 ôf 
grady f = 1.25 grad, f = gn g 


Terminons en signalant une condition nécessaire et suffisante pour qu'un champ de vecteurs YMD) soit le 
gradient d'une fonction scalaire : que sa circulation sur une courbe fermée quelconque soit nulle (ou, de manière 
équivalente, que sa circulation entre deux points 4 et B soit fonction seulement de ces points et pas du trajet 
les joignant, et le long duquel on calcule la circulation). 

Si cette condition est satisfaite, on peut noter : 


V(M) = — grad f(M) (Électromagn. 1-14) 


et on dit que Ÿ dérive du potentiel scalaire f. (L'introduction du signe — est ici traditionnelle.) 


RAPPELS MATHÉMATIQUES - 1 (SUITE) 


IV. DIVERGENCE D'UN CHAMP DE VECTEURS. 


Soit un champ de vecteurs M); si l’on dispose d’un repère orthonormé dans lequel M a pour coor- 
données x, y, 2, et F pour composantes K, V et F (qui sont trois fonctions scalaires de x, y et z), on définira le 


scalaire div Ÿ par : 


PRE: z 
es 7. 
div + à (Électromagn. 1-15) 


div F est bien un scalaire car la valeur donnée par la formule précédente, pour un champ de vecteurs donné, 
en un point M donné, est indépendante du repère choisi pour mener le calcul. 

L'introduction de la divergence est justifiée par le théorème suivant, connu sous des noms divers, et que 
nous appellerons théorème de la divergence : si F est une surface fermée quelconque qui enferme un volume #”, 


le flux sortant du champ de vecteurs Ñ à travers 4 est égal à l'intégrale triple de la 
divergence de V sur le volume # : 


j -dS = [ffa Ÿ : dr (Électromagn. 1-16) 
sg r 


(le flux sortant est calculé avec des normales dS orientées de l’intérieur vers l'extérieur de la surface. Le 
champ Ÿ est continu et dérivable dans #) 

Cette propriété fondamentale constitue une sorte de définition « intrinsèque » de la FEES indépen- 
dante de tout système d’axes. 

Supposons que (M) soit défini par ses composantes V, dans le système cartésien Mere M associé 
à un système de trois coordonnées curvilignes x}, x, et x. Consi- 


dérons le volume limité par les six surfaces : 
(1) (x; = x,(M) + dx;) et 09) Gi = x, (M); 


(2) (x, = xM) + dx,) et (2) (x: = x (M); 
(3) (x3 = x3(M) + dx3) et (3) (C3 = xM). 


Le flux de Ÿ à travers la surface 4 formée par ces six surfaces 
comporte : FIG. 1.8 

— une contribution de (1), produit de laire de cette face : h,h, dx, dx, par la composante normale 
de V : V; (orientée suivant la normale sortante à S), ces quantités étant calculées pour x, = x, (M) + dx., 
x.(resp. x3) = x:(M)(resp. x:(M)). à 

— une contribution de (l'), produit de laire h,h dx, dx,, par la composante de V dirigée suivant la 
normale extérieure à & soit — K; Vi, hz et h, étant ici pris pour x, = x, (M). 

Combinant ces deux expressions, et retenant les termes du premier ordre en dx,, nous obtenons 
ð 


(Vihohs | A 


"dx, dx, — (Vihihs |.) dx, dx, = 


sas) 


Ajoutant les expressions analogues correspondant aux faces (2) et (2°), (3) et (3°), nous obtenons le flux, 
élémentaire, de V à travers S: 


do = (J3 (hh h) dx, dx, dx. 
k 
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D’après le théorème de la divergence, cette expression peut être identifiée au produit de div Ÿ par le 
volume enclos par, soit à l’ordre le plus bas en dx; : h hh; - dx, dx, dx. Il vient finalement : 


. à 1 ð ô ð 2 
div Ÿ = hhohs LE (Vh;h;) + EPA (Vahh) + EA (Vahat) | (Electromagn. 1-17) 


(Pour des coordonnées cartésiennes, tous les A; sont égaux à 1, et l’on retrouve (1-15).) 
En coordonnées cylindriques, on aura : 


1 à 18V, ôV, s i 
div V = Sap O% + 5 4 Fe (Électromagn. 1-18) 
En coordonnées sphériques, on aura : 
div Ÿ = ii + 1n + le + se (Électromagn. 1-19) 
r? ôr r 00 r‘tg0 rsinm® à 


Terminons par une formule d'emploi très fréquent. Si V est un champ de vecteurs, f un champ scalaire, 
(fV) est un nouveau champ de vecteurs, et l’on a : 


div (FF) = f -div Ÿ + Ÿ:grad f (Électromagn. 1-20) 


ROTATIONNEL D'UN CHAMP DE VECTEURS. 


—— > 


Reprenons les notations du début du § IV. On appellera « rotationnel de V», noté rot Ÿ, 
vecteurs de composantes : 


(Électromagn. 1-21) 


= 5 
Remarque. Si V est un vecteur, rot V est urni pseudovecteur, et réciproquement. Si V est un vecteur, dans 


. . =i Là La d fak 
l'inversion par exemple, les V changent tous de signe; les dérivées z 


i 


ð ô ; 

placées par -~~ = — ——.et les composantes du rotationnel, au total, ne 
a- xi) ôx; 
changent pas : il s’agit bien d’un pseudovecteur). 

Le rotationnel satisfait à une importante propriété, que nous appellerons 
simplement le théorème du rotationnel. Si € est une courbe fermée quelconque, 
F une surface, quelconque, s'appuyant sur ®, la circulation du vecteur V sur € 
est égale au flux de son rotationnel à travers S, dont la normale positive est orientée par référence 
au sens de parcours sur € : 


Fic. 1.9 


| $ Ÿ -dM = f rot ÿ dS | (Électromagn. 1-22) 
d L 


Cette relation joue, pour le rotationnel, un rôle analogue à (1-16) pour la divergence : on peut, en parti- 
culier, la considérer comme une véritable définition intrinsèque du rotationnel, et lappliquer, par exemple, 
au calcul des composantes du rotationnel associées à des systèmes de coordonnées curvilignes. 


RAPPELS MATHÉMATIQUES - 1 aa 


Signalons deux propriétés fondamentales : 
1° La divergence d’un rotationnel est identiquement nulle 


— 


| div (rot V) = 0 | (Électromagn. 1-23) 


(la propriété se démontre en combinant (1-21) et (1-15)). 

Il en résulte que le rotationnel d’un champ de vecteurs quelconque est « à flux conservatif » : son flux à 
travers une surface fermée est identiquement nul; son flux à travers une surface non fermée, s’appuyant sur 
une courbe fermée, ne dépend que de celle-ci, et pas de la surface. On remarquera que cette particularité est 
nécessaire pour que le théorème du rotationnel ait un sens. 

Réciproquement, si A est un champ de vecteurs de divergence nulle, on peut l’exprimer comme le rota- 


tionnel d’un autre champ de vecteurs, soit Ñ: 


dv =0 + Å= | (Électromagn. 1-24) 


On dit que À dérive du potentiel vecteur V, 


2° Le rotationnel d’un gradient est identiquement nul. (Combiner (1-11) et (1-21).) Réciproquement, 
si le rotationnel d’un champ de vecteurs V est identiquement nul, V dérive d’un potentiel scalaire : 


rot (grad f) = 0 (Électromagn. 1-25) 
tř =0 e V= - grad f (Électromagn. 1-26) 


3° On démontre que tout champ de vecteurs Ÿ peut s'exprimer comme somme d’un gradient et d'un 
rotationnel : on peut trouver f et À avec 
Ÿ = rot À — grad f (Électromagn. 1-27) 


Si div Ÿ = 0, le dernier terme disparaît: si rot Ÿ = 0, il subsiste seul. (On remarquera que la décomposi- 


tion (1-27) d’un champ de vecteurs Ÿ donné ne saurait être unique : l'addition à f d’une constante, à 4 du 
gradient d’une fonction quelconque ne modifie pas le gradient ou le rotationnel qui apparaissent dans l’expres- 


sion de F.) Un champ de vecteurs Ÿ étant donné, le calcul immédiat de sa divergence et de son rotationnel 
s'impose : on saura ainsi si Ñ dérive d’un potentiel scalaire, ou vecteur, ou s’il a la forme générale (1-27). 
Terminons par une formule que nous aurons l’occasion d'utiliser : si À et B sont deux champs de vec- 
teurs, on a : 
div (Ë a 4) = À (rot B) — B- (rot À). (Électromagn. 1-28) 


VI. LAPLACIEN SCALAÏRE; LAPLACIEN VECTORIEL. 


f étant une fonction scalaire de point, on appelle « laplacien » de f, noté Afla divergence du gradient de f. 
De (1-11) et (1-15), en coordonnées cartésiennes, nous tirons : 


EE 2 2 2 3 
Af = div (grad f) = es + ui + Les (Électromagn. 1-29) 


ôx? | Op" 0 


De (1-18) et (1-19) et de l’expression des composantes du gradient nous tirons de même : 
— en coordonnées cylindriques : 


2 2 
Af = 12( of 4 LOS + (es (Électromagn. 1-30) 
ôp p? 80? ôz? 
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VIT. 


— en coordonnées sphériques : 
2 
Af = T r? ai + L Z (sin 0 21) + l nues (Électromagn. 1-31) 
r? ôr ôr r?» sin 0 ô0 20 r? sin? 0 ôg? 


P étant un champ de vecteurs dont les composantes, dans un certain repère cartésien, sont K, V, Vp 
on appellera laplacien (vectoriel) de’ P, noté AP ou A, le champ de vecteurs dont les composantes sont les 
laplaciens des composantes de P : 


[AP = AF, 


x 


; (AP), = AY, ; (AV), = AY. | (Électromagn. 1-32) 


(on vérifierait facilement que ces relations définissent effectivement un nouveau champ de vecteurs). 
L'introduction du laplacien vectoriel est justifiée par l’identité suivante, fondamentale en électromagné- 


tisme, et valable quel que soit le champ de vecteurs 4 auquel on l’applique : 


rot (rot V) = grad (div Ÿ) — AŸ | (Électromagn. 1-33) 


Nous arrêterons ici la revue des notions et relations, liées aux opérateurs différentiels, que nous aurons à 
mettre en œuvre dans la suite de cet exposé. Remarquons cependant que nous avons, partout, implicitement 
supposé que les grandeurs traitées étaient dérivables, voire continûment dérivables : une singularité d’un 
champ de vecteurs, par exemple, peut grandement perturber sa divergence; le « théorème de la divergence » 
ne sera appliqué, sans précautions spéciales, qu’à des domaines d’ intégration qui excluent de telles singula- 
rités. Il ne s agit pas pour nous d’une simple clause de style : le champ « en 1/r? » del’ Électrostatique présente 
une singularité à l’origine qui retentit de manière très importante sur le comportement de sa divergence; nous 
y reviendrons plus loin. 


FONCTIONS DE LA DISTANCE ENTRE DEUX POINTS. 


Nous aurons fréquemment, dans la suite, à manipuler des fonctions, par exemple scalaires, de la distance 
entre deux points P et M, soit (PM). Ce sera l’origine de transformations qui reviendront de manière suf- 
fisamment systématique pour motiver l’addition, à ce chapitre, de ce court paragraphe. 

La fonction f(r) = f(PM), (PM = r), sera tantôt prise à P fixe, M (de coordonnées xy, Ym, Zm) variable; 
tantôt à M fixe, P (de coordonnées x», Yp, Zp) Variable. Dans le premier cas, nous noterons le gradient de f 
grad,, f, pour marquer que les dérivées qui constituent les composantes sont d//0x,,, etc... Dans le second cas, 
nous noterons le gradient : grad, f [composantes : ES . De plus, P, quand il est fixe, peut être pris comme 

P == , 
origine du système de coordonnées : x}, Yy, Zm Sont alors les composantes de PM (etréciproquement, quand P 


varie, M, fixe, est pris comme origine et xp = (MP), etc). 
ôf df ôr Xy df pe 
Tirons alors parti du fait que f n’est fonction que de PM r o F T (on utilise (1-6), 


et la dérivation des fonctions de fonctions). 


De même : 
OF Xe ÎS 
ôxp r dr’ 
Or 
= (PM), = — (MP), = — Xp 
D'où 
f 2L 
ôxp  ÔXy 


Le même calcul peut être repris pour les autres composantes des gradients. Au total : 


grady f(PM) = — grad, f(PM) (Électromagn. 1-34) 


EXERCICES Rappels mathématiques - 1 


[1-1] Exprimer grad (f + g) et grad (fg), où f et g sont deux fonctions de points à valeurs scalaires, en fonction 
de grad f et grad g. 


a) grad (f + 8) : 
Calculons la composante suivant l’axe (cartésien) Ox : 


ô 
grad, (f +9) =z + g) Ars 
soit 
grad, (f + g) = grad, f + grad,g 


et, la même relation étant satisfaite pour les composantes suivant les deux autres axes Oy et Oz, nous pouvons écrire 
la relation vectorielle : 


grad (f + g) = grad f + grad g 
b) grad (f8) : 
On suit la même méthode : on trouve ainsi : 
2 -ofl p. , 
grad, (fg) = z (o) = ge t Saag ad, f + f:grad,g 


et 
grad (fg) = g: grad f + f grad g 


— _ > -> 


[1-2] Exprimer div (À + Ë) et rot (4 + B) (ou À et Ë sont deux champs de vecteurs) en fonction des diver- 


e 


gences et rotatiormels de À et 8. 


On calculera : 


diva = ZE (A + B) 


(où x, = x, y, z suivant la valeur de i) soit 
> ô ðA; ôB, 
A = = : 7 E pen à EE à 
ay Èx, (AER (5 a + (z a.) 
et finalement : 
div ( + B) = div À + div Ë 
La composante suivant Ox de rot (4 + B) vaut, par définition : 


à + > Ô _- =, 


rot, (A + B) = o (4 kB -= (4 + B), 
OA OA ôB, ôB -> 3 
— DR RER CAES = 
( 3y z7 ) + 5 z) rot, A + rot, B. 


Des relations analogues interviennent pour les composantes suivant les deux autres axes cartésiens Oy et Oz, 
et nous pouvons écrire : 


dag se > — 


rot (À + B) = rot À + rot B. 


EXERCICES | 1 | : 


appels mathématiques - 1 (SUITE) 


[1-3] Démontrer la relation (1-20) : 


div Y- P) = f -div Ÿ + Y- grad f 


Le vecteur f 4 (où f est une fonction, V un champ de vecteurs de composantes V, V, V, dans un système 
d'axes orthonormés), a pour composantes f K, fV, et fV. D'où 


TAE AEE TAER UM) 


Les formules habituelles de dérivation des produits permettent d’écrire : 


ð pi y. pN 
a UV) R T Fa ôx ôx 
ð ð ôV 
Zunes ry As K gad, f 
ð ôV, ôf _ ,. ôV, 
gz Fe à ÿz TEE- T J ôz 
et en additionnant, après regroupement des termes semblables 


; = ôV, ôV, ôv, 
div(fV}= f: (GE CPE E) + (a grad, f + V,: grad, f + V, grad, f) = 


= f: div Ÿ + Ÿ- grad f. 


[1-4] f est une fonction de point à valeurs scalaires; À et Ë sont deux champs de vecteurs; exprimer à Paide 
d'opérateurs différentiels portant sur ces quantités les trois expressions 


rot(f: 4) ; div(A À B) ; rot(4 À B). 


1° Calcul de rot (f À) : 
Nous aurons par exemple (dans un système d’axes orthonormé) : 


j- A E _ 24, 
= (grad, f ' A, — grad, f : A) + f : rot, À 
= (grad f À À), + frot, À 
Calculs analogues pour les autres composantes : on peut finalement écrire : 
rot (fA) = gad f À A + f rot À. 
2 Calcul de div (4 A B) | 


Toujours dans un repère cartésien, nous écrirons 


nB = HU A Bi + EU À 


LA 


B), + 


a 


R Ô , = a 
div (A z A B), 
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ðA ôB ðA ôB 
_ y. r Z A TIR et A 
ôx Beti ôx ð Bimi ôx 
ðA ðB ðA ôB 
=. B A: X x. aa z 
+ dy x T z ôy ôy z Ay Ô 
ðA ðB ðA ðB 
X. B dde reve e ` X 
i ög ? ai Ôz Ôz Er, ôz 
OA, A ; : : 
= B, T Ee + (termes analogues obtenus par permutation circulaire sur x, y, z) 
ôB ôB 
= CN à À 
af Iy z; ) + (id...) 
soit : 
div (£ À B) = B,' rot, À — A, rot, B + (termes analogues) 
et finalement : 


div(A À B) = B roA — A rot B 
3° Calcul de rot (4 A B) : 


La composante suivant Ox du champ de vecteurs cherché vaut : 


ð 
ar pie A B), = gy (48, — À, B,) - gz (4.8, z A,B,) 
ðA ôB ðA ôB 
e PE à Mecs x z 
dE HU 7 SE 7 
A6. E 04 


(sn de 2-0) + 4 (++): 04, _ , ÔB, 


ôx ðy z * (ôx dy z x 8x * ôx 
ðB ôB ôB ðA, ÔAÀ, ôA ôB ðA 
pe X x x : x \ _ ` X Fete z y xX . x 
(atA de eA ze) i pita T ôx 


Les deux termes ajoutés à la première ligne pour en rendre la forme plus symétrique sont précisément 
ceux qu’il convient de retrancher à la seconde pour obtenir un résultat analogue. Nous trouvons finalement : 


rot, (4 À B) = B- grad A, — À: grad B, + A, div Ë — B,: div À 
(on aurait de même : 


rot, (À À B) = B- grad A, — £- grad B, + A, div B — B, : div 4, etc...) 


[1-51 Démontrer les expressions (1-30) et (1-31) du Laplacien en coordonnées cylindriques et sphériques. 


1° Coordonnées cylindriques : 


Partons de la définition (f étant un champ scalaire) : 
Af = div (grad f) 
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avec ici 


- 1.0/ -L 
gad, f =G o gad f Say . grad, f = 
et 
ste: 1 ð 107, ôy, 
AFS À DT * à 
soit z : 
1 0 1 ð 
af = (ol) + ie f 
p ôp ĉp p? ôð ôz? 
2° Coordonnées sphériques : 
Le principe du calcul est le même avec maintenant : 
a _Lôf nf 
grad, f =. : gad fsp gg : grad, f = 30 
et 
arad pk LE L 2% 
r? ôr r 00 r‘tg0 r'sinð ôĉọ 
D'où 
1 ô of 1 f&?f  cosô ôf 1 of 
De (r 1) F2 e sin 60) r°sin° ooy 


as =(P) + ! Le + A 


r? : sin 00 r?. sin? 0 ĉo’ 


[1-61 Démontrer la relation (1-33) | 
rot (rot V) = grad (div Ñ) — AŸ. 


Étudions, dans un repère orthonormé, la composante suivant Ox : 
— ð z ô 4 

rot, (rot V) = — — — (rot, V 

à ob P) = 3 (rot, P) — (rot, P) 


L 2 (Va _ Wa) 2 (Va _ Wa) _ 2 (Wa 2v) (BVa , PYA, 
ôy\ ôx dy ôz\ ôz  ôx ôx\ ôy  ôz dy? 07? 
On peut ajouter au premier terme, retrancher au second : 


2 (VÀ _ V, 
ôx \ ôx ôx? 


Ô ,. » 
3x (div Y) —- AV, 


d’où 
rot, ÿ = 
qui est bien la projection de (1-33) sur laxe Ox. 
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EXERCICES Rappels mathématiques - 1 (SUITE) 


[1-7] Soit le champ de vecteurs : 
> xj — pi 
46,» = IX 
x+y 
dans le plan xOy où les deux axes cartésiens ont pour vecteurs unitaires i etj. 


a) Calculer rot À a 


b) évaluer intégrale | À: dM le long d’une courbe C fermée. (Il est recommandé d'utiliser les 
(e) 

coordonnées polaires (p, 0) et de distinguer le cas où l’origine © est entourée par C de celui où elle 

est à l’extérieur du domaine limité par cette courbe. 


c) interpréter le résultat de b) grâce à celui de a). 


a) Avec 
= y. ; xX 


x2 + y? x2 + y? 


il vient, pour la seule composante non nulle de rot À : 


orin M NE j- G f 2y? ) 
4 ôx êy x+y (x? 2») RS (xH yr 


soit 

a a E 

ty) Hy 
En dehors de l’origine, où ni le champ À, ni son rotationnel ne sont définis, nous avons : rot À = 0. 

d L'élément différentiel dM intervenant dans la circulation s’écrit 
dM = (dx}i + (dy)j 
d'où Nr 
A: dM = (x dy — y dx}: (x? + y?! 

En coordonnées polaires. il vient : 


X = pcosû  y=psin® x? + y? = p? 


dx = — psin 0: dû + cos 0: dp dy = p cos 0 : dû + sin 0 'dp 
x dy — y dx = p*(cos? 0 + sin? 0) d0 + p dp(cos 8 sin 0 — sin 0 cos 6) 


et la circulation cherchée s'écrit finalement : 


Come d db, 1 j 
M © ( FA 6-0 


La figure jointe montre que, pour une courbe 
qui entoure le point ©, 0 va croître d’une valeur 
initiale (par exemple 0 = 0, en A) et augmenter Fic. 1.10 
de + 2x au bout d’un tour complet (le signe 
dépendant du sens de parcours); au contraire, si # est à l’extérieur de C, 0 oscillera entre deux valeurs extrémales 0o 
et 0,, mais sa variation totale sur un tour complet sera nulle. 


c) D’après le théorème du rotationnel, la circulation calculée en b) est égale au flux, à travers le domaine 


limité par la courbe C sur le plan xOy, du rotationnel de A, Ce rotationnel est nul presque partout (exception 
faite d’une indétermination à l’origine). Si C n’entoure pas O, le domaine d'intégration pour le calcul du flux, 
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appels mathématiques - 1 (SUITE) 


est tel que rot À y soit partout nul : de même le flux, et la circulation de À; c'est bien ce que nous avons trouvé 
plus haut. Si C entoure O, nous devons tenir compte uniquement du voisinage immédiat de l’origine, où la 


discontinuité de Å et de son rotationnel peut jouer un rôle important. Nous savons en tout cas que le résultat 
est indépendant de la forme de C : c’est bien ce que donne un calcul direct. 


(Cf. ex. 2-3 pour une étude plus détaillé de la singularité du champ À) 


[1-8] Soit, dans un système cartésien orthonormé, les deux champs de vecteurs : 4[4, = — 
A,=0}et A'(A, = — B: y, 4, = 4, = 0). Calculer rot À et rot À’ et montrer que Č 


d’un potentiel scalaire dont on donnera la valeur. 


On a : 
rot, A = 0 ; rot, A =0 ; rot, À = B 
rot, d'=0 ; rot, À = 0 ; rot, A = B 
D'où 
rot C = Tot (A — 4') = rot À — rot À = 0 
et 
3 ôf By Of _ Bx of 
C = — grad f avec To era er 0 
; B 
d'où l'on tire : f = — 5 (xy) + (constante). 


; > 1 
En coordonnées sphériques, V n’a de composante non nulle que Ÿ, = ~. (1-19) donne alors : 
r 
; id 1 d 
div Ÿ = == (PY) ==> (1) = 0 
à 1 r? a 


La divergence est nulle. 7 

(Le calcul exclut qu’on se place à l’origine, où ni V ni V, ne sont définis). 

(Donnons, bien qu'il soit plus long, le calcul portant sur des composantes cartésiennes dans un système 
orthonormé. On a alors : 


x EA sr 
Mes ei Vi 
y, ôV. ôv, 
z= X + pn- 
ME a E T 


(en mettant à profit (1-6). D'où : 
3x2 + y? +779 3 3r° 


3 
r r° r? r° 


div Ÿ 


expression nulle elle aussi, sauf à l’origine où elle a une forme indéterminée.) 
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EXERCICES Rappels mathématiques - 1 (SUITE) 


> F 


[1-10] Faire de même un calcul direct de rot V avec Ÿ = —, 


Nous aurons par exemple, en coordonnées cartésiennes : 


Ve LP Rte SO = 2: _3\./» — } Le 
dy êz ð? êz \r? rj Nr r/ 
Ici encore, le problème reste posé du voisinage de l’origine, où Ÿ n’est plus défini. (On peut remarquer toutefois 
que si nous caleulons la circulation du champ coulombien Ÿ sur un petit cercle de rayon e, de centre O, tracé 
dans le plan yOz, nous obtenons un résultat nul car Ÿ, radial, est partout normal au cercle. Cette circulation est 


égale au flux de rot Y à travers le disque limité par le cercle, et vu la géométrie adoptée, ce flux fait lui-même 


intervenir rot, Ÿ au voisinage de O. La valeur nulle trouvée pour le flux montre qu’ici n ‘apparaissent pas en fait 
de singularités analogues à celle rencontrée dans l’exercice (1-7). Et nous pourrons écrire, sans restriction : 


rot V— 0.) 


Calculer sa divergence. 


Prenons k comme que unitaire dun axe Oz de coordonnées sphériques. W ma alors de composante 


sin 


non nulle que V, = et (1-19) donne : 


s 1 ô 
div V = —— Te = 0. 

rsin ĝo TA 
Ici encore, non-définition de Ÿ et de sa divergence à l’origine, mais sans difficulté particulière (le flux de Ÿ à 
travers une sphère de très petit rayon, centrée en O, est nul : il en est de même de la divergence à l’origine du 
champ étudié). 


[1-12] k est un vecteur unitaire; soit V= k À F. Calculer la circulation de Ÿ sur le cercle de centre O (origine 


du rayon vecteur F) et rayon À situé dans le plan P perpendiculaire à k en O. Calculer le flux de rot V 
à travers le disque limité par le cercle précédent, et une demi-sphère qui s'appuie sur ce cercle, Comparer 
ces différents résultats. 


Prenons Oz suivant #. V est partout tangent au cercle, et de module constant R. La circulation cherchée 
est RTR) = 27R°. 
Va pour composantes : V, = — y, V, = x, V, = 0. et rot, = 0, rot, V0, rot, V — 2. 
rot V est donc uniforme, parallèle à Oz et de module 2, son flux à er un cercle perpendiculaire à Oz, de 
surface nR? est 2rR? : on retrouve la circulation de FV, ce qui « vérifie » le théorème du rotationnel. 
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athématiques - 1 (SUITE) 


EXERCICES Rap: 


(Le flux $ de rot Ÿ à travers une demi-sphère (Z) de rayon R et d’axe Oz se calcule par l'intégrale : 


ni2 
= f (272R? sin 0 d0) - [rot | - cos 6. 
0 


(la parenthèse représente l’aire découpée sur la sphère entre les cônes d’axe Oz et de demi-angles au sommet 
respectifs 0 et (0 + d8); cos 8 est le cosinus de langle entre la normale à la sphère et rot Ñ, qui est parallèle à Oz. 
Quant à rot Ñ}, on sait qu’il vaut 2). Il reste : 
n/2 1 
P = (2n R?) 2 sin 0 cos 0 d0 = 2r R? [ 2x dx = 2rR?. 
0 JO 


C’est une nouvelle vérification du théorème du rotationnel.) 
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2] RAPPELS MATHÉMATIQUES - 2 


Nous regrouperons ici, sans lien logique, quelques considérations sur des raisonnements ou des méthodes 
de calcul que nous aurons à mettre en œuvre fréquemment dans la suite. 


I. UTILISATION DES NOMBRES COMPLEXES EN PHYSIQUE. 


Beaucoup d’attention est apportée, en Physique, aux grandeurs variant sinusoïdalement : nous consi- 
dérerons ici une variation sinusoïdale avec le temps 1; mais on traiterait sur le même modèle des variations 
par rapport à d’autres paramètres : coordonnées d'espace par exemple. 

Soit par exemple x(t) et y(t) deux telles grandeurs, caractérisées par une même pulsation œ (+); elles sont 
de la forme : 


x = x COS (at + p) y= yo Cos (ot + p). (Electromagn. 2-1) 


Xo Et Yo sont les amplitudes de x et y; (at + p) et (œt + Y) sont les phases de ces deux quantités. 

Comme l’a déjà montré le traitement des courants en régime sinusoïdal, il est possible, pour une vaste 
catégorie de problèmes, de remplacer les grandeurs réelles telles que (2-1) par des grandeurs complexes, que 
nous continuerons d’ailleurs à désigner par x, y, etc. pour ne pas alourdir les notations (il est évident que les 
grandeurs en question sont traitées, suivant le cas, comme réelles ou comme complexes). 

Le passage des grandeurs réelles aux grandeurs complexes se fait suivant les règles suivantes : 


1° La grandeur complexe a pour module l'amplitude de la grandeur réelle à variation sinusoïdale : 
2° La grandeur complexe a pour argument la phase de la grandeur réelle : 

3° La grandeur réelle est ainsi égale à la partie réelle de la grandeur complexe. 

A (2-1) seront ainsi associées les grandeurs complexes : 


x = Xo' Exp Ü@t + p) y= yo’ Exp wt + y) (Électromagn. 2-2) 


Le passage de (2-1) à (2-2) peut se justifier en remarquant que l’information contenue dans les notations 
réelles — variation sinusoïdale dans le temps, pulsation de valeur œ, amplitude et phase à l’origine (pour 
t = 0) données — se retrouve intégralement dans les notations complexes et réciproquement. 

D'autre part, si cela devient nécessaire, on peut toujours, à l'issue d’un calcul portant sur les grandeurs 
complexes, retrouver immédiatement les grandeurs réelles associées simplement en prenant la partie réelle 
des valeurs complexes. 

Mais l’utilisation des « complexes » présente de nombreux avantages : 

1° Comme un même calcul fait intervenir des grandeurs sinusoïdales de même fréquence w, toutes les 
grandeurs complexes qui y interviennent contiennent un même facteur exp iwt. On peut donc simplifier par 
ce facteur, ce qui allège d’autant les calculs. 

2° Entre grandeurs complexes apparaissent des relations de proportionnalité (avec des facteurs eux- 
mêmes complexes) qui sont dissimulées dans les calculs utilisant des lignes trigonométriques. Ainsi 
xX = Xo COS œt et y = y, Sin œt ne peuvent être décrites comme proportionnelles (car les fonctions cosinus 
et sinus ne le sont pas); au contraire les grandeurs complexes associées s'écrivent : x = x, exp (iœt) 


y = Yo EXP (or — 5 et l'on a 


y=- 29 x (Électromagn. 2-3) 


Xo 


Une relation telle que (2-3) condense sous une forme particulièrement commode deux renseignements : l’un 
sur le rapport des amplitudes de x et y, l’autre sur leur différence de phase. 

Cette possibilité du « calcul en complexes » a été abondamment exploitée dans l'étude des courants 
sinusoïdaux : elle est à l’origine de la notion d’impédance complexe. 

3° Les opérations de dérivation (ou d'intégration) des grandeurs complexes se traduisent par une simple 
multiplication (par iœ ou son inverse) : le passage d’un cosinus à un sinus est, ici encore, pris en compte par 


(1) Par abus de langage, on appellera parfois œ la « fréquence » de la grandeur considérée. On sait en fait que la fréquence y 
vaut œ/27. 
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H. 
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le facteur + i. Il en résulte que des calculs en fait analytiques (résolution d’une équation différentielle par 
exemple) peuvent se traiter comme des calculs algébriques, beaucoup plus simples. 
Terminons par une double mise en garde contre les pièges que recèle une utilisation hors de propos du 


calcul « en complexes » : 

1° Celui-ci doit être limité à l'évaluation d'expressions du premier degré par rapport aux grandeurs 
sinusoïdales. 

2° Les grandeurs sinusoïdales ainsi représentées doivent toutes être de même fréquence, cette fréquence 
étant fixe et indépendante du temps. 

Ainsi, lorsqu'il devient nécessaire de calculer le produit de deux grandeurs variant sinusoïdalement 
dans le temps, il est indispensable de revenir à un calcul réel : le produit de x par y, (2-1), vaut 


XoŸo COS (@t + p): cos (at + y) = Í xoyolcos (2œt + p + y) + cos (p — W)] 

ce qui fait intervenir un terme constant par rapport au temps, et un terme de pulsation 2w. Le produit des 

grandeurs complexes associées est en xoy, Exp i(2œt + @ + W) : seul le terme de pulsation 2% est présent, 

mais avec une amplitude double de celle calculée plus haut. Le terme constant n'apparaît pas. Grandeurs 

réelles et complexes cessent d’être simplement reliées, et, bien qu’on puisse rétablir la situation par des arti- 

fices, l’essentiel des avantages du calcul par les complexes disparaît, et il vaut mieux renoncer à cette méthode. 
De la même manière la somme de cos w,f et cos w,t vaut : 


w, + © W, — © 
COS œt + COS œ@,t = 2 cos (ee t) cos (A r) 


La somme des grandeurs complexes associées est 


exp i(œ@,t) + exp i(œ,t) = exp (ne rex (ue r) + exp (es | 


= 2exp (Pai): oos (217221), 


Ici encore, il ny a plus de relation simple entre les grandeurs réelles et complexes; il est préférable de se can- 
tonner à un calcul sur les premières. 


UTILISATION DES SYMÉTRIES EN PHYSIQUE. 


On invoque fréquemment des « raisons de symétrie » pour simplifier un certain nombre de problèmes. 
Tentons de préciser un peu cette notion. 

Supposons qu’un certain système physique Z crée un certain « effet » &. (Nous entendons ici « effet » 
dans un sens très général : toute manifestation de la présence de à quoi peut s'intéresser le physicien.) 
Ainsi peut être un ensemble de charges ou de courants, & le champ électrique ou magnétique créés par 
ces charges ou ces courants, etc... 

Soumettons F à une transformation géométrique T (translation, rotation, symétrie par rapport à un 
point ou à un plan), & subit la même transformation. Si nous notons TẸ le transformé de S, TF aura pour 
effet TS, le transformé de &. Ceci revient à dire qu’une distribution de charges électriques emporte avec elle 
son champ électrique, les surfaces équipotentielles associées etc...; de même pour une distribution de cou- 
rants et le champ magnétique associé. Nous exprimons ici, finalement, le caractère homogène et isotrope de 
l’espace où # crée &. 

Il peut se faire que, pour une certaine transformation géométrique T, TI coïncide avec S : S est « inva- 
riant » dans T. Plus généralement, il peut se faire que, pour certains systèmes 4 et certaines transformations T, 
on puisse passer de TF à F par une «transformation physique » (et non plus géométrique) P, dont l’effet 
sur & soit connu (& devenant PE) : par exemple, P correspondra au changement du signe des charges élec- 
triques, au changement du sens des courants, au « renversement du temps » (toutes les particules se meuvent 


HMI. 
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«à l'envers » sur leurs trajectoires, vitesses et moments cinétiques sont changés en leurs opposés), etc... Z est 
alors invariant dans le produit PT, qui transforme & en PTE. Comme ces deux effets sont, en dernière analyse, 
l’un et l’autre produits par le même système 2, ils doivent eux-mêmes coïncider. D’où la condition qui résume 
cette discussion : 


PEFS] 


l f (Électromagn. 2-4) 
PTE = 6 


Illustrons sur un exemple : soit à trouver le champ électrique E créé par une répartition uniforme de 
charges électriques sur un axe 4. Cette répartition est invariante dans toute rotation autour de 4; par suite 


du champ Ë en un point M on peut déduire par rotation autour de 4 le champ en 
tout point du cercle d’axe À passant par M (en particulier les modules des champs 
sont les mêmes). La répartition est invariante dans toute translation parallèle à À : 


on a donc même champ Ẹ en tous les points de toute droite parallèle à 4. Le plan 
(M, 4), où M est un point quelconque, est un plan de symétrie du système : la 
symétrie par rapport à ce plan laisse le système invariant; elle doit laisser le champ 
E(M) (où M est lui-même respecté dans la symétrie) inchangé. Ceci n’est possible 
que si E(M) est dans le plan de symétrie. Le plan H, mené de M perpendiculaire- 
ment à 4, est lui-même un plan de symétrie pour le système, donc pour le champ. ÆE(M) doit être 
dans H. Finalement le champ est porté par laxe intersection des deux plans de symétrie (i. e. : la perpen- 
diculaire à À menée par M). 

Que deviendrait ce raisonnement si À portait un courant constant, dont on cherche l'induction magné- 


tique Ë? Les résultats tirés de l'intervention de rotations et de rotations peuvent 
être repris sans changement. Mais quand on fait intervenir des symétries-plans, il 


convient de tenir compte de la nature de pseudo-vecteur de B : 
— dans la symétrie par rapport au plan (M, 4) qui respecte le système, 
B(M) doit être invariant. Or ce champ est transformé en l'opposé de son symé- 


trique par rapport au plan. La coïncidence n’est ici possible que si B est perpen- T 


FIG. 2.1 


diculaire au plan de symétrie; 

— la symétrie dans J change le sens du courant; pour retrouver le système 
initial il faut faire intervenir une opération P = inversion du courant. P change 
aussi le sens de Ë. Finalement B(M) doit coïncider avec « l'opposé (à cause de P) de l'opposé » (car B est 
un pseudo-vecteur) de son symétrique dans 11, donc avec son symétrique dans I. B(M) doit être dans II, ce 
qui confirme le résultat précédemment obtenu. 

Sur ces deux exemples très simples, l’analyse symbolisée par l'implication (2-4) est, évidemment, un peu 
lourde. Nous avons cependant eu ainsi l’occasion de voir comment il convient de la mettre en œuvre pour tirer 
de la symétrie des systèmes étudiés le maximum de renseignements. 


LES « DISTRIBUTIONS ». LA DISTRIBUTION DE DIRAC. 


L’électromagnétisme fera intervenir, de manière générale, des distributions continues de charge électrique, qu’on peut 
caractériser par une densité volumique, p, fonction de point à valeur scalaire : si dr est le volume d’un domaine infinitésimal 
de l’espace entourant un point M, ce domaine contient la charge : 


dq = p(M) : dr (Électromagn. 2-5) 


Cependant, il reste souvent nécessaire de considérer des charges ponctuelles, ne serait-ce que du fait de l'existence de 
particules élémentaires, portant une charge finie (au sens de non infiniment petite) alors que leurs dimensions restent 
inappréciables, même à l'échelle atomique. 

A priori, il semble difficile de concilier les formalismes respectivement applicables aux charges continûment réparties 
et aux charges ponctuelles. Nous verrons ici comment l'invention d’un nouvel être mathématique, appelé « distribution », 


29 


RAPPELS MATHÉMATIQUES - 2 (SUITE) 


30 


permet de résoudre la difficulté. L'utilité de ces distributions (+) dans des démonstrations ultérieures justifie cette introduction 
préliminaire au cours d’Electromagnétisme. 


1° La distribution de Dirac ô à une dimension : 


Commençons, pour simplifier, par considérer des distributions de charge sur un axe x’x (et non en volume). De telles 
distributions, si elles sont continues, peuvent être caractérisées par une densité linéaire A(x) telle que l'élément de x'x 
compris entre les abscisses x et x + dx porte la charge : dq = A(x) : dx. 

Inversement, on peut s'intéresser au cas d’une charge ponctuelle, par exemple de valeur unité, et placée à l’origine ©. 
Cette charge peut être considérée comme la forme limite d’une suite de distributions continues, 
mais qui seraient de plus en plus concentrées autour de l’origine, les densités À,(x) correspondantes À Al) 


devant : 
r us i À3 
a) ne prendre de valeurs appréciablement différentes de zéro que 
dans un voisinage de plus en plus étroit de O quand n croît A2 
b) prendre, dans ce voisinage, quand n croît, des valeurs croissantes ` 1 
de manière à (Electromagn. 2-6) J= A7 


+o 
c) garder un « poids » total égal à l'unité : f A(x) dx = 1 
LS = 0 14 x 


Par exemple, on peut considérer que la charge ponctuelle est limite, pour n tendant vers 
linfini, des distributions continues définies par les densités À,(x) « en créneau » : FIG. 2.3 


(Électromagn. 2-7) 


On voit donc, comme d'ailleurs il fallait s’y attendre, que la « densité linéaire » qu'on pourrait associer à une charge 
ponctuelle en O est une « fonction » hautement singulière puisqu'elle est : 


— nulle en dehors de O 

— «infinie » en O 

— et telle que son intégrale sur tout intervalle contenant O soit finie et égale à 1. (L'intégrale sur un intervalle ne 
comprenant pas l’origine étant évidemment nulle.) 


Cette densité ne peut plus être définie au sens des fonctions, mais on peut profiter du fait que son intégrale garde un 
sens pour convenir de la définition suivante : 


on appellera « distribution de Dirac », ou « distribution 8(x) » (?) une « fonction » : 


— nulle pour x #0 (Électromagn. 2-8) 
+w 

— vera | &(x) dx = 1 (Électromagn. 2-9) 

— limite d'une suite de fonctions satisfaisant à (2-6), par exemple les À,(x) (2-7) (Électromagn. 2-10) 


2° Quelques propriétés de la distribution 6 : 


Soit (x) une fonction régulière près de l’origine, où elle admet un développement de Taylor : 


(1) On ne confondra pas la « distribution » au sens mathématique, telle qu’elle sera introduite rapidement ici, et la distribution 
de charge électrique à quoi nous avons fait plus haut référence. 
(2) Par abus de langage, on dit encore « fonction de Dirac » ou fonction 6. 
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Les intégrales 
+ © 
n= eapex 
(où 1, est définie par (2-7)) peuvent donc se développer en 


+ oo 


1,= 00 | usé + #0 | 


Le] oo 


j 0 +o 
xå (x) dx + ro [ x? x) dx... 


La première des intégrales (condition (2-6)) vaut 1. Les autres tendent toutes vers zéro, quand n tend vers l'infini : elles sont en 
effet, en module, majorées par : 


+in re i 
[ bel A(x) dx < — [ À,(x) dx = — 
-im PI ipm 


+i/n l +i/n l 
[ x?l ; Ax) dx < Lf A(x) dx =— etc. 
in N J in n? 


A la limite, 7, tend donc vers @ (0). Or cette limite, puisque ô(x) est la limite des À,(x) peut s’interpréter comme l’inté- 


+o 
grale | (x): 6(x) : dx. D'où la formule fondamentale : 


7 


l AT dx = @(0) (Électromagn. 2-11) 


Très succinctement, nous résumerons en disant que la distribution de Dirac est une fonction nulle en dehors de l’origine, 
singulière à l’origine, mais sur laquelle nous pouvons continuer raisonnements et calculs à condition que cette fonction 
figure dans un intégrande, et qu’on tienne alors compte de la relation (2-11). 


Remarques. 


1° (2-9) n’est qu’un cas particulier de (2-11), obtenu pour ọ(x) = 1. 
2° Sila charge ponctuelle placée à l’origine a la valeur q (et non 1), elle correspond à la densité g. (x). On aura ainsi : 


NT a[ aprés 


3° Sila charge (disons de valeur 1) est à l’abscisse xo, il sera naturel de lui attribuer la densité 8(x-— xo), soit à(x) 
avec X = x — x,. Cette densité est en effet, comme il convient, 
— nulle pour X # 0 (soit pour x # Xo). 
+o +o 
-— de poids égal à l’unité | 8(x — xo) dx = [ ë&(X) dx = 1 


=o ~o 


Avec cette distribution de Dirac décalée, (2-11) devient : 


+o +o 


| plx)ô(x — xo) dx = [ p(xo + X)o(X) dX 


= om T ao 


= [p(X + xo)lx=0 
Soit 


+o 
| g(x)ô(x — xo) dx = (xo). (Électromagn. 2-12) 
-%0 
Il convient de souligner que la notion de distribution est beaucoup plus générale que ne pourraient le laisser supposer 
ces quelques notes, limitées au cas de la distribution de Dirac, particulièrement précieuse au physicien, et, surtout, que 


cette introduction a volontairement fait appel aux aspects les plus intuitifs et a laissé de côté toute exigence de rigueur. 
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RAPPELS 


Avant d'étendre à trois dimensions la notion de distribution, il est peut-être bon de suggérer, par un exemple simple, 
l'intérêt des développements qu’elle permet. 
Soit la fonction « échelon » : 


pour x > 0. 


Elle est discontinue, donc non dérivable à l’origine. Partout ailleurs sa dérivée Y'(x) est nulle. 
La fonction Y peut être considérée comme limite des fonctions 
continues Ÿ(x) : À 


Celles-ci sont dérivables (quoique non continûment dérivables) et leurs 33 12 
dérivées sont les 4,(x) dont la limite est ô (x). 

Il est donc naturel de dire que Y(x) n’est pas dérivable à l’origine au sens 
des fonctions, mais l’est « au sens des distributions », et que l’on a : 


Y'(x) = ë(x) 


FIG. 2.4 


(On peut dire encore que Y(x) est primitive de ë(x); on doit avoir : 


xX 
Y(X) = | 8(x) dx. (Électromagn. 2-13) 
— œ@ 
Effectivement, cette intégrale est nulle si l'intervalle d'intégration ne contient pas l’origine (d’où Y = 0 pour X < 0) et 
vaut 1 autrement (d’où Y = 1 pour X > 0).) 
Une relation telle que (2-13) peut être utilisée à condition que ô ou Y’ figurent dans une intégrale. Soit Ax) une 
fonction exempte de singularités et s’annulant pour x = + co. On peut chercher à calculer : 


+o 
) f(x): Y'(x) dx (Électromagn. 2-14) 
D’après (2-13) et (2-11), on doit trouver /(0). Si l’on en avait le droit, il serait tentant d'intégrer (2-14) par parties : 


+ © + + 
í fer-ds=l rY] _ f'(x): Y(x) dx. 
Le terme tout intégré est nul; l’intégrale vaut : 
+ 00 Le] 
= rral- 1] = f(0) 
-0 0 
C’est bien le résultat donné par (2-13). Autrement dit, l'introduction des distributions permet d'étendre certains procédés de 


calcul hors de leurs limites « normales » de validité. 
Nous nous en tiendrons ici à ces quelques indications. 


3 La distribution de Dirac à trois dimensions : 


L'extension des notions précédentes à trois dimensions ne présente pas de difficultés. Comme nous l’avons vu ini- 
tiaiement, il s’agit de définir un être mathématique susceptible de représenter la densité (ici volumique) de charge qu’on peut 
associer à une charge ponctuelle placée à l’origine. Commençons par une charge unité. La « densité » cherchée sera la 
distribution de Dirac 8(M) (à trois dimensions) qui satisfait aux conditions suivantes : 


— ô(M) est nulle pour tout point M distinct de l’origine ©. 


( Électromagn 2-15) (— on peut considérer (1) comme limite, pour n tendant vers linfini) de la suite des fonctions 


1 
1 4n 1 
P,(M) nulles en dehors de la sphère de rayon — et prenant à l’intérieur de celle-ci la valeur (=) ) 
n : 3n 
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— elle satisfait pour toute fonction © (M) régulière en O à la relation : 


[ | f S(M): g(M) : dt = (0) (Électromagn. 2-16) 
et en particulier à : 


il ô(M)dr = 1 (Électromagn. 2-17) 


Autrement dit, on peut traiter (M) comme une « fonction » nulle en dehors de l’origine O, singulière en O, mais sur 
laquelle on peut calculer à condition qu’elle figure systématiquement comme facteur dans Pintégrande d’une intégrale 
triple, et qu’on respecte (2-16) (£). 

Pour une charge q placée en O, la densité deviendrait q - (M). 

Pour une charge (unité) placée en un point P distinct de O, la densité sera la fonction de Dirac « déplacée » de O en P : 
nous la noterons 5 (M), ou (M, P); (dans le dernier cas, la « variable » est M; P joue le rôle d’un paramètre). 


4 La divergence du champ coulombien. Le laplacien de 1/r : 


Les relations (2.15) et (2.16) (ou 2.17) caractérisent la distribution de Dirac. Si une quantité satisfait ces relations, 
on peut l'égaler à 8(M). 
Considérons le champ coulombien centré à l’origine 


E(M) = = F-0M , r=0M (Électromagn. 2-18) 
r3 


Nous avons montré que sa divergence est partout nulle là où le champ lui-même est défini, c’est-à-dire en tout point M 
distinct de O. (ef. Exercice 1.9). div Ë satisfait à (2-15). Calculons f | f div Ë dr à l'intérieur de la sphère £ de rayon R 


centrée en O. D’après le théorème de la divergence, cette intégrale est égale au flux de Ē à travers X, soit, Ë étant radial 


et constant en module sur Z : 
a pe 2 1 
div E : di = 4r Ri :—— = 4n. 
R? 
(2) 


1 E 
T div É satisfait ainsi à (2-17). D'où l'identification : 


div = d4n-ù (Électromagn. 2-19) 
r 


On sait que le champ de Coulomb (2-18) dérive du potentiel coulombien en 1/r : 


à —— 
Peine grad E) 
r° r 
ES 1 
On a donc: div E = — a) = OM) et: 
1 z 
af) = — Arû (Électromagn. 2-20) 


(1) Dans (2-16) et (2-17) le domaine d'intégration est un volume quelconque comprenant O. De (2-17) on tire que ô est 
homogène à l'inverse d'un volume. 
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i 
Considérons maintenant la fonction (de M) : PM (où P est un point fixe quelconque). Son laplacien (calculé par 


rapport à la « variable » M, ce que nous préciserons en le notant À PM n'est autre que (2-20), mais décalé de O en P. 
D'où la formule fondamentale : 


1 2 
adsa) = — 4n: ôp(M) (Electromagn. 2-21) 


qu’on peut interpréter en disant que, g(M) étant une fonction régulière au voisinage de P quelconque, on a : 


(Dans 2-22, la variable d'intégration est M, point moyen de l'élément différentiel d'intégration que nous avons noté à 
cause de cela dr,,. Le domaine d'intégration est nimporte quel volume contenant P.) 
Remarquons pour terminer que nous démontrerons plus loin : 


— au chapitre 3, l'égalité : 


(Électromagn. 2-22) 


> P . MP : 
divy ==; = divp FE (Électromagn. 3-12) 
— dans le cours de l’exercice (3-4), l'égalité 
1 1 R 
zr] = Aio Elect . 3-24 
a ri a M z) (Electromagn. 3-24) 


(qui n’est d’ailleurs qu’une autre forme de (3-12). On a donc, aussi bien que (2-22) : 


(Il Aa (pen du = — 47: g(P) (Électromagn. 2-23) 


EXERCICES 2] Rappels mathématiques - 2 


see an 


[2-1] On considère un solénoïde, de section circulaire, et infini. Quelle orientation imposent les conditions de 
symétrie au champ magnétique Æ qu’il crée? 


On voit immédiatement que le solénoïde est de révolution autour de son axe x’x, et qu'il possède une 
symétrie de translation parallèle à cet axe. Il en résulte que Ë sera le même en 
tous les points d’une même parallèle à l’axe et que, si M et M’ sont deux points 
situés sur un même cercle axé sur x'x, on passe de BM) à B (M) par la rotation, 
autour de x'x, qui amène M sur M". 

De même, M étant un point quelconque de l’espace, le plan P mené de M 


perpendiculairement à x'x est un plan de symétrie pour le solénoïde. B(M) doit £ 
donc être invariant dans la symétrie par rapport à P : orce pseudo- -vecteur est 


transformé en l’opposé B" de son symétrique À” par rapport à P. B et B" doivent 
coïncider et À doit être perpendiculaire à P, et parallèle à x’x. Nous voyons que la 
symétrie nous impose ici des conditions très strictes : il ne reste plus, finalement, 
qu'à calculer une seule composante (parallèle à x’x) de B(M), fonction d’une 
seule variable (la distance p de M à xx). 


[2-21 Appliquer les raisonnements du (2-1) au potentiel vecteur À du solénoïde infini. 


On peut reprendre les mêmes arguments en ce qui concerne les rotations autour de x’x et les rotations 
parallèles à cet axe. 

Dans la symétrie autour de P, qui laisse le solénoïde invariant, AM) doit lui-même être invariant : 
il s’agit d’un vrai vecteur, transformé en son symétrique par rapport à P. Il n’est invariant que s’il est He 
dans le plan de symétrie. 

Dans le plan Q défini par M et laxe x’x du solénoïde, celui-ci est géométriquement invariant, mais le 


sens du courant qui le traverse doit être renversé, de même que celui de A(M). Ainsi, dans la symétrie par 


rapport à Q, suivie d’un changement de sens du courant, A(M) doit être conservé: par suite, ce vecteur doit coïncider 
avec l’opposé de son symétrique par rapport à Q :ilest tangent au cercle d’axe x’x passant par M (en coordonnées 
cylindriques d’axe x'x, À n'aura de composante non nulle que 4,). 

(On pourra montrer que la prise en considération de la symétrie par rapport à Q confirme les résultats de 
l’'Exercice 2-1.) 


[2-31 Reprendre lexercice (1-7) pour un champ de composantes 


Re TRE = 2 


ox + y? + 82 x? + y? + e7 


Montrer que rot, À, dans la limite & — 0, tend vers la distribution de Dirac à deux dimensions 8(M4) 
(à un facteur près qu’on déterminera) et achever d'interpréter les résultats de (1-7). 


Nous avons ici 
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soit 
2e? 
(x? + y? + 82) 


rot, À = 


Cette fonction a un maximum — en O (x = y = 0) et décroît rapidement à partir de l’origine pour 


s’annuler à l'infini : elle a donc Palus de fonction «en pic» requise d'une fonction qui doit tendre vers la 
distribution de Dirac à ou vers une distribution multiple de ô. 


Le « poids » total de rot, À (flux de rot À à travers le plan xOy) est commodément calculé en coordonnées 


polaires; il vaut : 
co 2 o -1 > 
[ 2xp dp ee 2xe? Í a e (=) 
o (p? + &} o (X Her X + 87/6 


i 1 ; à 
soit : 2ne? = 27, valeur constante, qui reste donc finie quand & — 0, 
E ; > ; 
Nous pouvons donc écrire, pour le champ de vecteur À de l'exercice (1-7) : 
rot, À = 216(M) 


ce qui explique que le flux de rot À à travers un domaine quelconque du plan xOy contenant l’origine 
(limité par une courbe fermée entourant ce point) soit égal à 2x, et que le flux à travers un domaine ne contenant pas 
le point O soit nul, comme l'avait montré un calcul direct. 


[2-4] Démontrer (2-19) en imaginant un champ de vecteurs, fonction d’un paramètre, qui tende vers le champ 
de Coulomb quand ce paramètre tend vers une certaine valeur limite, et dont la divergence soit une fonction 
«en pic», ayant son maximum à l’origine O, décroissant de plus en plus rapidement à partir de O quand 
le paramètre tend vers sa limite, et gardant cependant un « poids » total (intégrale sur tout l’espace) 
constant et égal à 4r. 


On cherchera par exemple un champ À radial à symétrie sphérique [en coordonnées sphériques, une 
; 3 i Ne 1 d 
seule composante 4, non nulle fonction d’un seul paramètre r; la divergence d’après (1-19) vaut ENT (Ah 


i i ; Jr” dr 4 
dont la divergence soit constante sur toute la boule centrée sur 0 de rayon e (e tendant vers zéro) et nulle’ à 
l'extérieur. La valeur k de cette divergence dans la boule, si son « poids » total est 4x, est telle que : 


3 
us ie Due 
3 e 
Dans la boule, on aura donc : 
«3 
ne k? ; PA, =k +e ; Re 
dr 3 3 r? 


On prendra désormais nulle la constante d'intégration c pour que le champ À reste partout défini. 
Pour r > 6, il vient au contraire : 


d + 
L(a) =0 4 =$. 
F r 
Imposons au champ À d’être continu pour r = e. Il vient : 
re soit d=k>=1 
3 e 3 
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. > à 1, ri i 
Finalement, le champ A a pour composante radiale A, = 7 à l'extérieur de la sphère de rayon e; autrement 


A 
dit il se confond avec le champ coulombien. Quand e tend vers zéro, la zone où il se distingue de ce dernier, 
c'est-à-dire la boule r < e devient de plus en plus petite : on peut donc dire que À tend vers le champ de 
Coulomb pour & — 0; dans le même temps div À tend vers 4. 


[2-5] Quel sens peut-on donner aux distributions 8(kx) et ô(kx — a) où k et a sont des constantes? 


Dans la relation (2-11) qui définit la distribution de Dirac quelle que soit la fonction œ@(x) : 


x est une variable d'intégration, une variable muette. Dans ô(kx), posons kx = X, d’où x = T Une intégrale 


telle que : | p(x)ô(kx) dx devient alors a (F) -ô(X)} dX soit d’après (2-11) : 


| p(x) (kx) dx = Z (0) et 8(kx) = = (x) 


Nous écrirons de même : ô(kx — a) = LE — i) = l g(x — A (a d’après (2-12) nous voyons que 


k 
Fitégrae | p(x)ô(kx — a) dx donne z3) ) 
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RAPPELS ET COMPLÉMENTS 
D'ÉLECTROSTATIQUE 


Revenons rapidement sur les résultats essentiels d'Électrostatique. On considère ici des charges au repos, 
c’est-à-dire des charges immobiles les unes par rapport aux autres, et qu’on prend dans un repère qui leur est 
lié. Les seuls changements de repère qui puissent intervenir ici correspondent donc à un changement d’orien- 
tation des axes de coordonnées ; il ne peut y avoir passage dun premier repère à un second qui serait 
en mouvement par rapport au premier: dans l’un au moins des repères les charges ne seraient plus au repos. 


LOI DE COULOMB. CHAMP ÉLECTRIQUE. 


Deux charges (immobiles, mais on cessera désormais de le préciser), q, et q, interagissent par des forces 


« de Coulomb », soit de q, vers q; : E 
F7, pa l da'a Th A 
Areo ri (Électromagn. 3-1) 
avec £o & (36710°)~t M 


he Lo 

FÈ = M M „ est le vecteur qui joint la charge (ici q, ) qui exerce la force au point 1 2 
où se trouve la charge (ici g,) qui la subit. FiS: 3:1 
On peut décomposer (3-1) en introduisant la notion de champ électrostatique (ou champ électrique) : 
une charge crée autour d’elle un champ électrique Ë, et si une autre charge est placée dans ce champ elle subit 


une force qui lui est proportionnelle. 


— q en P crée en M, avec F = PM, le champ y 
q4 FP ` M M) 
Ē(M) = — (Electromagn. 3-2) 
Ames r° que 
— Q en M subit une force D 
F = Q:E(M) (Électromagn. 3-3) Fic. 3.2 


On admettra la superposition linéaire des champs : si plusieurs charges q; sont présentes en des points P;, 
elles créent en M un champ qui est la superposition, c’est-à-dire la somme vectorielle, des champs Æ (Â4) que 
créeraient en M les charges q; si elles agissaient seules : 


+ 


(M) = DE(M) ; E 


LA 


(M) = di an ; F,= PM (Électromagn. 3-4) 
és ri 


Si la distribution de charge créant le champ est une distribution continue définie par la densité volu- 


mique p(P) non nulle sur un volume Q, on aura de même au « point d'observation » M le champ, obtenu par 
superposition : 


m P)- PM : 
E(M) = TAERA dtp (Electromagn. 3-5) 
(o 418o * PM? 


où nous notons dr} le volume élémentaire d'intégration pour bien marquer qu’il entoure le « point source » P 
où se trouve la charge responsable du champ. (Ce point source, dans l'intégration, décrit l’ensemble du 


volume Q où sont présentes les charges.) 


Remarque. (3-4) peut se ramener à (3-5) si nous utilisons les propriétés de la distribution ô définie par (2-15) et (2-16). 
À la charge q; au point P, correspond la densité q; : 6, (P). A l'ensemble des charges q; correspond la densité : 


p(P) = È qiôpi(P) (Électromagn. 3-6) 


RAPPELS ET COMPLÉMENTS 
D'ÉLECTROSTATIQUE (SUITE) 


(3-5) fera intervenir les distributions à dans une intégrale, où on peut les utiliser en respectant (2-16). (3-5) devient alors : 


> . PM 
Em) = Ji 1 
i Ar£ PM? 
c'est-à-dire (3-4). 
Disposant dun champ de vecteurs, ici le champ électrique, nous devons calculer son rotationnel (ou 


vérifier qu’il dérive d’un potentiel scalaire, ce qui entraînera l'identité à zéro dudit rotationnel) et sa divergence. 
C'est l’objet des deux paragraphes suivants. 


II. POTENTIEL ÉLECTROSTATIQUE (POTENTIEL ÉLECTRIQUE). 


Un calcul direct montre que le champ (3-2) d’une charge ponctuelle dérive du potentiel scalaire : 


VI) = 24.4 = PM) 


= Ares F 
E(M) = — grady V(M) (Électromagn. 3-7) 
{rot E = 0) 


V n’est, en fait, défini qu’à une constante près : l’expression (3-7) assure son annulation à linfini. 


(Prenons P, point fixe, comme origine des coordonnées. La composante suivant l'axe des x de — grad,, V est en : 
q ô /1 q 1 à H 
PRE EER el DER PEN — ; $ 
Anteo Ôxy Y Ame, r2) Xy 


(on utilise la formule de dérivation des fonctions de fonctions), d’où finalement (avec (1-6) : —— Ne E (M). 


Comme les champs, les potentiels peuvent se « superposer » (c'est-à-dire s’additionner); un ensemble 
de charges ponctuelles q; en P; donne donc en M : 


y — qi LA 
— peaa a 
(M) Z; Ineo P,M (Électromagn. 3 8) 


Une distribution continue définie par la densité p(P) sur un domaine Q donne en M le potentiel : 


N 1 p(P)-dr, i | 
V(M) = [Il fe PM (Électromagn. 3-9) 


Remarque. On peut repasser de (3-9) à (3-8) en prenant la densité (3-6) et en utilisant, de nouveau, la propriété fonda- 
mentale (2-16) des distributions de Dirac. 


Notons enfin que (3-8) et (3-9) définissent des fonctions V(M) continues : une discontinuité du potentiel 
correspondrait à un champ électrique infini, ce qui est exclu tant qu’interviennent des charges ou des densités 
de charge finies. 
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CALCUL DE div£ THÉORÈME DE GAUSS. 


1° Calcul du flux de Æ. Théorème de Gauss : 


(nous ne faisons que rappeler cette propriété fondamentale sur laquelle insistent systématiquement tous 
les cours d’Électrostatique). 

On démontre par des arguments géométriques que le flux sortant $ du champ (3-2) créé par une charge 
ponctuelle g à travers une surface fermée S vaut : 


a)  — 0 si q est à l'extérieur de S 


b) & — si q est à l’intérieur de S. 


Autrement dit, ® vaut (1/8) x (charge totale à l’intérieur de S}. Par suite si le champ Ë est créé par 
une distribution continue de charge (champ (3-5), on aura, Y étant le volume intérieur à S : 


ð= [fz -dŠ = 5 [fo d'y 


En transformant ® par le théorème de la divergence (1-16), il vient : 


il divy Ë: dty = (íf N: dre (Électromagn. 3-10) 
y y Êo 


dans (3-10), le volume d'intégration (#7), de même que la « surface de Gauss » S, est quelconque. L'identité 
des deux intégrales (portant toutes les deux sur ce même volume Ÿ” d'intégration) implique donc l’identité 
des intégrandes; d’où la relation : 


div Ë = = (Électromagn. 3-11) 


qui constitue la « forme locale » du théorème de Gauss. (Forme « locale » parce que (3-11) relie deux grandeurs 
prises au même point, et en un seul point. Au contraire, la « forme intégrale » (3-10) fait intervenir les valeurs 
du champ en tous les points de S, ou celle de sa divergence, et de la densité de charge, en tous les points de Ÿ”.) 


Nous avons donc pu, par wn biais, évaluer la divergence du champ électrique. Mais nous sommes armés pour en 
donner un calcul direct : ce sera un bon exercice de préparation aux calculs de l’Électromagnétisme. 
2° Calcul direct de div Æ : 


Partons de (3-5), le champ est calculé en M; c’est par rapport à ce point que nous prenons sa divergence : 
à P): PM 
divy ECM) = divy Ill POSE, à 
a ATE PM? 


Le calcul de la divergence fait intervenir des dérivations par rapport aux coordonnées de M. Celles-ci commutent avec 
l'intégration sur (Q) par rapport à P. On peut donc calculer en fait : 


divy EM) = (ff es divy pru | . a 
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Ps P): PM r s 
L’intégrande fait intervenir divy f Cana } soit, d’après (1-20) : 
PM 


PM PM — 
(P “diva = | +=: grady (p{P)) 
PF) | PM: e 


z5 ôp(P 
Le dernier terme est nul, car p(P) ne dépend nullement de M et grad,, (PCP) est nul (ses composantes : p(P) par exemple, 


x 
étant nulles). Il reste donc une intégrale où P est le point variable d'intégration et dans l’intégrande de”laquelle figure 
une divergence calculée par rapport à un autre point : M. On profiterait beaucoup plus facilement des propriétés de la 
divergence si la divergence figurant dans l’intégrande était calculée par rapport au point d'intégration P. (Rappelons que 
cette condition est évidemment nécessaire à l’application du « théorème de la divergence » lui-même.) D’où notre but : 


$ PM Š Le 4 b y, i R 
transformer div = en un div, (U) où U est un champ de vecteurs à déterminer. En fait, on trouve facilement qu’on 


di PM = di MP Él 3-12 
iva FR = divp TP) (Electromagn. 3-12) 


peut écrire, tout simplement : 


On peut démontrer (3-12) directement, par un calcul analogue à celui qui nous a donné (1-34), 
On peut aussi utiliser (1-34) précisément, et (1-20) : 


… {PM Fos ss à 
divu (E) = div (PT) + Ph grady (=) 
PM? PM? M3 


Or 


diva (PM) = = God + 2 Ou) + = Ex) = 3 


(l’origine étant prise en P, et cette valeur est aussi celle de div, (MP) (avec maintenant une origine en M). Par ailleurs, on a : 
EE 1 — 1 
grad =) = — grad ( SAREAT } 
“APM PNPM? 
PM 1 Sn one 
divy (©) = - dive (MP) + MP : gradp ee ) 
P 


PM? M? 


E 
= divp|—— |}. 
MP? 


D'où finalement : 


Il vient donc : 


à 1 MP 
divy E(M) = ji dtp’ divp ( 
4Teo JJJ, MP’ 


) y figure dans une intégrale); il vient : 


) * p(P) (Électromagn. 3-13) 


M 
Utilisons alors (2-19), ((3-13) s’y prête puisque ai E 


MP: 


divy E(M) = l [ | Jar PP 


Ares 


-> ff a HP): B4P) 


1 
soit d’après (2-16) : p(M) 


Nous retrouvons bien (3-11). 
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ÉQUATIONS DE POISSON ET DE LAPLACE. 
De (3-7), (3-11) et de la définition du laplacien, nous tirons immédiatement l'équation de Poisson 


AV = — 2 (Électromagn. 3-14) 
(e) 


Si l’on se trouve dans une région de l’espace où il n’y a pas de charges (les charges créant le champ et le 
potentiel étant ailleurs) (3-14) se simplifie; on a, avec p = 0, l'équation de Laplace 


AV = 0, (Électromagn. 3-15) 


Ces équations « aux dérivées partielles », comme les équations différentielles pour les fonctions d’une 
variable, ne suffisent pas à déterminer la fonction (ici V) sur laquelle elles portent : cette indétermination 
correspond d’ailleurs au fait qu’à un champ électrique donné correspond non pas un potentiel unique, mais 
une famille de potentiels, ainsi qu’on l’a déjà noté (chap. 1, § V, 3°). 

On achève de déterminer le potentiel en imposant des conditions aux limites : en général, on choisit de 
travailler sur un potentiel qui s’annule à l'infini : 


V(r => co) - 0 (Électromagn. 3-16) 


C'est, tacitement, cette hypothèse qui a été adoptée dans la formule (3-9). Notons cependant que (3-16) 
n’est conséquence de (3-9) que si les charges ne s'étendent pas à linfini : c’est bien le cas d'intérêt physique 
pratique, mais on exclut alors les idéalisations classiques que sont le fil infini chargé, ou le plan infini chargé. 
Dans de tels cas, des conventions particulières, et qu’on précisera soigneusement, devront être adoptées. 

Remarquons en outre que (3-14) et (3-15) sont des équations linéaires, ce qui est une conséquence du 
principe de superposition (ou une manière de l’exprimer) : si deux densités de charges p, et p}, créent, respec- 
tivement, les potentiels V, et F,, la densité p =å p, + up, crée le potentiel V = 47 + uV (å et u étant deux 
scalaires quelconques.) Cette propriété a, en particulier, été très souvent mise en application enélectrostatique 
des conducteurs. 


ÉNERGIE ÉLECTROSTATIQUE. 


1° Charge ponctuelle soumise au champ d’autres charges : 


Soit un ensemble de charges ponctuelles q4, g, … 4n, au repos en des points M4, M, … M, Amenons au 
point M une autre charge ponctuelle q. Comment définir, et calculer, son énergie poten- M; 


tielle électrostatique dans le champ des charges (immobiles) q;? 8G M 
Supposons que q soit amenée en M depuis l’infini (où elle n’est plus soumise à lac- P à 2 
tion des charges q;), et très lentement (de manière à rester, en tout instant, pratiquement M Go 
dans un état électrostatique). q Mp 
En chaque point P de sa trajectoire, la charge q est soumise au champ E(P) créé en P s dn 
par les g.. Elle subit donc une force F = q - E(P), et, pour que son déplacement soit effec- FIG. 3.3 


tivement très lent (que la situation des charges soit quasi-statique), il est nécessaire qu’une force extérieure F, 
non électrique, s'oppose à l’action de F. 
Le travail de cette force, F = Ë, sur l’ensemble du déplacement de q entre linfini et M, est : 
M M a M 
w=g F -dM = -$ F-dM = -a E- dM. 
co 


Le co 
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Si V est le potentiel dont dérive le champ Ë (Ë = — grad V), ce potentiel s’annulant à linfini, nous 
trouvons finalement : 


W = a[V(M) — V(o)] = q: V(M) 


Ce travail fourni au système (les charges q; et la charge q) par l'extérieur, emmagasiné par le système, et 
que l'extérieur pourrait récupérer si l’on procédait à la transformation inverse, constitue par définition l'énergie 
potentielle du système (par rapport à un seuil d'énergie correspondant à la charge q à l'infini et échappant à 
l’action des q) : 


E,=q'V (Électromagn. 3-17) 


Cette énergie potentielle résulte visiblement des interactions électrostatiques entre charges; on peut 
donc la nommer : énergie potentielle électrostatique de la charge g. 


Remarques importantes. 


a) Il est essentiel de noter que l’on calcule ici l’énergie potentielle liée à la mobilité éventuelle de la seule 
charge q. Les charges q; restent fixes; pour les amener en place, il a fallu fournir une certaine énergie, mais 
qui n’est pas prise en compte par ce calcul. En d’autres termes, on s'intéresse ici à l’énergie dépendant des 
seuls degrés de liberté du système qui correspondent aux mouvements de la charge q. 


b) Le potentiel V dans (3-17) est celui que créent les autres charges (les g;) au niveau de la charge q. Ce n’est pas le 
potentiel en M (point occupé par q) quand les q; et q ont atteint leur position finale, ce qui signifie en particulier que q est 
en M : d’ailleurs, un tel potentiel cesse d’être défini sur une charge ponctuelle telle que q, puisqu'il varie en 1/r, où r est 
la distance à cette charge. Le calcul a été mené en ne considérant que l’action des autres charges sur une charge donnée : 
comme c'était d’ailleurs — implicitement — admis depuis le début de ce chapitre, on considère qu’une charge magit pas 
sur elle-même. Il n’y a donc pas lieu de tenir compte de la « self-énergie » de la charge, c’est-à-dire de l'énergie correspondant 

.à une éventuelle interaction de la charge avec elle-même. (On verra plus loin la nécessité d'introduire le concept de self- 
énergie.) 


2° Énergie potentielle d’un ensemble de charges ponctuelles : 


Soit n charges q … qp en des points M, … M. Il est naturel de définir l'énergie potentielle électrosta- 
tique de l’ensemble de cette distribution comme le travail qu’a dû fournir l’extérieur, contre les forces élec- 
trostatiques qu’elles échangent, quand on a amené l’ensemble de ces charges, initialement à l'infini, et infiniment 
distantes les unes des autres (donc sans actions mutuelles), jusqu’à leurs positions finales A. 

Pour ce faire, on peut imaginer qu’on a 


— déplacé d’abord q, pour l’amener de linfini en M,, les autres charges 
restant à l'infini. D’après (3-17), il a fallu dépenser pour cela une énergie & Go 
WM=g n M, Mo 
où V est le potentiel noté en M, à la fin du déplacement, soit zéro, d’après (3-8), 
puisque les q;(i # 1) sont restées infiniment éloignées de M; M, 
— amené ensuite q, de l'infini en M,. Le potentiel que subit q, en M, est dû 
à la seule charge à distance finie : q, ; il vaut : Qa 


qda 1 ` 
V, = —  (oùr,, = M,M;) 
2 ATE Fiz ta 3 Ž FiG. 3.4 


L'énergie fournie par l'extérieur dans ce déplacement est W, = q, V3 
— amené ensuite q, en M, où elle subit l’action de q, et q,, qui créent le potentiel : 


= al + _42. t 
3 Anéoris 4MEo F33 
(avec r,; = M M3, F33 = M M3). L'énergie fournie par l'extérieur est ici W = q; © V4. 
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— continué des déplacements analogues jusqu’à la n°" charge. 
Au total, l'extérieur a fourni le travail : 


W =W +W, +W +a +W, 


Ce travail constitue lénergie potentielle électrostatique de l’ensemble des n charges, soit : 


n 1 ; z 
E,= È af 2) (r; = MM} (Électromagn. 3-18) 


On rendra cette formule plus symétrique en supprimant la restriction į > j et en sommant sur les i et les j 


de 1 à n en évitant seulement de les égaler : le terme en Ski qui n'apparaît qu’une fois dans (3-18) se manifeste 


; ; je ii à 3 ij o, ; 5 ass $ 
alors deux fois. L'expression symétrisée doit donc être affectée d’un facteur qui traduiten fait l'intervention, 
dans ce problème, d'interactions entre paires de particules : 


1 qq; 1 : 
= — l z2 — .° n - 
E =3 L'an, 5 L qi: V(M)) (Électromagn. 3-19) 


Dans (3-19), V(M,) est le potentiel créé, au niveau de q; par toutes les autres charges. On notera l'apparition du 


facteur 1/2 qui distingue essentiellement (3-19) de (3-17). 
(Dans (3-19) comme dans (3-17) les self-énergies des charges ne sont pas prises en compte.) 


3° Extension à une distribution continue de charge : 


Supposons maintenant les charges réparties continûment dans un volume Q, avec une densité p(M). 
L'élément de volume dt entourant un point M où règne le potentiel V( M), contient une charge p(M) dt qu'on 
peut assimiler à l’une des charges q; de (3-19) ; cet élément de volume contribue à énergie potentielle par : 

1 


dE, = 5 (M): di): V(M) (Électromagn. 3-20) 


et, au total, énergie potentielle électrostatique de l’ensemble de la distribution est 


E, = alll p' Vidt (Électromagn. 3.21) 
Q 


(Remarquons qu'ici notre calcul n'exclut pas la « self-énergie », où ce qwelle devient pour une distribution continue 
de charge : le V(M) qui figure dans (3-20) est le potentiel dû à l’ensemble de la distribution, y compris la charge contenue dans 


dr. 
Mais supposons que dr soit une sphère infiniment petite de rayon e centrée en M; elle porte une charge dq œ p(M) : &, 


qui crée un potentiel en = œ p(M) : #7. La contribution à V(M) de la charge présente dans dr tend vers zéro avec cet 
e 

élément de volume. En fin de compte, la « self-énergie » disparaît pour une distribution continue de charge électrique : 

c’est une notion qui n’a d'intérêt que pour des charges ponctuelles.) 

4 Densité d'énergie électrostatique : 


De (3-11} et (3-21) nous tirons : 


Ste dv Ë 
E, afff div E dt 


Or (cf. (1-20)) : 
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E, = $ fff av (VE) dr + [fl aE dr 


Remarquons d’autre part que le domaine d’intégration peut, sans inconvénient, être étendu du volume Q où 
sont présentes effectivement les charges à l’espace entier : ceci revient à prolonger p par une fonction nulle à 
l'extérieur de Q. 
La première intégrale obtenue plus haut devient jil div (VE) dz soit : 
esp 


lim | | | div (VË) dr 
R= © 


où l'intégration porte ici sur le volume intérieur à une sphère X de rayon R tendant vers l'infini et contenant Q. 
Le théorème de la divergence permet d’ailleurs de transformer en : 


lim [l VE : dÿ 
R>o JJ 


(flux du vecteur VË à travers X). Or pour des R très grands, V est, au plus, en 1/R, LÀ en 1/R2, et le flux de VE 
à travers X (desurfacetotale œc R?)en 1/R. Ce flux tend donc vers zéro quand R tend vers linfini. Au total ilreste 


i - (Électromagn. 3-22) 
(esp} 


Tandis que, d’après (3-21), l’existence d’une énergie électrostatique semblait intimement lié à la présence 
de charges électriques, nous obtenons une expression — équivalente dans le cas de l’électrostatique — qui 
suggère au contraire que l’essentiel est la présence d’un champ électrique, et que, là où un tel champ existe, 
il apparaît une densité d'énergie électrostatique 


et (avec (3-7)) E, devient : 


d 2 5 
Le = fo | (EÉlectromagn. 3-23) 


Quel est le point de vue le plus fécond ? 

On verra plus loin que, pour les régimes variables, c’est (3-23) qui se généralise le plus directement : en 
fait, on peut concevoir des situations où subsiste un champ électrique alors que les charges qui l’ont créé ont 
disparu. (3-21) donnerait alors un résultat nul, alors qu’il y a encore énergie électromagnétique, en accord 
avec (3-22) ou (3-23). 
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[3-1] Démontrer la relation (3-19) en imaginant que la répartition finale de charges ponctuelles qu’on considère 
(a, en Mis eg; en M) est constituée en amenant (très lentement) mais en même temps toutes les 
charges depuis l’infini jusqu’à leur position finale de manière qu’à chaque instant f, la figure formée par 
les positions instantanées M, (t), M (8 … M,(9 … M (£) des charges gj, 4, oo Gj e 4, Soit homothétique 
de la figure formée par les positions finales par rapport à un point fixe O. 


A l'instant ż, la force qui s'exerce sur la charge q; et que doit contrebalancer la force extérieure F;(r) vaut : 


=. F (t) = di. aMDM{) 
es fi MOMAY 


On posera OM (O = AD: OM » À() variant de linfini à 1 au cours 
du processus. Dans le temps dr, À subit l’accroissement di, et M(t) vient 
en M(t + dé) par un déplacement de vecteur : 


dM, = OM, di. 


MES M. =. 90) MM, ÖM, di 
Ameo jazi A1) M M? 


et le travail sur l’ensemble des charges que fournit l’extérieur vaut : 


1 MM, OM;\ di 
dW = aw, = (- 1 qa) it 
> | mE Ga) M M? | 2 
if) 


La somme sur i et j j peut être « symétrisée » : on y verra intervenir le couple GP par MM. OM, puis par 
M M ;: Où, = — MM M; OM. Au total, il apparaît M, 1 M, - (OM, — Où J= MM .M?, associé à chaque couple 
Ge D (i À D, ce qu'on peut représenter par une somme sur i et 7 indépendamment (l'égalité restant exclue) 
à condition de diviser par 2 l’expression obtenue. Au RE le travail élémentaire fourni par l'extérieur vaut 


dW = — 1 y 4. a då 
Ares 24 MM, 92 
et le travail total au cours du processus (À variant de q à 1) vaut : 
1 
w= 1 .l.s 4 | M 1. s 1. 44 
4neo 2 Fi MM; Ja 4 2 {Fj ATE MM, 


L'énergie potentielle £, du système de charges est opposé de ce travail, soit : 


1 1 
S a n 
ý er a A 


on retrouve bien (3-19). 


[3-2] On considère la sphère de centre O, rayon a, portant la densité volumique uniforme de charge p. (La 
charge totale de la sphère est g.) Il wy a pas d’autres charges que celles présentes dans la sphère : 


1° Calculer le champ électrique Ë et le potentiel V en tout point de l’espace, intérieur et extérieur à la 
sphère. 
2° Calculer l’énergie potentielle #7 de la sphère chargée en sommant sur la sphère Pexpression à Vo àr 
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(où dr est l’élément de volume). 
3° Retrouver W, en sommant sur l’espace la densité d'énergie électrostatique exprimée en fonction de Æ. 


1° Le champ est partout radial (il doit être invariant, au point M, dans toute rotation autour de l’axe OM, 
qui laisse invariante la sphère chargée). Il ne dépend, d’autre part, que de la distance de M à O, soit r. En 
coordonnées sphériques, Ë wa donc qu'une composante non nulle : E,(r). Le flux de ce champ à travers la 
sphère de centre © et de rayon r, soit (r), est donc simplement égal au produit de E (r) par la surface de cette 
sphère, soit 
D(r) = 4nr?- Er). 


D’après le théorème de Gauss, (r) vaut 1 Q(r) où Q(r) est la charge totale présente à l’intérieur de la sphère 
0 


de rayon r, soit : 
r>a Q(r) = (=F) 
Arr r? 
r<a: r) = = q'— 
Q(r) A re 
Il vient donc : 
r>a: Er) = boga 
Antep 7’ (0) 
r<a Ej{r) = ataa 
âne, aè 


(Le même résultat pouvait être obtenu à partir de la forme locale du théorème de Gauss (3-11); le calcul, 
en fait, a déjà été présenté à l’occasion de l’exercice (2-4).) 


ôV 
En intégrant (1), eten lui imposant d’être nul à linfini, nous obtenons le potentiel V(r) G que Er) = — F) ; 
; -1.4 
r>a:  Vfr)= dre + 
1 2 
r<a: Vir)= e a 
AT, aè 2 


La constante d'intégration c est choisie pour assurer la continuité de V en r = a; d’où: 


__qg B-r 
ue A 24 } 


W, = Il T p: V(r)dr 
(sphère r < a) 


La contribution de l'élément de volume compris entre les sphères r et (r + dr) (égal à Axr° dr) est 


dW, = 3 Grr dr): p- V(r), d'où 


2° Soit à calculer : 
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et en réexprimant en fonction de g et a : 
EUR 
W, = 
£ 207eça 
2 


te 5 E&E 
3° La densité d'énergie est en 2> soit : 


E r? 
— pourr <a: F (ze) Home dont l'intégrale sur le volume intérieur à la sphère vaut : 
af 


4TEo 
E E a Y 1 a 
“0. > 60. PP tu 
5 E) + 1 anr° dr x r 5 G) 2e 4n z 
— porr >a: Žo, (2 sk dont l'intégrale sur l’extérieur de la sphère vaut : 
2 \4re rê 


to (a Ÿ. ("am dr eo_a Y gnd 
2 4T£o a r° 2 ATEo a 
La somme de ces deux expressions est : 
2 2 
nes 2a | 5 207604 


On retrouve bien l’expression obtenue plus haut. 


[3-3] Une distribution à symétrie sphérique autour d’un point © crée un potentiel électrostatique 
q e” (> 0) 
Are, F q > 0) 
en un point situé à la distance r de ©. 
1° Calculer la densité de charge p(r) et préciser son signe. 
2° Par application directe du théorème de Gauss, trouver la charge Q(r) à l’intérieur d’une sphère de 
centre © et rayon r. Montrer que la comparaison des résultats des 1° et 2° exige qu’on suppose placée en © 
une charge ponctuelle dont on précisera la valeur qo. 


1° (3-14) nous permet d'écrire immédiatement pr) = — £ AV avec (grâce à l’expression (1-31) du 
laplacien en coordonnées sphériques). 


= ripe- q aleia) 
r? dr dr amer? dr r rP 
q4 d -Ar —àr — Àr 
== — je 1 + Ar = — — + Àr)e + àe A 
4negr? dr LA 1 re n ) 


Au total, nous obtenons 


Il s’agit donc d’une distribution de charge négative. 
2° Le champ électrique n’a de composante que radiale : 
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EXERCICES | 3 Rappels et compléments 
d'électrostatique (SUITE) 


Le flux de ce champ à travers la sphère de rayon r (le champ étant partout normal à la sphère et constant 
en module) vaut simplement : 


et la charge Q(r), d'après le théorème de Gauss, vaut : 
QC) = q(t + Arje-# 
charge positive, qui ne peut être produite par la seule densité (négative) de charge p(r). 


Or celle-ci introduit dans la sphère de rayon r la charge : 


La 


O= [am ar: dr = ga [rec ar 
0 


(e) 
2 I -À 1 SE 
fem) [re LR +) 6" dr 
(e) 


ree7# 1 e z 
- w3 f- s ara = -a + al +ine 


Q(r) diffère donc de Q’(r) par la constante q, quel que soit r. Ceci signifie bien qu’à la charge en volume décrite 
par p(r), il convient d’ajouter une charge ponctuelle, placée à l’origine, et de valeur qo = q(> 0). 

(Cette charge ne pouvait apparaître dans le calcul du 1°, limité au domaine de l’espace où V(r) est défini, 
ce qui exclut l’origine. La méthode du 1° ne nous fournira d'indications sur la présence en © de la charge q 
que si nous l’étendons à l’origine en utilisant des distributions. 


On y parviendra en écrivant V(r) = 2A + V'(r) avec 
ANEo F 

; -4(1- e7”) EE E 

V'(r) Ar. > (non singulier à l’origine). 
On a : 
— air — Àr 
AY! 2 q A d f (1 e PR p | 
Ane? dr y7 r 
St L —-1+{(+are "} = — E A + Àr)e”#] 

amer” dr 4re,r° dr 


(la constante — 1 disparaît dans la dérivation). Pour r # 0, AV” se confond ainsi avec AV, et la densité p(r) 


1_ où Pon reconnaît le potentiel 
AREF q 


d’une charge ponctuelle q placée en 0. Ce terme se manifesterait par une contribution en — = (M) à AV 
0 


crée en fait la contribution V'à V seulement. V comporte en plus un terme 


et en qô(M) à p(r). Cest ce terme que le calcul du 1° ne prenait pas en compte.) 


[3-4] Démontrer par un calcul direct que, quelle que soit la fonction d’une variable f@), on a : 


Auf (PM) = Apf (MP) 


Appliquer à la démonstration de (3-24). 
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HüLIN. - Électromagnétisme ~ 1 


EXERCICES 


ppels et compli 
d'électrostatique (SUITE) 


On a par exemple, M étant pris comme origine momentanée et, x», Yp, Zp désignant alors les coordonnées de P, 
r la distance PM, en appliquant systématiquement les formules (1-6) : 


afir) ft) usa 


ÔXp r y 
& 2 1 | 1 
i 1? =P Ea + f'E) xe’ z 
P 
O om atya 3 S Rt yetz 
ApJ =f E > Pg | 3 = 
= f + : E 


Cette expression reste symétrique par rapport à M et P, puisqw’elle ne dépend, comme f(r) elle-même, 
que de la distance entre ces deux points. On aurait obtenu le même résultat en dérivant par rapport à M à 


P fixé : d’où la relation demandée. En spécialisant au cas où f(r) = + nous obtenons : 


a u(7r) = asd Fr) (Électromagn. 3-24) 


relation qu’on trouvera de nombreuses occasions d'utiliser en Electromagnétisme. 


[3-5] Démontrer par un calcul direct qu’un potentiel V défini à partir d’une densité de charge p par la relation : 


=} p(P) dtp 
VM) = res ji PM m 
satisfait à équation de Poisson A,,V = — 202 (3-14) 
0 


Il vient 


IEOR 


(on peut permuter les dérivations par rapport aux coordonnées de M que suppose le calcul du laplacien et 


p(P) 


l'intégration par rapport à P. Dans les dérivations par rapport à M qu'implique le calcul de aG ; 


p(P) qui ne dépend pas de M ni de ses coordonnées peut être considéré comme une constante et mis en 


facteur. Il reste : 
1 1 


ire | | | p(P) (7) TA 


soit, d’après (3-24) : 
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EXERCICES 3 Rappels et compléments 
d'électrostatique (SUITE) 


et,en utilisant (2-21) ou (2-22) : 
ue C.Q.F.D. 


Au = AE £o 


[3-6] Deux charges ponctuelles opposées (+ q et — q) sont immobilisées à une distance a Pune de Pautre. 
Calculer leur énergie électrostatique par (3-19), puis (3-22). Les résultats sont-ils compatibles? Comment 
peut-on interpréter leurs divergences éventuelles? 


(3-19) nous donne 
2 


P? Anega’ 


cette expression est négative, comme il est normal : les charges, de signesopposés, s’attirent; elles fournissent 
du travail à l’extérieur quand on les laisse s’approcher jusqu’à une distance a l’une de l’autre; le système a 
cédé du travail pour arriver dans sa position définitive : l'énergie emmagasinée (énergie potentielle) est négative. 

Le calcul de E, par (3-22) conduit nécessairement à un résultat différent : il faut sommer sur l’espace une 


.…, [RE 5 > i r s y 
quantité ( 5 ) positive et énergie calculée ne peut être, elle-même, que positive. 
2 


De plus, E varie en r~ ^ au voisinage de chacune des charges 
à E, du voisinage de chaque charge est- en 


E&E 


2 


2 -4 


est donc en r * et la contribution 


lim a da a lim 
e+ 0 ; e=>0 £ 
E, est positive, mais infinie (ou, si l’on préfère, égale à une intégrale divergente). 

La difficulté provient du fait que, comme on l’a noté dans le texte, le calcul par (3-19) donne une expression 
de l’énergie potentielle d’interaction qui exclut les self-énergies des charges (si elles sont ponctuelles). Au 
contraire (3-22) a été démontré à partir d’une expression (3-21) valable pour des distributions continues de 
charge, pour lesquelles le problème des self-énergies, ainsi qu’on l’a noté, ne se pose pas. L'application de (3-22) 
à des charges ponctuelles fait réapparaître les self-énergies de celles-ci, qui sont positives, divergentes, et 
masquent l'énergie d'interaction. 

On peut d’ailleurs remarquer que le champ total Ë est superposition des champs E €t Ë, créés indépen- 
damment par les deux charges. La densité d’énergie électrostatique peut se décomposer en 

EE? £o 


2 = (Ë F E,) = Ei + < E2 sa 80Ë1 E, 


Les deux premiers termes représentent les énergies créées respectivement par + q et — q isolées et sans inter- 
action; dans l'intégrale sur l’espace, ces termes conduisent à des divergences; ce sont eux qui représentent 
les self-énergies des particules considérées. 

Seul le dernier terme décrit l’interaction entre les deux charges, associée à l'énergie prise en compte par (3-19). 
On peut vérifier par un calcul direct que son intégration sur l’espace redonne 


q? 


ee ARE a 
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RAPPELS ET COMPLÉMENTS 
D'ÉLECTROCINÉTIQUE. 


On suppose maintenant que les charges électriques qu’on étudie sont mobiles : nous donnerons ici la 
définition et les principales propriétés des grandeurs qui caractérisent les courants électriques qui apparaissent 
alors. 


DENSITÉ DE COURANT. 


Soit un flux de particules chargées qu’on supposera identiques : même charge q, même vitesse Y. On en 
compte 
dN = n' dr (Électromagn. 4-1) 


dans un volume drz : ici n est une fonction, en général, des coordonnées d'espace et du temps. 
A ces charges correspond la densité volumique de charge 


p=n'q (Électromagn. 4-2) 


et un champ de vecteurs, la densité de courant : 


|= pi=n"q'À (Électromagn. 4-3) 


Si plusieurs espèces de charges sont présentes, avec des charges q, et des vitesses 7, différentes, et en nombre 
n; par unité de volume, (4-2) et (4-3) deviennent simplement : 


P = 2 n;°g; 
> P Électromagn. 4-4 
J 2 ni 'di' vi l 4 ) 


L'introduction du vecteur j défini par (4-3) ou (4-4) se justifie par son intervention dans l’ensemble des 
formules qui décrivent les propriétés des courants électriques. En voici un premier exemple. 

Soit une surface élémentaire d.S immergée dans le milieu où se déplacent les charges mobiles. On l’oriente, 
et on lui associe, selon le procédé habituel, le vecteur d5 orienté suivant la normale positive. Quelle est la charge 
d?Q qui traverse dS entre les instants f et t + dt? 

d?Q est somme des charges des d?N particules qui ont traversé dS dans (t, t + dÀ soit (pour une seule 
espèce de charges) : 


d’ = q` dN 


Les d?N particules intéressantes sont celles qui, à Pinstant « ini- 
tial » ź, se trouvent dans le cylindre de base dS, de génératrices parallèles 
à T et longues de 7 : dz : les particules se déplacent de v dt dans le temps 
considéré, et seules celles qui sont initialement à l’intérieur du cylindre ((1) 
par exemple) ont leur trajectoire qui coupe dS ((3) passera à côté de dS; 
(2) n’atteindra pas dS, etc...) Le volume du cylindre est : 


dS- (dd et N =n: dS- dt. 
D'où finalement 
d2Q = ngë: dS- dt 


soit 


= 


dQ =f- dS; dt (Électromagn. 4-5) 


C’est la combinaison j de n, q et 8 qui nous permet de calculer d2Q : premier exemple de lintérêt de cette 
grandeur. (On verrait facilement que (4-5) reste valable quand interviennent plusieurs types de charges.) 


RAPPELS ET COMPLÉMENTS 


D'ÉLECTROCINÉTIQUE. (SUITE) 


Remarque : il convient d'interpréter (4-5) algébriquement. Un d?@ > Ocorrespond soit à un flux de parti- 
cules chargées positivement qui traversent d§ dans le sens de la normale positive, soit à un flux de particules 


négatives qui « remontent » dS (Tet d5 faisant un angle obtus, avec un produit scalaire < 0). Au contraire 
d?Q < 0 correspond à q > 0, des charges qui remontent dS, ou à q < 0, des charges qui « descendent » dS. 

Soit maintenant une surface S finie (non infiniment petite), orientée, à l’intérieur d’un milieu oùsemeuvent 
des charges. Dans le temps dé, elle est traversée d’une charge dQ qui est la somme des d?Q correspondant à 
tous les éléments dS de cette surface : 


dQ = i dS : dt 
(5) 


On appellera intensité du courant à travers S la grandeur 


I= € = | | Feds (Électromagn. 4-6) 
(5) 


C’est le flux de j à travers la surface considérée, en même temps que la charge (comptée algébriquement) 
qui la traverse par unité de temps, conformément à la notion traditionnelle d'intensité. (En particulier, Pinten- 
sité du courant dans un circuit est l’intensité à travers une section de ce circuit, c’est-à-dire le flux de la densité de 
courant à travers une section du circuit.) 


RELATION DE CONTINUITÉ. 


Parmi les principes de conservation sur quoi repose la Physique (conservation de la quantité de mouvement, 
de lénergie...) figure un principe de conservation de la charge électrique. Conformément aux observations qu’on 
peut mener dans des conditions expérimentales très diverses, on admet qu’il ne peut y avoir ni apparition, 
ni disparition spontanée de charge électrique dans un système isolé; même dans le domaine des particules 
élémentaires, qui peuvent s’annihiler, on constate qu’une particule chargée ne disparaît jamais seule, mais 
simultanément avec une particule de charge opposée : le système, en ce qui concerne les deux particules en 
question, a une charge nulle après leur annihilation, mais il avait initialement une charge totale déjà nulle, 
et il y a bien conservation de la charge électrique totale du système. 

Soit alors une surface fermée S immergée dans un milieu où se meuvent des charges. 

Soit p (comme on l’a noté une fonction des coordonnées d’espace et du temps) la densité de charge. La 
charge totale Q(t) contenue dans S à l'instant ż vaut : 


Q(t) = IEG t) dtp 
y 
où y est le volume limité par S. 


p variant avec le temps, il en est de même de Q. A l'instant ż + dé, cette charge est devenue : 


Q(t + dt) = Q(t) + dQ = [re + dt) dtp 


On a donc, au premier ordre en dt : 


dQ = [fee + dt) dt» — [Lee t) dtp 


2 If [AP t + de) — pP, 8] des = Me (P, t) de: dtp. 


et par suite 
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RAPPELS ET COMPLÉMENTS 


D'ÉLECTROCINÉTIQUE. (SUITE) 
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Mais comme il ne peut y avoir, dans Ÿ; apparition ou disparition de charge électrique, si Q varie c’est que 
des charges quittent le volume #; ou y pénètrent, dans le temps dé. Le principe de conservation de la charge 
électrique totale nous permet de ne considérer, comme cause des variations dans le temps de Q, que les dépla- 
cements des particules chargées : si, au total, dans le temps dt, une charge dO’ a traversé ($) vers l’extérieur, 
il y aura un déficit de dQ’ au niveau de Q, soit : dQ = — dQ’, relation dont il est évident qu’elle est algébrique. 


Or (4-6) nous donne : ag = E d5. D'où l'égalité : 
S 


Le flux de j peut être transformé en une intégrale triple portantsur % (théorème de la divergence). D'où : 


IL + aiv) -dr = 0 


Le volume Ÿ, de même que la surface S, sont quelconques : l'intégrale ne peut s’annuler systématiquement que 
si son intégrande est nul; d’où l’équation : 


divj + e =0 (Électromagn. 4-7) 


qui porte le nom d’équation de continuité : elle traduit en effet l’absence de transformations discontinues de la 
charge électrique (par création ou annihilation). Dans le cas d’un régime stationnaire, donc indépendant du 
temps, la charge Q(f) considérée plus haut est, pour un volume #” donné, une constante : il doit entrer dans #° 
autant de charges qu'ilen sort, etjest à flux conservatif. On le voit d’ailleurs directement à partir de (4-7) : toutes 
les dérivées par rapport au temps, telles que dp/üt, s’annulent; il reste : 


divj = 0 (Électromagn. 4-8) 


EXERCICES Rappels et compléments 
d'électrocinétique 


[4-1] Un courant de symétrie sphérique est décrit par une densité de courant j qui ma de composante non nulle 
que radiale : cette composante j, n’est fonction que de la distance r au centre © de la distribution de 
courant. Peut-on ainsi réaliser un régime stationnaire? Que se passe-t-il si l’on place en © une source de 
courant, ponctuelle, qui laisse échapper, de manière isotrope, une intensité 7 constante? 


On doit, pour que le courant soit stationnaire, lui imposer de vérifier l’équation (4.8) : 
divf = 0. 
En coordonnées sphériques, on en déduit : 


i =o 


soit, k étant une constante : j, = 2 

(Le champ j est donc « coulómbien ».) Or on sait (grâce à l’Électrostatique) qu’un champ coulombien 
n'est à divergence nulle qu’en dehors de l’origine : ici, cette restriction s’interprète très simplement en 
remarquant que le courant radial permanent qu’on cherche à établir vient soit accumuler des charges à l’origine, 
soit en retirer. Il convient donc de placer à l’origine une charge ponctuelle variable qui puisse s’accroître du 
courant ou l’alimenter. Stricto sensu, on ne peut donc réaliser de régime stationnaire. 

Par contre si l’on admet, en O, la présence d’une charge ponctuelle variable, il devient possible d’établir le 
courant stationnaire « coulombien ». Le flux de ce courant à travers une sphère de centre O et de rayon r, 
c'est-à-dire l’intensité 7 à travers la sphère, vaut 

I = (4xr°) ue 


P = 4rnk 


x 


Finalement, le courant cherché est en : 


[4-2] On considère un milieu de conductivité y limité par deux plaques infiniment longues qui se coupent sur 

un axe Oz. Ces plaques font entre elles un angle œ petit; elles sont interrom- 
pues au niveau de deux cylindres, qui ont pour axe commun Oz, et pour rayons € ,-77// T j TT 7 J) Dai 
respectifs a et b. On admettra que les lignes de courant sont des circonférences A E f > 
(ou des portions de circonférences) d’axe Oz. a 
Les plaques jouent le rôle « d’électrodes » ; leur conductivité est infinie; elles 
sont maintenues aux potentiels V = 0 et V = K. 
1° Décrire le potentiel V et le champ électrique Æ dans le milieu conducteur. 
En déduire la densité de courant j. 

2° Calculer la résistance R de ce milieu entre les deux plaques pour une 
hauteur. comptée parallèlement à Oz. 

3° Calculer la puissance totale P dépensée par effet Joule dans le milieu 
conducteur, pour une portion de celui-ci limitée par deux plans perpendiculaires 
à laxe et distants de h. 

4° Calculer la puissance dépensée par effet Joule dans un élément de 
volume dr du milieu conducteur, en fonction des coordonnées qui fixent sa position et par intégration sur le 
volume du 3°, retrouver la valeur de P obtenue plus haut. 


1° Par hypothèse, il y a relation linéaire entre densité de courant et champ électrique : 


J=y È. 
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appels et compléments 
d'électrocinétique (SUITE) 


EXERCICES 


Le régime de courant étudié ici est stationnaire et j est à flux conservatif. On a ainsi 
div E = Q. 


divf = 0 et 


En utilisant des coordonnées cylindriques d’axe Oz (p, 0 et z), et en remarquant que, les lignes de courant 
étant des cercles d’axe Oz, E et j n’ont de composante non nulle que E, et jẹ, nous tirons de (1-18) les conditions : 


8=a 


E, et ją sont fonctions de p seulement. La circulation sur une telle ligne de courant, entre les deux électrodes, du 
E({p): p d9 = «' p> Elp) si p est le rayon de la ligne de courant considérée; 


champ électrique Ë vaut : f 
8=0 
quel que soit p, cette expression redonne la d.d.p. V; entre les deux plaques. D'où : 
V, yy, 
Edp) = —— ilo) = = 
o(p) jp) x- p , 


ap 
Le potentiel a même valeur en tous les points d’un plan passant par Oz et faisant langle 8 avec la plaque 


8 
Elp): p d8 = v, £. 


de potentiel nul; ce potentiel est : 
0 
| 0 


2° L'intensité 7 qui traverse le milieu peut être calculée par exemple au niveau d’une des électrodes : 


Ë et j lui sont perpendiculaires et 7, flux de j à travers la plaque, vaut simplement 
b 


p=b p=b 
I | h: dp == S Er loge 


De la loi d'Ohm V = RI nous tirons la valeur de la résistance : 


(On remarquera que R, comme il convient, diminue quand augmentent les dimensions transversales h et 
(b — a) du milieu conducteur, ainsi que la conductivité y de celui-ci.) 


3° La puissance cherchée est simplement P = — soit 


yh b 
P= V} 0e (2) 


4 Dans l'élément de volume dt = p d0 : dp : dz entourant un point de coordonnées (p, 0), est dépensée 


la puissance (cf. chapitre 6, Équation (6-11)) : 


dP=j Éd =7y"E?:dr 


soit 
2 
ap=2Vi.d 
a? p’ 
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EXERCICES | 4 | Rappels et compléments 
d'électrocinétique (SUITE) 


Sur le volume de conducteur étudié, ceci correspond à une puissance totale : 


h p=b 8=a& $ 2 à 
p= [a] ap | pa 24.1 
vo /p=a 489=0 x P 


2 b 2 
-Wina [8 2 08) 
a P x a 


Nous retrouvons ainsi le résultat du 3°. 
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RAPPELS ET COMPLÉMENTS 
DE MAGNEÉTOSTATIQUE 


i. 
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L'existence de courants stationnaires (indépendants du temps) entraîne l'apparition de champs magné- 
tiques (eux-mêmes indépendants du temps), et dont l’étude constitue la magnétostatique, dont nous rassem- 
blons ici les résultats essentiels. 


CHAMP D’INDUCTION MAGNÉTIQUE D'UN COURANT STATIONNAIRE. 


La présence, au niveau d’un point P, point moyen d’un volume élémentaire dtp, d’un courant stationnaire 
de densité j(P), entraîne l’apparition au point d'observation M d’un champ d'induc- 
tion magnétique (*) élémentaire : 


j = do ÑP) ^ PM 


dB(M) = — - dr Électromagn. 5-1) 
4x PM? i l 


(Ici les longueurs sont en mètres, j en Ampère/mètre carré; la constante Ho, per- 
méabilité magnétique du vide, vaut 
Ho = 4n x 1077 (Électromagn. 5-2) 


ce qui définit l'unité d’induction magnétique : le Tesla.) 

On admet (ce que suggère l'expérience) qu’il y a superposition des champs magnétiques créés par des 
éléments de courants distincts. S’il y a des courants stationnaires, définis par le champ de vecteur j{ P) à l’inté- 
rieur d’un volume Q, l'induction magnétique totale qu'ils créent en M sera donc : 


B(M) = Il le HORS (Électromagn. 5-2) 
a4t PM? 


Un cas particulier important est celui où les courants sont confinés à des circuits filiformes de dimensions 
transverses négligeables (négligeables devant les dimensions longitudinales des circuits, ou les distances point 
du circuit-point d'observation). La répartition du courant sur une section d’un de ces circuits perd son intérêt : 
seule subsiste la notion d'intensité Z parcourant le circuit. 

Le volume élémentaire dtp de (5-1) est en fait un petit cylindre de longueur d/, taillé dans le circuit, et 


dË(M) 


entourant le point P. Si s est laire d’une section droite du fil, on a : dtp = dl: s. . 
Le vecteur j{P) est nécessairement longitudinal au fil : les charges mobiles suivent le IP. ) 
circuit. Ce vecteur est donc perpendiculaire à la section droite s. Son flux à travers s est 
simplement JP): s, et c’est l’intensité I dans le fil. 
On voit ainsi que le module du vecteur J: dr qui intervient dans (5-1) peut se 
transformer en : 
j'd=j'sdi=1" dl 


Si au segment dl, nous associons un vecteur d/ élémentaire de longueur d/, orienté dans 
le sens dej, nous aurons, vectoriellement cette fois-ci 


1 di =F de 
(5-1) se transforme alors en : 
to 1: di A PM 


ae Électromagn. 5-3 
An PM? l i ) 


dB(M) = 


(1) Nous utiliserons, concurremment avec la dénomination légale d’induction magnétique, le terme plus simple de champ 
magnétique. En fait le champ magnétique est une grandeur distincte de B, mais qui lui est proportionnelle dans le vide. Nous n'aurons 
donc pas de confusion à craindre. 


RAPPELS ET COMPLÉMENTS 


DE MAGNÉTOSTATIQUE (SUITE) 


HE. 


et (5-2) en 
B Ho Íi d, a PM 5 
B(M) = LR Électromagn. 5-4 
on L $ 4x PM? l ) 


où C; représente le À circuit, de point courant P,, parcouru par l'intensité J.. (Remarquons qu’un circuit, 
en régime stationnaire est nécessairement fermé, d’où la notation de l’intégrale curviligne : .} (5-3) et 


(5-4) constituent la loi de Biot et Savart. Nous étendrons cette dénomination aux relations plus générales (5-1) 
et (5-2). L'induction magnétique étant définie, nous entreprenons le calcul de sa divergence, et de son rota- 
tionnel. 


CALCUL DE div Ë POTENTIEL VECTEUR. 


Calculons directement la divergence du champ élémentaire (5-1) dans un système d’axes centré en P, 
l'axe Oz étant colinéaire à j; x, y et z sont les coordonnées de M; r = PM. Il vient : 


et 


TETE ORTE 7 vs 


Pourr # 0, ona donc div (dB) = Q0. N’y a-t-il pas une singularité à l’origine (comme pour le champ électrique 
d’une charge ponctuelle), avec une div (dB) finalement proportionnelle à ôp(M)? 
Dans ce cas, l'intégrale triple de div (dÆ) sur un volume limité par une surface fermée S entourant P serait 


non nulle : or cette intégrale est égale au flux de dB à travers S. Prenons pour S une sphère de centre P ; dB est 
tangent à la sphère en tout point de sa surface et son flux à travers elle est nul. Nous ne rencontre- 


rons pas ici les anomalies vues pour le champ coulombien ; on a partout : div (GB) = 0, et de même, pour 


l'induction magnétique totale : 
(Électromagn. 5-5) 


Ainsi apparaît une propriété fondamentale : l'induction magnétique en régime stationnaire est à flux 


conservaiif. ` M 
On sait alors qu’ilexiste, pour le champ B, des potentiels vecteurs À dont il dérive, c’est-à-dire qui satisfont : 


B= À (Électromagn. 5-6) 


C’est ainsi que le dB de (5-1) dérive du potentiel vecteur 


-> _ Ho JP): de 
ee PM `’ 
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Dans le système d’axes utilisé plus haut, dÅ n’a de composante non nulle que dA, = Fe -j dr: G) du 
rot, (dA -4 Bye) otya dB, 

dy An r? r 


— ĉA Bos 1 x 
ot, (d4 = — — = {#20}; dt | —- =) = dB 
EA ôx (és e ) ý 


rot, (dA) = 0 = dB. 


et 


En sommant sur tous les points sources, nous obtenons un potentiel vecteur pour le champ À de (5-2) : 


A(M) = fl A (Électromagn. 5-7) 
Q 


Comme nous l'avons noté, la condition (5-6) ne suffit pas à déterminer À pour À donné : le passage d’un 


potentiel vecteur à un autre porte le nom traditionnel de transformation de jauge. La jauge (5-7) bénéficie des 
propriétés suivantes : 


1° ÄM) s’annule quand M tend vers l'infini, au moins dans les cas (seuls réalisés de fait) où le volume Q 
est borné. 


2° La divergence de À est elle-même nulle : 


(Électromagn. 5-8) 


Le potentiel vecteur (5-7) porte le nom de jauge de Coulomb. 
Démontrons (5-8) : 


Oublions la constante m Il faut calculer : 


vu [Pare [f a) 


(permutation de la dérivation par rapport à M et de l'intégration par rapport à P). Or : 


TELS 


1 
(re terme supplémentaire qu’on attendrait, d’après (1-20), soit : FM -divy P), disparaît parce que ÑP) n’est 


pas fonction des coordonnées du point M | De (1-34), on tire que l’intégration sur Q porte finalement sur : 


M (zm) SeN OR: + py dive GP) 


(On a réutilisé (1-20) mais de manière ici à faire apparaître les divergences calculées par rapport au point 
courant d'intégration.) A priori, les deux termes voulus par (1-20) sont présents; cependant div, (RP) est 
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DI. 


nul, d’après (4-8) et parce qu'on est en régime stationnaire, (c’est donc très différent de l’annulation, toujours 
vraie celle-ci, de div, ÑP)). Au total, div À est proportionnelle à 


[ (Il “ai, (EP). à, 


(PM) 


et, comme la divergence est calculée par rapport au point d’intégration P, on peut transformer par le théorème 
de la divergence en Je) -dS (où Z est la surface fermée qui borne 9). 
JVS 

Nousavonssupposéquelescourantssontlimités au domaine borné 9 ; c'est dire que les chargesr mobilesne 
traversent pas la surface X : entoutpointde Z, j's’annule ou est tangente à X. Il en est de même dej} PM, dont 
le flux à travers X s’annule. C’est aussi le cas de div A; d’où (5-8). 

(Ce calcul, évidemment un peu long, donne un premier exemple des manipulations de champs de vecteurs 
que requiert l’ Électromagnétisme. Ilimporte de bien posséder les résultats généraux : théorème de la divergence, 
formules (1-20) et (1-34), etc. De plus il faut prendre garde à rester maître des hypothèses physiques préalables : 
ici, la stationnarité des courants qui créent le champ magnétique. Le reste est affaire d'attention et de soin...) 

Remarquons pour terminer que dans le cas de courants limités à des circuits filiformes, on peut transfor- 
mer (5-7) en : 

I di, z 
Ho. Electromagn. 5-9 
=a $ 4n PM ( mE 


CALCUL DE rot Ë THÉORÈME D’AMPÈRE. 


Toujours pour des courants stationnaires, calculons rot B. De (5-6) et (1-33) nous tirons 
rot B = rot (rot À) = grad (div À) — AA., 


Si nous utilisons la jauge de Coulomb (5-7)-{5-8), il reste — AŽ. On a ainsi : 


rot, B(M) = — AyA{M) = >i 1 adh A dtp 


(une fois de plus, on permute l’ordre de la dérivation par rapport à M, et de l'intégration par rapport à P). 


Dans À (2 Aa): jP) est une constante vis-à-vis de M : on peut la sortir du Laplacien. Il reste à calculer : 
1 
-ge [lee Ad prr) 


= z | í | k JAP): [~ 4n 8U(P)] = Hoj (M). 


D’après (2-23), il vient : 


Un calcul analogue montrerait que les deux autres composantes de rot B sont de même proportionnelles 
aux deux composantes de même nom de j (tous ces champs de vecteurs étant pris au point M). D'où la 


relation fondamentale : 
rot B = uoj (Électromagn. 5-10) 


qui constitue la forme locale du théorème d’ Ampère. 
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(Celui-ci s’obtiendrait, sous sa forme « intégrale », en prenant les flux des deux membres de (5-10) à 
travers une surface S limité par une courbe C. On aurait ainsi : 


jera-ajfre 
S S 


Le théorème du rotationnel permet de transformer le membre de gauche en la circulation de B sur la courbe C. 
L'intégrale qui figure au membre de droite, représente, d’après (4-6), Pintensité Z qui traverse S. D’où légalité : 


é B. dM = mI (Électromagn. 5-11) 


Je 


Contrairement à (5-10), (5-11) fait intervenir le champ B en tous les points de C : il ne s’agit donc pas d’une 
relation locale. On retrouve uñe distinction analogue à celle qu’on a notéeentre (3-10) et (3-11).) Notons enfin 
que, de (5-6) et (5-10),nous tirons une équation reliant directement le potentiel vecteur à la densité de courant : 


rot B = rot (rot À) = grad (div Å — A4 = ug). 
Pour la jauge de Coulomb, le terme en div À disparaît : il reste simplement 


AA = —pj | (Électromagn. 5-12) 


dont on remarquera l’analogie formelle avec l'équation de Laplace (3-14). 


CHAMP D’INDUCTION MAGNÉTIQUE DÛ A UNE CHARGE PONCTUELLE EN MOUVE- 
MENT. 


La magnétostatique, concernée par des courants stationnaires, ne considère que les champs magnétiques 
totaux (5-2) ou (5-4) produits par l’ensemble des courants ou des circuits : c’est en effet à ce niveau que se mani- 
feste la stationnarité du régime de déplacement des charges. 

On dit par suite assez souvent que les champs élémentaires définis par (5-1) ou (5-3) n’ont pas de « réalité 
physique » : ilfaut reconnaître que, puisque les éléments de volume dx ou de circuit d/ ne peuvent être isolés 
du reste du système de charges en mouvement sans en interrompre le caractère stationnaire, on ne mesure 
jamais que les champs totaux (5-2) et (5-4). Il serait cependant tout à fait contre nature de supposer que le 
champ dB de (5-1) n’est pas effectivement la contribution au B de (5-2) des charges mobiles présentes dans 
l’élément de volume dr, (charges qui ont une réalité physique parfaite) et qu’en ajoutant des champs élémen- 
taires dépourvus de réalité, on obtienne, par le jeu de compensations miraculeuses, le « bon » champ total. 

Nous considérerons donc que les dB de (5-1) ou (5-3) représentent effectivement les champs des charges 
présentes dans les dt ou d/ considérés. 

On peut même aller plus loin. Raisonnons sur (5-1) et supposons que des charges q, de vitesse ÿ, au nombre 
de n par unité de volume, sont celles qui, dans dtp, transportent le courant. De (4-3) nous tirons 


3 Lo 4 A PM 
dB = LE (ndr 
4x PM We 
Or n dr, est le nombre de charges mobiles dans dtp. Cette proportionnalité de dB à ce nombre nous conduit 
à admettre que chacune des charges q, de vitesse D, placée en P (ou infiniment voisine de ce point) produit en M le 
champ : 


| Bu) = £e. ZA PM 


7 = (Électromagn. 5-13) 
rm PM 
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Effectivement, on verra ultérieurement que (5-13) représente une approximation du champ magnétique 
produit par une charge ponctuelle en mouvement valable tant que sa vitesse v reste très inférieure à la vitesse c 
de la lumière ce qui est toujours le cas en Électrocinétique. 

Remarquons que, comme l’a d’ailleurs suggéré la discussion qui précède, l'examen d’une charge ponc- 
tuelle isolée en mouvement nous fait complètement sortir du cadre de la Magnétostatique. Cependant l’analo- 
gie formelle de (5-1) et (5-13) nous permet de donner très rapidement certaines propriétés analytiques de b : 


b est à flux conservatif divb=0 (Électromagn. 5-14) 
b dérive d’un potentiel vecteur : 


(Électromagn. 5-15) 


Il est par contre plus délicat d'évaluer rot È et nous reviendrons ultérieurement sur ce point. 

(5-13)et (5-15) ont, de plus, l'intérêt de poser pour la première fois le problème des changements de repère : 
ces formules sont écrites pour un repère où la charge a la vitesse T; si nous passons à un autre repère mobile par 
rapport au premier (ce qui. ne nous est plus interdit, comme en électrostatique), © devient une autre vitesse ÿ” 
et les valeurs de b et de & sont elles-mêmes modifiées. Les modifications du champ magnétique dans un change- 
ment de repère seront abondamment décrites dans la suite : elles sont en fait plus complexes que ne laisserait 
supposer (5-13) ou (5-15) dont nous rappelons que ce sont des formules approchées. Mais nous voyons bien 
qu'il y a là un problème qu’on ne peut éluder. 
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{5-1} Un cylindre conducteur d’axe Oz, de rayon a, est siège d’un courant de densité J uniforme parallèle à Oz. 


Calculer le champ magnétique Ë et le potentiel vecteur À qu’il établit dans l’espace, tant à l’intérieur du 
cylindre qu’à l'extérieur. 


agnétostatique 


La symétrie du problème impose d'utiliser des coordonnées cylindriques d’'axe Oz, et qu’on notera p, et z. 

Les raisonnements de l’Exercice (2-1) peuvent être repris ici pour orienter les deux champs B et À 

Soit M le point où l’on cherche ces champs, P le plan (M, Oz), Q le plan perpendiculaire à Oz mené par M. 
Le système est invariant dans la symétrie par rapport à P : le pseudo-vecteur B(M) doit donc être confondu 
avec l'opposé de son symétrique par rapport à P; il est donc perpendiculaire à ce plan. Le « vrai » vecteur A(M) 
doit coïncider avec son symétrique par rapport à P; il est donc situé dans ce plan. La symétrie par rapport 
à Q change le fil en un fil identique mais parcouru par un courant de sens opposé; il faut procéder à une 
inversion du sens du courant pour récupérer le système initial, ce qui change le sens à la fois de B et de À. 
Par suite, B(M) doit coïncider avec son symétrique par rapport à P : c'est effectivement le cas. D'autre part, 
A(M) doit ici coïncider avec l’opposé de son symétrique par rapport à Q, autrement dit être perpendiculaire 
à ce plan. 

Au total Ë n’a de composante non nulle que Be fonction de la seule coordonnée p, et À wa de composante 
non nulle que 4,(p). 

L'application du théorème d'Ampère à un cercle d’axe Oz et rayon p, le long duquel la circulation 
de B est : 


2n 
| BP): p d9 = 2xpBj(p) 
0 


permet de calculer B;(p) (B(p) en abrégeant). Cette circulation est en effet égale à l'intensité traversant le 
cercle multipliée par Lo, soit 


a)p <a: 2xpB(p) = noj’ Ho ; Bo) ==: p 


; 2 
a i z a 
b)p>a 2xpB(p) = na?j,' uo ; Ble) = de S 


On en tire la fonction A(p) en appliquant la formule fondamentale : 
B = rot À. 


Plutôt que d'utiliser, dans ce cas très simple, les équations donnant les composantes du rotationnel 
en cylindriques, appliquons le théorème du rotationnel au circuit élémentaire ABC D schématisé ci-contre. Le 
flux de B à travers ce rectangle est, algébriquement, 


d’ = — dp : dz: B(p). 
Il est égal à la circulation du pctentiel vecteur sur le périmètre soit 


dz[Alp + dp) — A(p)] = de” dp -$$ 


D'où 


ôA 
m — Bl) = -55o pour p<a 
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PE 


k £d p° (en choisissant d'annuler le potentiel sur l’axe 


Dans l’intérieur du cylindre, on prendra A(p) = 
u conducteur). 
i Fextérieur, le potentiel vecteur devient 
PR oN 
-7 do (à) +K. 


| 


Da 


„a constante K est choisie pour assurer la continuité de A (qui est essentielle à l’obtention, par B = ro 
Pun champ magnétique dépourvu de singularités, tel que celui dont on dispose ici). On a ainsi : 


Alp) = - Hd à fi + 210g (£ e) (p > a). 


On remarquera qu’il n’est pas ici possible d'utiliser la formule générale (5-7) qui conduit à un résultat 
ivergent du fait que les courants considérés s’éloignent à linfini. Pour la même raison, les constantes 
“intégration figurant dans l’expression de À ne peuvent être définies en imposant à ce potentiel d’être nul à 
infini. (On voit d’ailleurs qu'il y présente en fait une divergence logarithmique.) 


5-2] Un solénoïde infini de section circulaire, d’axe z’z, comporte, par unité de longueur comptée suivant Paxe, 

.N spires jointives que parcourt une intensité constante 7. Son rayon est a. 

1° Rappeler l’orientation du champ magnétique B et du potentiel vecteur À, et la valeur, suivant les 
points de l’espace considérés, de B. 

2° Trouver la forme générale des potentiels 4 dont dérive 8 à l’intérieur, et choisir parmi eux celui qui 
vérifie la condition div 4 = 0 (jauge de Coulomb), respecte la symétrie du système, et reste fini dans son 
domaine de définition. 5 

3° Trouver la forme générale de À à l'extérieur; déterminer ce potentiel vecteur en lui imposant de 
satisfaire aux conditions de symétrie et à la condition de Coulomb, de s’annuler à linfini et de se raccorder 
sans discontinuité avec le potentiel à l’intérieur du solénoïde. 


1° L'exercice (2-1) a permis de montrer que Ë est partout parallèle à l’axe, l'exercice (2-2) que, s’il respecte 
ı symétrie du système physique, À est partout normal au plan « méridien ». En coordonnées cylindriques, ces 
hamps de vecteurs ont pour seules composantes non nulles B,(p) et Alp). 


On démontre d’ailleurs classiquement (+) que B est en fait partout nul à l'extérieur, et, à l’intérieur, 
niforme et de valeur B = po WT. 


> Appliquons le théorème du rotationnel au circuit élémentaire schématisé ci-contre. Le flux de Ë à travers 
[BCD est p dp dð - B. On l’égale à la circulation de A soit : 


(p + dp) d8A4(p + dp) — p d0 Aj(o). y 
vient finalement : XX AA W 40 p-p? 


ô ô 
Bp dp d8 = dp: de Le Ap) et z CoA(p)] = KoMT: p 


oit 
AXE Ado) = ET -o 


(1) Cf. par exemple nos fascicules Électromagnétisme (MP. PC, F° année). 
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(annulation de la constante d'intégration permet de ne pas introduire de terme en . qui rendrait À singulier 


(SUITE) 


sur z'z). La condition de Coulomb est automatiquement satisfaite. 


3 Onintégrera ici 


ô c 
=— [pA =0 d’où Alp) = —. 
do L° (p)] (p) x 
La continuité est réalisée pour : 
E hNi, où = HW a 
a ro i d'où Alp) 5 À 


On remarquera que nous n’adoptons pas, à l’extérieur du solénoïde, la solution la plus évidente : le 

+ . . => $ . . 4 i . . 7 
champ B étant nul, on pourrait prendre simplement Æ = 0. Ce faisant, on introduirait une discontinuité au 
niveau du solénoïde entre les potentiels vecteurs à l’intérieur et à l'extérieur; Æ cesserait d’être dérivable et 


B = rot À d'être défini. Or ce champ magnétique, s’il est discontinu sur le solénoïde, y demeure défini. D’où la 
nécessité d’adopter un potentiel vecteur plus complexe, comme nous l’avons fait. 


[5-3] Une spire circulaire de rayon a, de centre ©, d’axe Oz, est parcourue par un courant d'intensité 7. En un 
point M (OM = R, avec OM > a) elle crée un potentiel vecteur dont on admettra qu’il a la forme : 


où Mest le «moment magnétique » de la spire, soit : 


A = KI(RR?) (È: vecteur unitaire de Oz). 


1° Calculer les composantes cartésiennes de À, et ses composantes dans un système de coordonnées 
sphériques d’axe Oz. S $ 

2° Calculer les composantes du champ magnétique B(M) qui dérive de 4. 

3° Montrer que Æ dérive aussi d’un potentiel scalaire V qu’on calculera. 


1° Si x, y, z sont les coordonnées de M, on a immédiatement 
=W. A sfa. ES 
4x R? ” 4x R? 
et en sphériques : 
___ bo. M: sin ð 
? Ar R? 


pastaa -Ż)G a h 
i ôz 4 Rê) \R 4x R° 


B, = += + 


ÀA,=0 ; 4,=0 ; À 
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p -4 _ 0A, _ mer + 
Z ôx ôy An R RŽ R? RŽ 
_ Ho# 3z — R? 
4r R5 


B, et B, sont, toutes choses égales d’ailleurs, proportionnelles à x et y : B(M) est donc contenu dans le plan 
méridien (M, Oz) ; soit B, et B, ses composantes, en coordonnées polaires, dans ce plan. On a par exemple : 


R-B=R-B,= Hok. I f 3x27 + 3y2 + (322 — R3} 
4n RŠ 
soit 
HomM 2 cos 8 
B, = . ; 
An R? 
Ona alors : 


B, = B,cos 0 — B, sin 0 


(par projection sur Oz) 


d’où 


2 29 2 
PR ee 2 cos” 8 _ 3cos 8 — 1\ _ Mko sin 0 
R? R? 4x R 


et finalement 


3° On aura reconnu l’analogue du champ électrostatique d’un dipôle électrique : on est ainsi assuré qu’il 
peut prendre la forme 


B = — grad V. 
avec 
= -W p BoM. 2 cost 
ý ôR 4r R? 
o1 V , pM sing 
0 ROO m R° 


Ces deux équations ont bien une solution commune (définie à une constante près) : 
_ KoM cos à 
4n R? 


y 


[5-4] Un solénoïde infini « réel » est formé par un enroulement hélicoïidal d’un fil conducteur de section carrée 
; APSS a E E E E 
(côté a) sur un manchon cylindrique de rayon { R — zj enroulement est à spires jointives. On a a < R, 


et on désigne par œ langle que fait avec l’axe z'z du solénoïde la tangente en chaque point à l’enroulement 
hélicoïdal. 
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Le fil est parcouru par un courant de densité j: en chaque point du conducteur J est tangent à Phélice 
formée par enroulement, et garde un module constant. 

1° Montrer que le champ magnétique B créé par cet ensemble est superposition du champ B, d’un 
solénoïde infini « parfait » (c’est-à-dire sans hélicité) et du champ È, d’un système de courants dont on 
précisera la constitution. y 

2° Calculer B, et B, en tout point M de l’espace, en fonction de sa distance p à z’z, de j, 4, et a. 

3° Réexprimer B, et B, en remplaçant j par Pintensité totale 7 à travers le solénoïde et évaluer le 
rapport B,/B.. 


1° On peut (au moins en principe) calculer B par l'intégrale 


, A + ri + , , 
Si l’on décompose partout j en deux champs j; et ja, on pourra écrire : 


5 =e || (A EA RE TE: à 
T r? 


jL aF 3 _ Ho j aF 
rs ele 


B, -4f 


Ici, nous projetterons j sur l’axe (en jz) et sur la tangente au « parallèle » du solénoïde (en J,), avec : 


avec 


ji =j' cosa Ja = j'sin a 
En jı, nous retrouvons la distribution du solénoïde « idéalisé » de l’exercice (5.2); en jẹ, une distribution 


; ; f a wLa 
tubulaire, comprise entre les deux cylindres d’axe z'z et de rayons {| R + 3) uniforme et parallèle à z'z. Vu 


, id A , 4. 7 
la petitesse de w, j} apparaît comme une correction à /1. 


2° B, est uniforme, parallèle à z'z et égal à : B, = HoN T si I est l'intensité qui traverse le solénoïde et N 
le nombre de spires par unité de longueur. Avec un fil de section carrée, dont laire est a?, on a : I = ja?. 


(cf. fig. 5-3). On en trouve 


: : eni a 
Chaque spire occupe d’autre part, parallèlement à z z, une longueur : b = 
donc “M = 


par unité de longueur. Au total, il vient 


“j'a = Wa cos? a 


B,, champ d’une répartition « tubulaire » de courant, garde les propriétés de symé- 
trie du champ du fil cylindrique rectiligne de l’Exercice 5-1. En coordonnées cylindriques, 
il n’a de composante non nulle que B,,(p). On le calcule en appliquant le théorème d’Ampère fo 
à un cercle axé sur z'z et de rayon p. La circulation du champ magnétique vaut : 2rpB,(p). Fio. 5.6 
L'intensité (du courant tubulaire) qui traverse le cercle d'Ampère est : 


— nulle pour p < R — S (à l’intérieur du solénoïde), d’où B, = 0. 


2 2 
— égale à j, X (a(r + — {R — 3 } pour p>R +5 (à l’extérieur du solénoïde). (On pourrait 
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compléter par un examen de la traversée de l’épaisseur du solénoïde.) Il vient alors : 


B __Uf2Ra _ ug’ Rasing 
De Na 2 een e 


2xpB, = Hoja’ Ra et 


p p 
3° Avec 1 = ja?, nous pouvons poser : 
nn ae Lol cos? à 
— intérieur du solénoïde : B, = Fe et B, =0 
— extérieur du solénoïde : 
Ising R 
E E A 
a p 


< ; B ; j 
d’où en prenant pour B, sa valeur maximale (p & R), un rapport i ~ tga très petit devant l’unité avec 
4 


les hypothèses admises ici. 


[5-5] On considère deux fils cylindriques identiques à celui de l’exercice (5-1), parallèles, et dont les axes 
sont distants de d. Ils sont parcourus par des courants identiques, de densité f uniforme et correspondant, 
dans chacun des fils, à la même intensité totale 7, mais les sens des deux courants, dans les deux fils, sont 
opposés. a : 

Donner un potentiel vecteur Æ pour le champ magnétique Æ créé par ces deux conducteurs, en lui 
imposant de s’annuler dans le plan de symétrie qui les fait se correspondre. Peut-on, ici, utiliser la formule 
générale (5-7)? Quel résultat dorme-t-elle ? 


A est somme de deux champs À, et À, créés séparément par les deux conducteurs, et dont l’expression 
a été trouvée dans l'exercice (5-1). La superposition est simple puisque À,, À, — et par suite 4 — n'ont de 
composantes non nulles que parallèles aux axes des deux conducteurs. On devra néanmoins réintroduire dans 
les expressions de (5-1), des constantes d'intégration œ, et %3 qui auront ici à être fixées par la condition 


d'annulation imposée à À. Il vient par suite : 


“2 
A= A, +A = alt + 2 Log (=) — (e + 2 Log (2))| (pour p, et p, > a), 


soit 


À >z oe 


A l’intérieur d’un conducteur, le conducteur n° 2 par exemple, il convien- 
drait de remplacer : 


ia? 
À, = — Eol (a + 2 Logf2) 


Fc. 5.7 
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EXERCICES f 
de 


(SUITE) 


par : 
: 2 
A = = FE + B; 
avec la condition de continuité pour p, = a : 
12 e 
KoJa Hoja 
4 + B2 = g Qa 
ce qui donnera 
2 
u a 
À; = Éd (a? - p3) - Fo %2 


et pour le potentiel vecteur total : 


qu’on peut récrire : 


On obtiendrait, pour l’intérieur du conducteur n° 1, une expression tout à fait analogue. 
Le calcul direct par (5-7) : 


est ici possible. La présence de courants de sens opposés fait que les contributions aux potentiels dues aux points 
des conducteurs les plus éloignés du point où À est calculé se compensent : les dénominateurs r deviennent voisins 
et les numérateurs en j dr sont opposés. Cette circonstance fait disparaître les divergences notées en (5-1). (On 
peut aussi dire que l’ensemble des deux conducteurs est la limite d’un circuit rectangulaire dont la longueur 
croîtrait indéfiniment, et l’on sait que le calcul du champ magnétique d’un tel cadre est possible, et conduit 
à un résultat fini de même que pour le potentiel vecteur.) 

On obtiendrait ainsi un potentiel À’ parallèle à z’z (comme la densité de courant j qui permet de le calculer) 
et indépendant de la cote z du point M où on le calcule (du fait de la symétrie de translation du système de 
courants parallèlement à z'7). 4’ s’annule, de plus, dans le plan de symétrie des deux fils : chaque contribution 
élémentaire (jdt) du premier conducteur peut être ici exactement annulée par une contribution opposée 
venant du second. 

Or on sait que À’, qui conduit au même champ magnétique B que À, diffère de ce potentiel vecteur d’un 
gradient : 


+ + 


A = À + grad f. 
À' et À sont parallèles à z'z, de même que grad f; f n’est donc fonction que de z, de même que grad f. Or 


ce gradient est identiquement nul dans le plan de symétrie : il est donc nul dans tout l’espace, et À et A’ sont 
identiques. 
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FORCE DE LORENTZ. 


En électrostatique, une charge q placée dans un champ électrique (électrostatique) Ë subit une force 


F=q4 Ë (Électromagn. 3-3) 
Si cette charge est mobile, dans un repère où sa vitesse est V, elle subit une force de la forme : 
l F=qE+rA B) | (Électromagn. 6-1) 
appelée force de Lorentz, qui se décompose donc en : 
— une force indépendante de 7 : qË (Électromagn. 6-2) 
— une force proportionnelle à 7 : q? À B (Électromagn. 6-3) 


(6-1), (6-2), (6-3), relations obtenues à partir d'observations expérimentales, définissent deux champs de 
vecteurs : E et B. 

Si la particule est immobile, seule subsiste la force (6-2). Si l’on est dans une situation d'équilibre électro- 
statique, cette force (6-2) doit coïncider avec (3-3) et le champ Æ avec le champ électrostatique. Par suite nous 
conviendrons 

— d’appeler le champ E défini par (6-2) champ électrique 

(— d'appeler la force (6-2) : qË, la force électrique) 

— d'i imposer au champ électrique de se confondre avec le champ électrique de l’Électrostatique quand les 
charges qui interviennent dans le problème, c’est-à-dire q et les charges avec lesquelles elle interagit, sont immo- 
biles. 

Si d’autre part les charges présentes dans le système sont mobiles et créent des courants stationnaires, on 
constate que le champ B défini par (6-3) coïncide avec l'induction magnétique (magnétostatique) créée par ces 
courants. Par suite nous conviendrons : 

— d’appeler le champ E défini par (6-3) champ (ou induction) magnétique 

— (d’appeler la force q? À B de (6-3) force magnétique (*)) 

— d'imposer au champ Æ de coïncider avec l'induction magnétostatique (5-1) ou (5-3) quand on est en 
régime de courants stationnaires. 

Ces définitions ont, nécessairement, un caractère bien formel pour le moment : cependant c’est (6-1) 
qui nous permettra de trouver les transformations relativistes des champs électrique et magnétique. Pour le 
moment, nous ne pouvons tenter de calcul explicite. Mais ce qu’il importe de bien comprendre, c’est que le 
passage aux régimes variables remet en cause la définition même du champ électrique (tel qu’on l’a étudié en 
électrostatique) et du champ magnétique (de la magnétostatique); a fortiori l'expression mathématique de 
ces champs (champ de Coulomb ; induction magnétique de Biot et Savart) doit-elle être reprise : ce n’est que 
dans ces cas limites qu’elle réapparaîtra. 

De la même manière, les propriétés de E et B (la valeur des divergences ou des rotationnels par exemple) 
seront à réexaminer : ilest essentiel de ne pas prolonger indûment les relations qui les déterminent hors de leur 
domaine réel de validité, En fait, et comme on le verra bientôt, les équations de Maxwell, pour une part, servent 
à préciser ce que l’on conserve des résultats de l’Électrostatique et de la Magnétostatique quand on passe aux 
régimes variables les plus généraux. 

Notons enfin que (6-1) pose aussi le problème des changements de repère. Si onse restreint à des référentiels 
galiléens, f doit être invariante. Or la vitesse ÿ de la particule chargée change avec le référentiel : il doit donc 


(1) C'est cette force qui est à l'origine des effets mécaniques que subissent les circuits placés dans une induction magnétique, 
autrement dit, de la force de Laplace : par extension de langage, on pourrait encore l'appeler force de Laplace. 
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en être de même au moins de l’undes deux champs È et B. En fait, une modification de B seule est insuffisante. 
Comme nous le verrons, E et B prennent des valeurs différentes dans des repères différents. Les formes exactes 
des transformations accompagnant un changement de repère ne seront obtenues qu’après des développements 
bien longs : pour le moment, nous retiendrons l’affirmation d'une dépendance du « champ électromagnétique » 
(c’est ainsi que nous désignerons désormais l’ensemble champ électrique + champ magnétique), pour un 
système de charges donné, vis-à-vis du repère qui sert à décrire ce système. 


LE PHÉNOMÈNE D’INDUCTION ÉLECTROMAGNÉTIQUE. 


Nous rappelons très brièvement les résultats essentiels ($). 


1° Circuit mobile dans une induction magnétostatique : 


Soit (T) un circuit électrique (c’est-à-dire une chaîne de conducteurs) filiforme, que nous orientons, 
d’ailleurs arbitrairement. Il est plongé dans une induction magnétique Ë créée par des courants stationnaires : 
une induction « magnétostatique ». On le déplace dans ce champ d’induction. 

Un porteur de charge libre présent dans le conducteur, de charge q, (par exemple un électron) aura une 
vitesse de la forme (Ñ + #) où : 


— à est la vitesse de déplacement du porteur dans le circuit ; le circuit (1°) (r) -i> 
étant filiforme, # est porté par la tangente à (T°) KE 
— Test la vitesse d'entraînement du circuit au point où se trouve le porteur. 
En l’absence, supposée, de champ électrique, le porteur subit une force LA y 
magnétique : 


F= qü +a Bsq Ë, FIG. 6.1 


où Ë n Joue le rôle d’un champ électrique : nous l’appellerons, conformément à l’usage, champ électromoteur 
d’induction. Par définition, on appellera force électromotrice d’induction sur (F) la circulation de E,, sur (I) : 


| e= Ë, dM = Ga B)-at +Ÿ @ À B) dM 
a) 


D T) 
Le produit mixte qui figure dans la première intégrale est nul : 4 et áM sont colinéaires. Il reste 


ë= f (À dj Dis de (Électromagn. 6-4) 
T) 


où d®, est le flux coupé par le circuit dans le temps dz, du fait de son déplacement. Le champ B est à flux conser- 


_vatif (nous sommes en Magnétostatique) : (6-4) est équivalente à 


e= — (S) | (Électromagn. 6-5) 
myt 


3 dë 
où & est le flux de B à travers (T`) et (T) caractérise la variation de ce flux qui résulte du mouvement de (T`) 
myt 


dans le champ B fixe. 


(1) Pour plus de détails, cf. par exemple, dans la même collection, le fascicule Électromagnétisme, chapitre 19: 
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2° Circuit immobile ans une induction magnétique variable. 


Le même circuit filiforme (T`) est maintenant fixe. Par contre È varie dans le > temps : ce n’est plus une 
induction magnétique de la Magnétostatique. (Par exemple les courants qui créent B ne sont plus des courants 
stationnaires. Ou encore, les circuits qui sont le siège de ces courants, dont linten- 
sité est maintenue constante, sont déplacés par rapport à (T), d’où la variation (r) g 
de À au niveau de (T).) L'expérience montre alors que (T) est encore le siège d’un KE 
courant induit, traduisant le déplacement dans (T) des porteurs libres qu'il 
contient, et qu’on peut caractériser par une force électromotrice d’induction : 


nn (T) (Électromagn. 6-6) Fic. 6.2 


(où (d®/dr)nagtraduit la variation dans le temps du flux d’induction magnétique à travers (1°) qui résulte de la 


variation de B). 
e ne peut plus résulter de l’action sur les porteurs de la force magnétique : avec les notations du 1°, celle-ci 


est en g ^ B ((T) étant fixe, 8 s’annule) et sa circulation sur (I) est nulle. Si les porteurs libres sont entraînés 
autour de (T`) c’est qu’une force électrique est apparue, produite par un champ électrique Æ tel que 


e= $ Ē: dM = — (S) (Électromagn. 6-7) 
dt 
(T) mag 


(Si nous retenons la notion de champ électromoteur, nous voyons qu’il se confond ici avec le champ £.) 


Remarques. 
1° Dans le cas où la variation de ® est produite par un déplacement des courants qui créent Ë, ces courants 
restant stationnaires, l'expérience du 2° se ramène à celle du 1° par un changement de repère : dans le premier 


cas, on se place dans un référentiel lié aux courants qui créent B et le déplacement relatif de (T) et 
deces courants setraduit simplement par un déplacement de (T`) ; dans le deuxième cas, on choisit un référentiel 
lié à (T), et le même déplacement relatif de (r) et des sources de champ magnétique se traduit par un 


déplacement de celles-ci, et une variation de B au niveau de (T). Le phénomène « intrinsèque » (indépendant du 
repère) est le même dans les deux cas : c’est le déplacement relatif de (I) et du reste du système étu- 
dié. Il n’est donc pas surprenant qu’à cette même cause corresponde un même effet : l'apparition d’une f.e.m. 
d’induction donnée dans les deux cas par 


mais avec une interprétation différente dans les deux cas de l’origine de la variation de ® avec le temps £. 
L'expérience nous confirme donc dans l'idée avancée plus haut d'une dépendance du champ électromagné- 


tique, c’est-à-dire de l’ensemble (Ë , B),vis-à-vis du repère utilisé pour le décrire : 
— dans le repère où (I`) est mobile, E est nul, seul existe B. 
— dans le repère où (T) est fixe, À varie dans le temps et il apparaît un champ électrique. 


2° Divers effets peuvent concourir à la variation de ® avec le temps : leurs influences se superposent. 


. dë re RS Re e ; 
Si Pr la dérivée par rapport au temps qui décrit l’ensemble de cette variation de ®, la f.e.m. d’induction 


nous est finalement donnée par la loi de Faraday 


d® 


DT (Électromagn. 6-8) 
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3° Rappelons enfin un aspect qualitatif du phénomène d’induction électromagnétique, traduit par la 
loi de Lenz : il agit toujours pour s'opposer à la cause qui l’engendre. 


LA FORME LOCALE DE LA LOI DE FARADAY. 


Dans (6-7) ou (6-8), comme dans (6-5) nous avons fait intervenir le flux ® de l'induction magnétique Ë 


à travers le circuit (T). Or nous savons qu’une telle expression n’a de sens que pour autant que B est à flux con- 
servatif : c’est le cas en magnétostatique (d’où (6-5)), mais, a priori, nous ne sommes pas sùrs que la propriété 


se conserve dans le cas général des régirnes variables. Cependant toute l’étude expérimentale très succinctement 
décrite plus haut suggère très fortement qu’on peut extrapoler au cas général la relation fondamentale : 


divB=0, (Électromagn. 5-5) 


(sans quoi, la signification même du terme « flux du champ À à travers le circuit » serait remise en cause). 
Cette hypothèse étant faite, (S) étant une surface quelconque s'appuyant sur (T), (6-7) donne : 


Pa = -4| ff 5- d5 
a t (S) 


Dans le membre de droite, on peut permuter la dérivation par rapport au temps et l'intégration par rap- 
port aux coordonnées spatiales. Si 7 désigne la dérivée de l’induction magnétique par rapport au temps en 


un point fixe (d’où le recours au à des dérivées partielles), nous obtenons : 


pem (eines 
a) (5) (S) 


(grâce au théorème du rotationnel). (S) étant, comme (T), quelconque, l'identité des intégrales suppose celle 
des intégrandes. D’où la relation fondamentale reliant en régime variable les deux composantes du champ 
électromagnétique : 


rot È = —— (Électromagn. 6-9) 


Remarque. 


Le champ E que nous avons considéré jusqu’ici apparaît au sein d’un circuit conducteur filiforme. Nous 
pouvons très naturellement lever cette restriction et accorder à (6-9) une validité beaucoup plus générale. En 


particulier, Ẹ, satisfaisant (6-9), peut apparaître au sein d’un conducteur massif, où l’on sait que le phénomène 
d’induction électromagnétique se manifeste par des courants de Foucault. Ë peut aussi apparaître dans le 
vide : une charge ponctuelle isolée (un électron, un ion), initialement immobile, placée dans un champ B 
variable « sentira » ce champ électrique et se mettra en mouvement (ce ne serait pas le cas si seule intervenait 
la force magnétique, proportionnelle à la vitesse). 

Comme nous avons admis que Ë reste à flux conservatif quand on passe aux régimes non stationnaires, 
on peut continuer à définir un potentiel vecteur À avec 


se + 


B = rot À (Électromagn. 5-6) 
(6-9) donne alors : 
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(on a permuté les dérivations spatiales du rotationnel et le a D'où 


ôt 
rot (2 + S) = 0. 
4 Á ,,. , ; ; Ja 
Le vecteur E + E dérive donc d’un potentiel scalaire V; d’où : 
+ Lot 3 r 
E = — grad V — F (Electromagn. 6-10) 


formule qui donne la forme générale du champ électrique : à un champ de « style » électrostatique (dérivant 
d’un potentiel scalaire) s'ajoute un champ qui ne satisfait pas à cette propriété, et qui traduit l’intervention 
du phénomène d'induction; on pourra appeler cette composante ( _— 2A le champ (électromoteur) d’in- 
duction. 

L'ensemble nous apprend qu’une charge ponctuelle mobile, qui crée donc, d’après (5-13) un champ 
magnétique non constant, est l’origine d’un champ électrique dans lequel figure un tel champ d’induction et 
qui ne dérive plus d’un potentiel scalaire. En particulier, le champ électrique d'une charge ponctuelle mobile 
ne peut plus être un champ de Coulomb. Comme nous l’avions indiqué au début de ce chapitre, le passage aux 
régimes variables remet en cause les relations que nous pouvions considérer jusqu'ici comme les plus fon- 
damentales. 


ÉNERGIE POTENTIELLE MAGNÉTOSTATIQUE D'UN SYSTÈME DE COURANTS. 


Avant de poursuivre plus loin l’étude des régimes variables, effectuons un retour sur un problème de 
magnétostatique, mais pour lequel la connaissance du phénomène d’induction électromagnétique est nécessaire, 

De même que le rassemblement d’un certain nombre de charges met en jeu un travail, emmagasiné par 
le système sous forme d’une énergie potentielle électrostatique, de même l’établissement de courants station- 
naires dans un ensemble de conducteurs, à partir d’un état de spos des charges, nécessite un travail, qui est 
emmagasiné par le système sous forme d’une énergie potentielle magnétostatique. Pendant la phase d’éta- 
blissement du courant, si lente qu’elle soit, apparaissent des effets d’induction électromagnétique dont la 
prise en compte est nécessaire, 

Soit une particule de charge q, vitesse 7 participant au transport des courants étudiés. Dans un champ 


électromagnétique (E, B) elle subit la force de Lorentz (6-1) 
F=qgËE+8a B). 
Le régime restera stationnaire si l’extérieur exerce sur la charge une force $ = — F Dans le temps dż, 
la particule se déplace de 7 di, et la force extérieure fournit le travail : 
dW = F- ọdt = — gË-vdt 
(la force magnétique, perpendiculaire à la vitesse, ne travaille pas : la force extérieure fournit un travail opposé 


à celui de la seule force électrique). 
Dans un volume dr, on trouve n dr particules de ce type, et le travail total fourni par l'extérieur au niveau 


de drt est : 


-> 


dW = — ngÿ°Ë dr = -j É-dr:dr (Électromagn. 6-11) 
en faisant intervenir la densité de courant donnée par (4-3), 
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1° Sinous sommes en régime permanent, le champ électrique (6-10) se réduit à un champ « de type élec- 
trostatique », dérivant d’un potentiel scalaire V. (6-11) donne, pour valeur de la puissance totale fournie par 
l'extérieur aux charges en mouvement : 


S j) J- grad V dr (Électromagn. 6-12) 
e 


où l'intégration porte sur le volume auquel est limité le déplacement des porteurs de charge. Avec (1-20), 
l’intégrande de (6-12) devient : E . 
div (V - j} — V -div (j) 


Le dernier terme disparaît car on est en régime stationnaire et 
divj=0 (Électromagn. 4-8) 


Pae S 


où Z est la surface qui] limite Q; comme les charges ne franchissent pas £ J n’a pas decomposantenormaleàcette 
surface et le flux de j (ou de Vi) à travers elle est nul. Au total : 


ce (Électromagn. 6-13) 
dt 

Ainsi, en régime permanent, il n’y a pas de transfert d'énergie, au total, entre l'extérieur et le système 
de courants. Typiquement, nous considéroris ici des circuits dans chacun des- 
quels figurent un ou plusieurs générateurs et un ou plusieurs récepteurs (des résis- 
tances par exemple). Au niveau des générateurs une puissance P, > 0 est fournie 
par l'extérieur aux courants; au niveau des récepteurs, une puissance P, > 0 
est rétrocédée à l'extérieur (on doit compter dans P, l’ensemble de la chaleur appa- 
rue dans les conducteurs par effet Joule). D’après (6-13), P, = P, : on récupère à 
l'extérieur l'intégralité de l’énergie fournie aux courants. Un circuit électrique en 
régime permanent permet donc le transport et la transformation éventuelle de l’éner- 

gie, sans plus (ce qui présente d’ailleurs un intérêt pratique déjà considérable!). 
2° Supposons maintenant que nous soyons dans une phase d'établissement De 
du régime stationnaire ; nous supposerons cet établissement ducourant suffisamment 


lent pour qu’à un instant donné le champ d’induction magnétique À ou le potentiel 
vecteur 4 dont il dérive soient encore donnés, en fonction de la densité de courant J, par les relations linéaires 
de la Magnétostatique : (5-2) et (5-7). 

Dans ces conditions, nous supposerons que si KP) est la densité de courant au point de rayon-vecteur F 
une fois le régime stationnaire établi (pour les temps £ > T), la densité j(F, i, au même point, à un instant £ 
de la phase d’établissement du courant (soit pour 0 < ż < T) est en : 


JP, t) = À(t) JP) (Électromagn. 6-14) 


Il reste : 


R 


où la fonction scalaire de £ A(t) croît de O à 1 entre £ = 0 et t = T. 

Si B(F}et A(F) sont les induction et potentiel vecteur en régime permanent, les mêmes quantités, pendant 
la phase d'établissement de ce régime, et du fait de leur dépendance linéaire vis-à-vis de la densité de courant, 
vaudront : 


BE, 0) = AA BO ; AE, t) = Ale A) (Électromagn. 6-14) 
Reprenons (6-11), mais avec, maintenant, un champ électrique ayant la forme générale (6-10) : 
go — GA — > då 
Er, t) = — grad V RRAES, grad V — A() a 
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L'énergie fournie, entre les instants £ et (t + dt), par l'extérieur au système de charges mobiles vaut : 


dW= il -ira |- grad V — an haa 


Q 


= AD- de ||| AP: grd v ar 


k ao D de [RA AA ar 


(Ici encore Q est le volume où s'établissent les courants étudiés.) IH et AD sont les grandeurs relatives au 
régime permanent final; l'intégrale qui figure dans le premier terme n’est donc autre que (6-12) : ce terme est nul. 

Il reste uniquement le terme introduit par le champ électromoteur d’induction, dont la prise en compte 
s'impose : 


L'énergie fournie par l'extérieur pendant toute la phase d'établissement des courants, et qui constitue 
l'énergie potentielle magnétostatique Æ, de ceux-ci, est donc égale à : 


t=T 1=1 sn 
s,=w=| dW = HEC 
t=0 1=0 a 


| je Adr (Électromagn. 6-15) 


Soit 


3° Densité d'énergie magnétostatique. 


La densité de courant J et le potentiel vecteur A de (6-15) sont, rappelons-le, ceux d’un régime station- 
naire, de même que le champ magnétique B correspondant. Nous avons ainsi : 


ro B = ao (5-10) et J î= A roi B 
0 
Or: div (Ë À A = À ro B— B rot À (Électromagn. 1-28) 


D'où : 


nee div (B À Ā dt + [Bride 
240 JJJ, 20 JJJo 


Le domaine d’intégration de (6-15) peutêtre,sans inconvénient, étendu à l’espace entier (disons à une sphère £, 
de rayon R tendant vers l'infini, et qui englobe Q) : il suffit de proionger j par un champ de vecteurs partout 
nul hors de Q. 

La première intégrale de E, devient donc, grâce au théorème de la divergence : 


| [ë PEET 


Pour R suffisamment grand, B| varie au plus en 1 /R? et Zil en 1/R. Le flux de (B A À) sur £, dont l’aire totale 
est en R?, décroît donc en 1/R, et disparaît quand R tend vers l'infini. 


TL 
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Il reste finalement : 


E, = Il Zm dt (Electromagn. 6-16) 
(esp) 


Ce résultat nous suggère que la présence d’une induction magnétique B, indépendante du temps, se 
traduit par la présence d’une énergie répartie dans l’espace avec une densité 


2 
dE, _ (Électromagn. 6-17) 
dr 2 


Le para!lélisme de toutes ces formules avec celles de l’Électrostatique est évident : 

— dans les deux cas nous obtenons une première expression liée à la présence de charges fixes (resp. 
mobiles) dans un potentiel scalaire (resp. vecteur). Notons au passage, pour la première fois, la dualité entre 
les deux couples : 


e+p V 


densité de charge potentiel scalaire 


Hoj A 
densité de courant potentiel vecteur 


qui s’affirmera, puis se précisera, dans la suite. 

— Dans les deux cas, cette expression se transforme en une autre où c’est la présence d’un champ qui 
apparaît comme primordiale. (Il faut pareillement transférer le calcul des intégrales de la région où se trouvent 
les charges à l’espace entier, ou, en tout cas, à la région où règne un champ, électrique ou magnétique.) Comme 
on le verra plus tard, c’est la deuxième forme qui se généralise le plus directement aux cas non stationnaires. 

LUE aS Los &0E? B? 

— Les densités d'énergie ont des allures très voisines : Far : DS . (Notons cependant que les fondateurs 
de systèmes d'unités auraient été mieux inspirés en introduisant, Eet constante fondamentale, non pas Ho 
mais son inverse, disons Yọ = #5 '. Le parallélisme entre Électrostatique et Magnétostatique, tant au niveau 
des formules fondamentales de définition des champs, (3-5) et (5-2), que des densités d'énergie potentielle, se 
serait trouvé sensiblement renforcé.) 


EXERCICES | 6 | Première approche 
des régimes variables: 
l'induction électromagnétique. 


-> 


[6-1] Une spire plane de surface Set résistance À est dotée d’une normale positive de vecteur unitaire #. 
# tourne dans le plan xOy, autour de l’axe Oz, avec la vitesse angulaire uniforme œ. La spire est plongée 
dans un champ magnétique Æ qui, à un instant donné, est uniforme et parallèle au plan xOy, et qui tourne 
dans ce plan à la vitesse angulaire uniforme œ, en gardant une amplitude Z, constante. 

1° Calculér Pintensité Z(t) dans la spire. 
2° Montrer qu’un couple s'exerce sur la spire, dont la valeur moyenne dans le temps est non nulle. 
Donner brièvement la variation de cette valeur moyenne avec œ. 


1° Le flux ®(f) à travers la bobine vaut : Bo S- cos (#, B()). Les angles polaires de # et B(r) sont en 
ot et wat. D'où : 


S(t) = BaS cos (os — wt 
Apparaissent dans la spire une f.e.m. d’induction : BoS(wg — w) sin (w — w)t, et un courant induit : 


i= BE Co — w) sin (@ — œt. 
2° La spire prend ainsi un moment magnétique À parallèle à # et proportionnel à Set i : 


RTE 


M= + (ao, — ©): [sin (a, — okt]: R Z 


Placé dans le champ B, ce moment magnétique est soumis à un couple 
=AnB 


parallèle à l’axe Oz, sur lequel il a pour mesure algébrique : 


__ BaS? 
— R 


= 1B? 
La valeur moyenne dans le temps est l = ZR 


est même maximal au démarrage (pour w = 0, quand la spire ne tourne pas encore). Ce dispositif donne 
le principe des moteurs asynchrones; signalons qu'une description plus fidèle aurait dû prendre en compte 


la self de la spire, ce qui aurait assuré une décroissance moins rapide de l avec œ. 


r 


(o9 — œ): [sin (ù — ok]? 


(@ — w)}qui reste positive tant que w est inférieur à œw : le couple 


[6-2] On considère un conducteur semi-infini homogène, de conductivité 6, limité par un plan xOy. Il est le 
siège d’un courant sinusoïdal, de pulsation œ, uniforme et parallèle à Ox dans tout plan parallèle à xOy, 
mais dont l’amplitude varie avec z. Calculer la dépendance en z de la densité de courant. 


La seule composante de j est j, qu’on notera j(z, t), de la forme j,(z)e” , en dénotant la dépendance 


sinusoïdale en temps par le facteur e~". La loi d’'Ohm nous donne un champ électrique E de même orientation, 
avec 
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EXERCICES ière approche / 
des régimes variables; 
l'induction électromagnétique. (SUITE) 


De la loi de l'induction 
= & ôB 
ðt? 
et en nous restreignant au champ magnétique de régime stationnaire sinusoïdal, qui varie lui-même en e7 +°! 
nous tirons que B n’a de composante non nulle que parallèle à axe Oy avec : 


: = dE. = j0(2) —iot 
HS EE 
d’où : 
B, EN ; jou) Tiot, 
O'o 


Appliquons le théorème d'Ampère tel que nous l’avons rencontré en Magnétostatique (c’est-à-dire pour des 
courants continus et des champs magnétiques indépendants du temps). 


= -> 


rot B = Hoj. 
Il vient ici : 
| EP ôB JA) -i 
ot ne Due Q lot 
Hojo(z) € ôz ti do 


Ceci nous laisse une équation du second ordre à coefficients constants pour déterminer la fonction j(z) : 
jo + icoo)jo = 0. 
V2 


Les solutions sont, si nous remarquons que les racines carrées de i valent + (1 + -7 


La présence de ces exponentielles d'arguments complexes (qui ne sont ni réels, ni imaginaires purs ce qui 
indiquerait une dépendance sinusoïdale avec la cote z) doit s’interpréter comme traduisant la présence d’une 
variation sinusoïdale de jọ(z) dont l’amplitude varie exponentiellement; c’est ce que nous trouvons en repassant 


, , Lie ow MaS A ; ons , 
aux parties réelles; ainsi, eU*% avec à = + PE s'écrit e*'e", soit e“(cos x + isin x), et la partie 


réelle est la sinusoïde « amortie » ou « amplifiée », suivant le signe de « : 


e” COS x = exp +( EIF z) cos o.z 


Si le milieu conducteur s'étend vers les z < 0, comme nous devons exclure l’apparition de champs 
électromagnétiques ou de densités de courant croissant exponentiellement vers linfini, nous ferons B = 0 dans 
l'expression générale de jọ(z) qui se réduit ainsi à | 


j(Z) = exp (1 + i) FRE z=e"lrlè. ir 


où 
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EXERCICES Première approche 
i des régimes variables; 


l'induction électromagnétique. (SUITE) 


constitue la « profondeur de peau » pour le conducteur caractérisé par sa conductivité g à la fréquence œ. 
On voit que le courant parallèle à la surface se répartit en nappes successives, alternativement d’un sens, puis 
de l’autre. L’amplitude du courant décroît exponentiellement, avec une longueur caractéristique qui est préci- 
sément égale à ô : pour z = 6, j, est réduit à la fraction 1/e de sa valeur initiale. (Il est recommandé de vérifier 
l’homogénéité de ô à une longueur.) 

On voit que l'essentiel du courant se concentre au voisinage de la surface (c’est « l'effet de peau ») sur 
une épaisseur de l’ordre de ô Contrairement à ce qui se passerait en courant continu, on ne diminuerait pas 
la résistance du conducteur en augmentant son épaisseur au-delà de valeurs de l’ordre de 6. 


(N.B. Notre raisonnement suppose qu’on peut adopter telle quelle la forme locale du théorème d’Ampère 
valable en magnétostatique. On verra au chapitre suivant que cette transposition n’est pas rigoureuse : il con- 


viendrait d'ajouter à la véritable densité de courant j un terme complémentaire en e, le « courant de 
déplacement ». Notre calcul n’est donc valable que si ce courant de déplacement est négligeable vis-à-vis du 
courant de conduction j, c’est-à-dire si nous avons : 


= 


ÔE 


tor | = E00 IE| <[]=olË, soit eo < o. 


2 
On peut transformer cette condition en introduisant l'épaisseur ô : il vient : 8? < a ~ (5) ) où À est la 
longueur d’onde, dans le vide, des ondes électromagnétiques (cf. chapitre 9) de pulsation œ. L'exercice 


suivant montrera que cette condition est aisément satisfaite pour de bons conducteurs à des fréquences pas 
trop élevées.) 


[6-3] Évaluer l'épaisseur de peau, à 50 Hz et 1 MHz, dans du cuivre de résistivité p = 1,7 x 10 Q-m. 
Peut-on effectivement négliger le courant de déplacement? (Cf. Exercice 6-2.) 


Pour (= 2) = 50 Hz (4 = 6 000 km), nous avons 


#2 2p A ne - (D) 
4m x 1077 x 2nv \2r 1077 x v 2r 50 x 1077 


—2 
ô = Es x + 34 = 4,65 cm. 


soit 


Pour v = 1 MHz (10% Hz), nous trouvons ô # 0,33 mm (la longueur d'onde À est alors de 300 m). 
Dans les deux cas, la longueur d’onde reste très supérieure à l'épaisseur de peau, et l’approximation de 
l'exercice précédent est très valable. 


[6-4] Un solénoïde infini, comptant ^ spires circulaires par unité de longueur comptée suivant laxe, est par- 
couru par un courant d'intensité 7. Evaluer l’énergie magnétique contenue dans une tranche d'épaisseur k 
suivant laxe; le rayon des spires sera noté a. 
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HULIN. - Électromagnétisme — i 6 


EXERCICES 


des régimes variables; 
l'induction électromagnétique. (SUITE) 


cree 


Application numérique : évaluer l’énergie contenue dans une tranche de 19 em de long, 5 cm de rayon, 
où règne un champ de 5 Teslas. 


Il n’y a de champ que dans l’intérieur du solénoïde : ce champ vaut B = Ho Mi (cf. exercice 5-2), il est 


uniforme (et parallèle à l'axe). D’après (6-16) dans la tranche d’épaisseur A et section za? du solénoïde apparaît 
donc une énergie : 


1 


~ nah: B? nah: N? z 
Wa = a] 


Numériquement, et avec les chiffres proposés, nous obtenons : 


marx (5: 107° x (0,1) x 25 


no à T = 0,78 x 104 Joule 
x 4n x 


(Cette énergie n’est pas considérable, mais on peut remarquer que pour concentrer, dans le même volume 
la même énergie sous forme électrostatique, il faudrait y maintenir un champ E, uniforme, de module tel qu’on 


ait £? = T 2B?, soit E = cB = 3 : 10° x 5 = 1,5 10° V7! -m ce qui est beaucoup plus difficile, sinon 


impossible, à réaliser, par comparaison avec le champ magnétique, déjà considérable, sur lequel nous avons 
raisonné ici.) 


[6-5] On reprend le système de deux conducteurs étudié à lExercice 5-5, Calculer Pénergie magnétique emma- 
gasinée dans l’espace entre deux plans perpendiculaires aux fils et distants de h. 


(On s'appuiera sur Pégalité : 


2n 
| Log {(a? + b? — 2ab cos 6)!/?} d0 = 2x Log { sup (a, b) +) 
0 


Nous calculerons ici l'énergie magnétique W,„ par la formule : 


1 > > 
"= [ff Adr 


étendue à l’intérieur des cylindres conducteurs. On fera le calcul pour le seul conducteur n° 2, (les rôles des 
deux conducteurs sont parfaitement symétriques et leurs contributions à W, égales), et écrirons par conséquent : 


h a 2n + 2 z 
W, = | dz| do, | p40 de} 1 — £2 t 2Log|% }- 
i à ò 4 a? a 
2q? a a 2n 
= EE on | dpa P 1 — AE 2f apa | a8- p, Log (2) | 
4 o 0 o a 
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EXERCICES | 6 | Première approche 
des régimes variables; 
l'induction électromagnétique. (SUITE) 


cn 


La première intégrale donne immédiatement : 


NL a \ _ me 
CA2 dé) 2 


Dans le triangle 0,0, nous pouvons poser : 
P = É + d? — 2p,:d: cos 0 
et 


a CE À CU 


avec (à l’intérieur du conducteur) p, < d. La formule du texte nous permet-donc d'écrire la deuxième intégrale 


sous la forme : 
Se d 2 d 
4r | dp,’ pa’ Log = = 2ra“ Log al 
0 


Au total, nous aurons, en introduisant l'intensité = j : ra’): 


soit 


[6-6] 1° On rappelle qu’un champ magnétique B uniforme peut être dérivé d’un potentiel vecteur en : 


A0 = 4B a OM =B AF 


Da | mt 


(O étant une origine arbitraire pour le rayon vecteur F). Justifier ce résultat. 


2° Un électron (masse m, charge q = — e) décrit dans ce champ uniforme et constant dans le temps une 
trajectoire circulaire de rayon À, perpendiculaire à B. Calculer la vitesse angulaire œ, de Pélectron sur cette 
orbite. 


3° Le champ £ reste uniforme à chaque instant, mais on fait lentement croître son module, à orientation 


fixe. Montrer que l’électron est soumis à un champ électrique Ë dont on précisera les caractéristiques; 
en déduire la dérivée de la vitesse v de lélectron par rapport au temps. En supposant que la particule 


demeure sur une trajectoire circulaire, montrer que son rayon décroît, le flux de B à travers cette trajectoire 
demeurant constant. 

4° On cherche à maintenir la particule, malgré son accélération, sur une trajectoire de rayon constant. 
Montrer qu’on peut y parvenir en adoptant un champ Z non plus uniforme, mais fonction de la distance p 


à l’axe de la trajectoire circulaire (on supposera que Ë reste parallèle à l’axe de la trajectoire) et en 
imposant une condition qu’on définira entre le flux magnétique & à travers cette trajectoire, le rayon Æ 
de celle-ci, et le champ À, au niveau de l’électron (condition « bétatron »). 
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EXERCICES 


Première approche 
des régimes variables: 
l'induction électromagnétique. (SUITE) 


5° On supposera que B(p) est en Bop”. Déterminer # pour satisfaire à la condition bêtatron. L’en- 
semble des hypothèses faites en 4° et 5° vous semble-t-il cohérent? 


1° Le calcul demandé est celui de l’exercice (1-8). 


2° La force « magnétique » subie par l’électron est ici centripète et de module g-vB = qgRBw,. Elle assure 
une accélération centrifuge Rœ? avec : 4RB«w, = mRuw?, soit 


aB un eB\ 
w, = — en valeur absolue : — }. 
m | m j 
(C’est l'expression de la « fréquence cyclotron ».) z 
3° Il apparaît un champ électrique d'induction E=- 24 soit 


? 1 /0B\ 
E = 3 ()1 


(il s’agit partout d’un champ « tangentiel », dont la seule composante non nulle, en coordonnées cylin- 
driques, est en E). E reste donc sensiblement tangent à la trajectoire (à peu près circulaire, et qui, comme le 
champ B lui-même, n’évolue que lentement), et communique à la particule une accélération tangentielle : 


dv > eļôb , e dB 
Ma Ss Rails RE 
Nous en déduisons, avec : 
v = @,R 
dede kinen Lee 
2 m 2 
©, dR + 3R dw, = 0  R°«, = k (une constante). 
D'où rxR°B = & = Jen (le flux du champ magnétique à travers la trajectoire reste invariable : l'augmentation 
e 


de B est compensée par une diminution de R). 
4 On garde maintenant un rayon R invariable. Seules varient, avec B, les vitesses angulaire «, et 
linéaire v de l’électron. Le champ d’induction E a pour circulation, sur le cercle décrit par cette particule, la 


force électromotrice d’induction + Si le champ Ê garde la géométrie (champ tangentiel à symétrie 
cylindrique) des § précédents, il aura ainsi pour module : 
E = (2rR) ! - \d@/di| 


et entraînera, dans le temps dż, une modification dv de la vitesse du mobile en : 


La force centripète subie par l’électron est définie par le champ B, au niveau de l’électron, et le raisonnement 
du 2° nous donne une vitesse angulaire 
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EXERCICES Première approche 
des régimes variables; 


l'induction électromagnétique. (SUITE) 


D'où une nouvelle expression de dv: dv LE ‘dB, Par comparaison, nous obtenons la condition 
; m 
cherchée : 
1 dë 


dB, = 
2 nR? 


e 


On peut remarquer que /nR? constitue une sorte de champ moyen sur l’ensemble de la trajectoire, soit 
B,,, et que la « condition bêtatron » s'écrit finalement : 


Il faut que le champ cesse d’être uniforme et que sa variation moyenne soit le double de sa variation au 
niveau de la trajectoire électronique. 


5° Nous aurons 
R 
B, = Bo: R" et D — f 2rp dp ` Bop” 
0 


n+2 
soit ® = 2rB, 433 (formule valable pour n > — 1), et la condition du 4° nous donne alors : 


n+2 
(28, R E = 2B,R". 


n +2) nR? 
D'oùn + 2 = l,n = — let 
B 
Blo) œ =. 
(P) 


Ce champ n’est acceptable que s’il satisfait à la condition div Ë = 0. C’est effectivement le cas, comme 
on le voit en appliquant la formule (1-18). B n’a qu’une composante, B,, qui n’est fonction que de p. Reste la 
difficulté de la divergence de B(p) pour p = 0 : la variation en 1/p de B(p) ne peut être, au mieux, qu’approchée. 


[6-7] La sphère uniformément chargée de l’exercice (3-2) est animée d’un mouvement de translation uniforme 

de vitesse Ÿ. 

1° En traitant cette charge comme ponctuelle, calculer le champ magnétique et l'énergie magnétique W, 
qu’elle crée à l’extérieur de la sphère de rayon a. 

2° W’, est proportionnelle à v”; en l’égalant à l’énergie cinétique classique de la particule (dont 
la masse est m), calculer le rayon a. Quelle est la valeur numérique de a pour l’électron, dans ce modèle. 
(On prendra q = 1,6 : 107 !° Coulomb, m = 0,9 x 107 °° kg.) 

Que devient alors l’énergie électrostatique W, calculée en (3-2) ; et l’énergie W; présente à l’extérieur de la 
sphère chargée ? Que suggère ce résultat? 

Calculer maintenant lénergie magnétostatique totale créée, tant à l’intérieur qu’à l’extérieur de la 
sphère chargée, et sans refaire, comme au 2°, l’approximation d’une charge ponctuelle. Soit W, la valeur 
trouvée: sa substitution à W; dans le calcul de a, Wy et W, apporte-t-elle des changements notables? 


17 8 
Ho. OTi module #2 - qv“ aa i 


1° Le champ magnétique en M, la charge q étant en O, vaut : B = 
4n r° án r 


La densité d'énergie magnétique en M est : 
dW,, B _ Ho . 2. sin? 0 


dt is. 32r? r+ 
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EXERCICES 


Première approche 
des régimes variables; 
l'induction électromagnétique. (SUITE) 


L'intégrale sur tout l’espace extérieur à la sphère est : 


œ 2% 02,2 
Wa = | r? dr | 2r sin 8 - d0- n £ x (A), 
a 0 


L'intégrale en r vaut : 


Celle en 8 donne : 


# 
i 3 
Au total, nous obtenons : 
wW m4. LE 1 t 
3 a 32? 3 4x a 
2° Si nous posons Wp =; mv?, nous obtenons une estimation de a : 
a = 2 xX Ho pa 
3 nm 
Pour électron, il vient : 
— 1912 
a -2 x 077 LE = 1,9 x 10715 m 
9 x 


Ce « rayon classique de l’électron » est de l’ordre de grandeur des dimensions nucléaires. 
Avec 


Re FANS 
Ro a ? HE sn * a’ 
nous trouvons : 
- ot m = x mc 
2078 2 Ho 7 10 
W = a a C E 


Dans les deux cas, le résultat est voisin (mais non identique) à « l'énergie au repos » de la particule, telle que 
la donne, en Relativité, la formule d’Einstein : 


W = me. 


3° L'énergie magnétostatique totale W, sera commodément calculée, sans approximation, en intégrant 


sur la sphère chargée mobile le produit ic J À. 


sphère sphère 


7 Ho 
isg 


+ 


jest ici en p#; À, potentiel vecteur, qu’on tire de : 
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EXERCICES Première approche 
des régimes variables: 


l'induction électromagnétique. (SUITE) 
; : | 1 p dt 
est multiple du potentiel scalaire V = TE “———. Au total, nous aurons 
ĉo J sphère 4 
l oe >2_1 3; 4 
zX A= 5 (PEKo40Ë)* V 
v? 1 
= |—] x -pr 
z) 2 £ 
et 
iz v 
n- Marai a pE 
sphère 2’ sphère 2 
nii — o 3ko Pe 12 , us Ru 
soit très précisément : so ah Si c’est W que nous égalons à l'énergie cinétique, celle-ci prend pour 
rayon : a = e x bo. Nous aurons d'autre part : 
10 x a 
c 2 f1 1 
e =M] Wm = (mt) ra 
5 5 2 
wW, = 6 Wa = 72% 


Les ordres de grandeur sont conservés, mais la modification des facteurs numériques n’a pas suffi à 
rapprocher W; et W, de l'énergie d’Einstein mc’. 


[6-8] Un fil conducteur cylindrique, d’axe z'z, de rayon R, de conductivité 6, est placé dans un champ magnétique 


uniforme, parallèle à z'z, et variant sinusoïdalement dans le temps : soit B, — = B, ei, On négligera, 
dans tout cet exercice, l’influence éventuelle du courant de déplacement (ef. Exercices (6.2) et (6.3)}, 


1° I apparaît, dans le conducteur, un champ électrique d’induction ÈE Le 

a) déterminer son orientation 

b) calculer sa valeur maximale, et son déphasage par rapport à B. 

2° À E, est associé un courant de densité j. 1 € Courant de Foucault »). Quelle est la puissance 
moyenne P dépensée, par effet Joule, par ce courant dans une portion du conducteur longue de 1 mètre? 

3° En fait, j, crée un champ magnétique B,, lui-même variable, et qui engendre, par induction, 


un champ électrique È, Le conducteur sera supposé infiniment long et B, nul à l’extérieur du conduc- 
teur. Calculer B, et Ë, dans le conducteur. 

4° A quelles conditions aura-t-on El < 1 | et |B,| < |Ż,|. Calculer P et |, |/|B,] pour R = 1 mm, 
œw = 100 z, B, = 0,1 Tesla, e = 10° 


P La symétrie du problème impose d’utiliser des coordonnées cylindriques d’axe z'z, soit p, 6 et z. 
Pour le Le uniforme B,, nous respecterons cette symétrie en lui associant un potentiel 


vecteur À 1 (M) = = = À ^ OM (O étant un point quelconque de z’z pris comme origine) (ef. Exercice (6-6)) 


dont la seule oué non nulle est 4,, = 3 B, peT 
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EXERCICES le Première approche 
des régimes variables: 
l'induction électromagnétique. (SUITE) 


mr 


f ; i a 04, ,à ; 
Il en est de même du champ électromoteur d’induction : E, = — FE (E, est partout tangent à la surface 


du conducteur : il ne se pose pas, grâce à un choix convenable de À :. de problème de conditions aux limites.) 
La composante non nulle de Ë, est en : 
Eio = wA = “Fe emor 
déphasé de x/2 par rapport à Ë i 
2° j, est proportionnelle à Ë, Q, = cE). La puissance dP dépensée dans le volume dr entourant 
le point M est dP = j, : E, drsoit 


o? B? ire m B 


oE? dt = © : p’ cos? wt: dt 


(en revenant à des expressions réelles puisque nous raisonnons ici sur les carrés E? ou B? de grandeurs à variation 
sinusoïdale dans le temps). Dans un cylindre de hauteur h, disparaît par effet Joule la puissance 


h R 24 2p2 2p2 R 
Í dz: f de | p dû: PL cos? oi] = 2rh Le cos? œt x f p? dp 
0 0 o 0 


. nRê 
soit g ou? B? cos? wt, dont la valeur moyenne dans le temps est : 
P= nRtow? B$ 
16 


3° Le courant de densité j, se répartit par nappes cylindriques, toutes axées sur z'z, et dont chacune 


constitue l’équivalent d’un solénoïde. Il doit donc créer un champ magnétique À, parallèle à z'z, et fonction 
de la seule distance p à l’axe. 
Le théorème d'Ampère : 


nous donnera : 


qui s'intègre en : 


où K est la constante d’intégration par rapport à p. On sera assuré de l'annulation pour p = R en prenant : 


B;{p) = r a (R? — p°): Bo e™™ 


On retrouvera le même résultat, par un raisonnement plus « physique », en exploitant à fond l'analogie 
avec la suite de solénoïdes emboîtés les uns dans les autres. La nappe cylindrique comprise entre les cylindres 
de rayons p et p + dp sur une hauteur élémentaire dz, transporte une intensité i = j,(p) < dp dz. Sur une longueur 

ba , f , 1 : k de 
unité, comptée suivant l'axe z'z, nous trouvons M =I telles « spires » D’après les propriétés connues 
du solénoïde (infini), celui-ci maintient à l'extérieur un champ dB, nul, à l’intérieur un champ dB, uniforme, 
parallèle à laxe, et de module : dB, = noi ajil) dp. 
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EXERCICES Première approche 
j des régimes variables: 


l'induction électromagnétique. (SUITE) 


A la distance p de laxe règne donc un champ B,(9) qui est la superposition de tous les champs 
élémentaires de ce type, que créent les solénoïdes de rayon plus grand que p, soit : 


8,6) = mnt) w= | 


p 


[i Posa -B, el dp' = Eaa (R? e p?)Bo eTiot, 


C'est bien le résultat trouvé plus haut, avec, en plus, la justification du fait que B, est nul à l'extérieur du 


conducteur. E : 
(Remarquons aussi que la présence du champ B,, qui vient se superposer au champ variable B, du fait 


des phénomènes d’induction, confirme la loi de Lenz : B, (de même orientation que B,) a en effet un signe 
tel qu’il « s'oppose à la variation de B, ». En notations réelles nous aurons en effet : 
dB, 


B,®% B, cos wt Mag oB, sin wt 


; . ,, dB 
B, X B, cos (o — 3) = B, sin wt, de signe opposé à t) 


Nous tirerons Æ, de : 


Ē, a même orientation que Ë,. Sa composante non nulle est E,, avec 


Š (Enp) = inBA(o) = - jun B, e= lee = p) 


D'où 


2 3 
PAU) ee p 
PEzdp) = — ( E Bo e (rg a Z) 


(nous prendrons nulle la constante d'intégration pour que E,,, après division par p, reste fini sur l’axe z'z), et 
2 
ow 
En = — Fe Bo’ eTR? — p’). 


4 Nous aurons B,/B, et E,/E; < 1 pour : HoswR? < 1 soit R < ô (où ô est l'épaisseur de peau 
du matériau considéré, cf. Exercices (6-2) et (6-3)). 
Avec les données numériques de l'énoncé, il vient : 


P % 1,95 x 107? Watt:m°! 


2 
BB(* Horon. ) æ 107?, 
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7] GÉNÉRALISATION DU THÉORÈME D'AMPÈRE 


AUX REGIMES NON STATIONNAIRES. 
COURANT DE DEPLACEMENT. 
EQUATIONS DE MAXWELL. 


L 


IL. 
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REMARQUES PRÉLIMINAIRES. 


Dans notre quête des équations générales qui règlent le comportement des champs électrique et magné- 
tique, redéfinis par la force de Lorentz (6-1), nous avons, avec l’aide des observations expérimentales du 
phénomène d’induction électromagnétique, été conduits à admettre, fût-ce à titre d’hypothèse plausible, les 
deux premières relations : 


divB=0 (Électromagn. 5-5) 
rot Ë = — (Électromagn. 6-9) 


Comme on le voit, les champs Ē et B semblent devoir être déterminés par leur divergence et leur rota- 
tionnel; de telles grandeurs jouent d’ailleurs, comme on l’a vu, un rôle important dans l’ expression de théo- 
rèmes fondamentaux dans les cas limites de P Électrostatique et de la Magnétostatique. Ce qui nous manque, 
pour le moment, c’est donc des expressions de rot B et de div Ë. 

Nous savons cependant : 

— que dans le cas limite de l’Électrostatique, nous devons retrouver la forme locale du théorème de 
Gauss : p 


div Ë = re (Électromagn. 3-11) 
0 


— que dans le cas limite de la Magnétostatique, nous devons retrouver la forme locale du théorème 
d'Ampère 
rot Ë = moj (Électromagn. 5-10) 


Nous pouvons aussi remarquer que cette expression de rot Bne peut sûrement pas être transposée telle quelle 
de la Magnétostatique dans l’Electromagnétisme des régimes variables, comme nous avons admis que c'était 
le cas pour la relation div B=0. 
Prenons en effet la divergence des deux membres de (5-10). Celle de rot Ë s’annule (cf. (1-23)); il vient 
donc F 
divj=0 


Or cette relation est impérativement limitée au cas des courants stationnaires; dans le cas général il convient 
de la remplacer par léquation de continuité : 

div j + ôp/ôt = 0 (Électromagn. 4-7) 

Pour obtenir ce qu’on peut appeler le théorème d'Ampère généralisé, c’est-à-dire une expression générale 


— > . . . , . rA . . . [A [A 
de rot B, nous serons en particulier guidés par les considérations suivantes, relatives au champ électromagné- 
tique d’une charge ponctuelle. 


LE CHAMP ÉLECTROMAGNÉTIQUE D'UNE CHARGE PONCTUELLE MOBILE. 


N.B. On doit comprendre qu’on vise à suggérer une généralisation du théorème d’Ampère, et pas à la 
démontrer. Dans les raisonnements qui suivent, les expressions manipulées sont approximatives; à la limite, 
on pourrait dire que tout est faux dans ce paragraphe, sauf le résultat, grâce à Maxwell! (Il sera d’ailleurs 


GÉNÉRALISATION DU THÉORÈME D'AMPÈRE 
AUX REGIMES NON STATIONNAIRES. 
COURANT DE DEPLACEMENT. 
EQUATIONS DE MAXWELL. (SUITE) 


intéressant de revenir sur ces calculs lorsqu'on disposera des expressions rigoureuses des champs électrique 
et magnétique créés par une charge ponctuelle.) 

Soit une charge ponctuelle q, passant à l'instant d'observation (t = 0) à l’origine du système d’axes par 
rapport auquel le « point d'observation » M a les coordonnées x, y, z (on pose OM = F; OM = r). La vitesse 
de la charge en O est Y qu’on prend parallèle à Oz. 

Nous avons vu que, pour des vitesses v petites devant c, vitesse de la lumière, le champ magnétique b 
créé en M par q doit avoir la forme approchée (5-13). Il dérive du potentiel-vecteur : 


Ho q? 


à = I (Électromagn. 5-15) 
Calculons rot b, c’est-à-dire : 
otb = rot rot à = grad (div à) — A (Électromagn. 7-1) 
(7-1) se décompose en deux termes : 
r — Af = — Eqn: af) (ä na de composante, comme v, que sui- 
vant Oz; il en est de même de AG.) 
Or af) = — 4rô (2-21). Ce premier terme n’a de sens que comme distribu- 
tion; il vaut : x 
Lo ‘ q' TSM) (Électromagn. 7-2) 
Or qô(M) constitue l'équivalent d’une densité de charge p pour une charge Fic. 7.1 


ponctuelle de valeur q. Le produit gô(M) > 5 est donc lui-même l’analogue, tou- 
Jours pour une charge ponctuelle, d’un produit pô, c’est-à-dire d’une densité de courant j j (d’après (4-3)). Ce 
premier terme dans rot b, c’est donc le produit Hoj du théorème d'Ampère de la Magnétostatique. 


2° Mais il apparaît maintenant un terme supplémentaire : grad (div à). (Il n'existe pas en Magnéto- 
statique parce qu’on raisonne alors sur des potentiels vecteurs de divergence nulle.) Ici, on a : 


me? ft _ _ BW 2 
4r Ôôz\r än r? 


et grad (div 4) a pour composantes (avec (1-6)) : 


div à = 


— suivant Ox = | #222 . 3x2 
dr r5 
— suivant Oy = (us) ee L bon 7-3) 
T r 
Z aia 
— suivant Oz = | £2 6z- =y $) 
år r5 


On voit donc qu’une composante supplémentaire non nulle vient s'ajouter à la composante « magnétosta- 
tique » OA de rot Ë, dans cet exemple particulier, mais simple, de régime variable. 

Il nous faut cependant interpréter la valeur ainsi obtenue, de manière à disposer d’une expression de 
portée générale. 

Supposons que le champ électrique de la charge ponctuelle, bien qu’elle soit mobile, soit encore un champ 
coulombien; en M, il vaut : 


+. (Électromagn. 7-4) 
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(En fait nous savons qu'il s’agit là au mieux d’une expression approchée, puisque le champ électrique vrai 
comprend une composante, due à l'induction électromagnétique, qui ne dérive pas d’un potentiel scalaire, 
alors que c’est le cas de (7-4).) 

La charge q se déplaçant, la distance r et le rayon-vecteur F entre la charge et le point M varient au cours du 
temps, de même que le champ (7-4). Dans le temps dr, la charge a parcouru 


O0" = ÿdt 
et r devient 


D'où V'dt F 


Fic. 7.2 


4re, Lr? rt 
-4l a 2 | de 
ÅRE r? r5 


Avec 7 suivant Oz, il nous reste finalement : 


a(o) dE, _f{_ a \3vzy 
Ares) TŽ dt 


dE, _ G) i E _ 2 
dt Areg RE 


De la comparaison de (7-5) et (7-3) nous tirons que le terme supplémentaire apparu dans roi b serait : 


ôE 
EoHo Er 


(Électromagn. 7-5) 


(nous reprenons ici une notation de dérivée partielle pour rappeler que È est considéré à point d’observation 
fixe : sa variation dans le temps est donc « authentique », et uniquement due au déplacement de la charge qui 
le crée). E 

Outre la simplicité des calculs, avantage de raisonner sur une charge ponctuelle est que, comme tout 
courant est finalement dû au mouvement de charges ponctuelles, le résultat doit avoir une valeur générale. 
Souvenons-nous cependant que ni ë (5-15), ni E (7-4) ne sont des valeurs exactes : notre propos est seulement 
de justifier une expression que des expériences ultérieures, au niveau des prédictions qu’elle aura permise, 
permettront seules d'imposer. 


LE COURANT DE DÉPLACEMENT. LE THÉORÈME D'AMPÈRE GÉNÉRALISÉ. EXTENSION 
DE LA VALIDITÉ DU THÉORÈME DE GAUSS. 


Ce que nous suggère le calcul précédent, c’est qu’au terme « habituel » go j, vient, dans l'expression de 


fs =. | ; ig ÔE j 
rot B, en régime variable, s'ajouter la quantité &,u, =. Remarquons que la présence de ce terme nouveau 
> okoa; 
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IV. 


est effectivement limitée au cas non stationnaire : en régime permanent, les grandeurs électromagnétiques 


cessent de dépendre du temps et 6Ë [ðt s’annule. 
Nous aboutissons ainsi (mais c’est toujours une hypothèse!), à l'équation 


à > ôE 
rot B = ro + 0%) 


A la densité de courant « vraie », j, (correspondant à une conduction effective d'électricité par des por- 
ei 


(Électromagn. 1-6) 


teurs de charges mobiles) s'ajoute ce que Maxwell a baptisé le « courant de déplacement », de densité £o F 


(L'origine du terme est peu claire : en tout cas, ce courant ne correspond pas à un déplacement de charges.) 
Nous avions noté plus haut une difficulté liée à équation de continuité : prenons la divergence des deux 
membres de (7-6). Il vient : 2 


divf + £o div LE 0, 


soit, en comparant à l'équation de continuité (4-7), et en permutant les dérivations par rapport à t et par rap- 
port aux coordonnées : 3 a 

-= (eo div Ë) = P. 
Fto t 

Si nous acceptons d'ignorer l'influence éventuelle d'une constante qui apparaîtrait dans l'intégration par 
rapport au temps, nous retrouvons le théorème de Gauss (3-11), et les deux expressions manquantes, celle 
de div E et celle de rot Ë, nous sont simultanément suggérées par cette suite de raisonnements. (Remarquons 
qu’au chapitre précédent, nous devions, de même, adopter en bloc les expressions de div À et rot Ë, à moins 
de se résoudre à des complications considérables.) 

Notons bien au passage que nous gardons le théorème de Gauss, alors que nous devrons renoncer à ce que 
le champ électrique d'une charge ponctuelle soit, dans le cas général, c'est-à-dire 
quand la charge est mobile, coulombien. “ 

Ainsi, Si à un instant donné, nous calculons le flux du champ E dû 
à la charge q, à travers une surface fermée (S), nous trouverons 0 ou g/8, sui- 
vant que q est, à cet instant, à lextérieur ou à l’intérieur de (5). 

Remarquons en outre que si l’on démontre, en Électrostatique, qu'un champ 
électrique dû à une charge ponctuelle et satisfaisant au théorème de Gauss est 
nécessairement coulombien, la même démonstration ne peut être reprise dans le 
cas général, et ne vient doncpas invalider les conclusions précédentes (cf. Ex.7-1). 


TABLEAU RÉCAPITULATIF DES ÉQUATIONS DE MAXWELL. 


D’ après la théorie de Maxwell, les champs électrique È et magnétique B sont déterminés, dans le cas 
général, à partir des densités de charge p et de courant j, par les quatre équations fondamentales, ies Équations 
de Maxwell : 


div Ë = £ (3-11) (Théorème de Gauss) 

0 
divB=0 (5-5) (B est à flux conservatif) 
Tot Ë = — Z (6-9) (Loi de Faraday) 


rot Ë = u(i + o3) (7-6) (Théorème d'Ampère généralisé). 
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On peut leur ajouter l'expression de la force de Lorentz : 
F= qË +UA B) (6-1) 


(A elles seules, cependant, les équations de Maxwell et (6-1) ne nous permettront pas de définir la situation 
électromagnétique dans un milieu où sont présents des porteurs de charge, bien définis physiquement et qui 
déterminent p et j, mais que rien, jusqu'ici, ne permet de caractériser (alors qu’il est naturel de penser que 
leur charge, leur masse, les actions non électromagnétiques qu'ils subissent de l’extérieur, etc. contribueront 
à déterminer leurs réactions au champ électromagnétique.) De nouvelles relations doivent donc venir complé- 
ter les Equations de Maxwell, en fonction des propriétés spécifiques des milieux auxquels on les appliquera. 
De telles relations sont appelées « équations constitutives » : un bon exemple est la loi d'Ohm (pour un conduc- 
teur ohmique de conductivité y) : 


f= VE (Électromagn. 7-7) 


Cette question sera reprise au chapitre 11 (deuxième fascicule). 


EXERCICES 7] Généralisation du théorème d'Ampère 
aux régimes non stationnaires. 
Courant de déplacement. 

Equations de Maxwell 


[7-1] On montre, en Électrostatique, que, si le champ électrique satisfait au théorème de Gauss, c’est un champ 
de Coulomb, radial, variant en -7 Comment se fait-il que le champ d’une charge en mouvement puisse 


r 
satisfaire au théorème de Gauss sans être automatiquement coulombien ? 


Le champ d’une charge au repos est à symétrie sphérique, centré sur la charge immobile : toutes les 
directions de l’espace, à partir de celle-ci, sont en effet équivalentes. Le champ est par suite radial, et son 
module n’est fonction que de la distance r à la charge. Le flux de ce champ à travers une sphère de rayon r 
centrée sur la charge vaut donc 4xr? : E(r), où AU) est le module du champ à la distance r; d’après le théorème 
q 


de Gauss, il est aussi égal à 2. D'où E(r) = z: il s’agit d’un champ coulombien. 


| ES a , o, . RAN SAEN au ES La 

Si la charge est animée d’une vitesse Ÿ, la symétrie sphérique est réduite à une symétrie cylindrique d’axe & : 

toutes les conclusions précédentes tombent du même coup; le champ peut continuer à satisfaire au théorème 
de Gauss et cesser d’être coulombien. Nous verrons ultérieurement quelle devient sa forme exacte. 


[7-2] On considère un condensateur C, plan, de surface X. I est lentement chargé par un courant d’intensité 7 : 
la charge q des armatures croît donc suivant la loi dg — 7 di au cours du temps. Le courant 7 est 
quasi stationnaire : le champ magnétique B qu'il crée doit fort peu différer d’un champ magnéto- 
statique, en particulier au niveau d’un circuit d'intégration l entourant juste le fil d’amenée du courant, 
assez loin du condensateur. Sa circulation sur 1, d’après le théorème d’Ampère, doit être égale 
à l'intensité qui traverse toute surface s'appuyant sur F, multipliée par #ọ Or si la surface S est 
effectivement traversée par un courant, il men est pas de même de $. Montrer que lintroduction 
du courant de déplacement permet de lever le paradoxe. 


La circulation cherchée, si l’on se fie au courant de conduction qui traverse S, vaut ul. 

A travers S’ ne passe aucun courant de conduction. Entre les deux armatures existe par contre un champ 
électrique Ë variable avec le temps, et par suite un courant de déplacement. Simplifions le calcul en admettant 
que le champ du condensateur est uniforme, sans effet de bord, et par conséquent limité strictement à la région 
en regard des armatures, et que, dans cette région, S’ est plane et parallèle aux armatures. 

Si ø est la densité de charge sur les armatures, |E | vaut Z; son flux à travers S’ est oÈ = +. 
Le flux de = est 1a = L, et celui du courant de déplacement I. Nous retrouvons donc une circulation 
de B sur F èn ol. i i 

(On s’affranchirait des hypothèses supplémentaires introduites ici en remarquant que S et S$’ forment 
à elles deux une surface fermée ; si l'orientation de la normale positive à $’, définie par une certaine orientation de T, 


est extérieure à cette surface, celle de S est intérieure. Le flux sortant de E 


à travers l’ensemble S + S’, d’après le théorème de Gauss vaut q/£ọ ; (dans le I 
calcul précédent, seul intervenait le flux sur S’). D'où : r 
+ — -= — ao ee 
[es-es a 
(5) (5) 9 E: 
a à DES EE 
CRE ST 
KAT a ĉo gy C -q 
6) ®© 


[ee 
(5) i FIG. 7.4 
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énéralisation du théorème d'Ampère 
aux régimes non stationnaires. 

Courant de déplacement. 

Equations de Maxwell. (SUITE) 


se 


EXERCICES | 7 | 


(où Jest la densité de courant de conduction dans le fil d’amenée du courant). Au total, et à condition de faire 
intervenir systématiquement courant de conduction et courant de déplacement, nous voyons qu'il est équivalent 
de faire intervenir l’une ou l’autre des deux surfaces S ou S' pour évaluer la circulation de B sur T. 


[7-31 Un conducteur est limité par un plan P chargé et doté d’une densité uniforme de charge, variant au cours 
du temps, (£) En effet, les charges initialement « collées » sur P fuient pro- + 


diculaire à P ). ; 
1 LY 
3° Si, cependant, nous appliquons le théorème d'Ampère de la magné- 


gressivement ce plan : elles créent ainsi une densité uniforme de courant j, À H (c) A 
constante dans le temps, dans le demi-espace surmontant le plan P. (j est perpen- à | TA | 
1° Donner l'équation d'évolution de e(t). | 
2° Montrer, par des arguments de symétrie, que le champ magnétique B À ; 
créé par le courant j est nécessairement nul. 
tostatique à un circuit d'intégration (C), de forme quelconque, situé dans un FIG, 7.4 
plan parallèle à P, nous trouvons que la circulation de B est différente de zéro. Montrer que la prise en compte 
du courant de déplacement lève le paradoxe. 


1° Soit une surface S du plan P. A l'instant ż, elle contient la charge Q = ø(t)' S; à l'instant (t + dr), 
la charge n’est plus que Q + dO = o(t + dif} S : o a varié de do et dQ = S : do. Cette charge dQ est celle qu’a 
emportée le courant j dans le temps dr. D'où : 


dQ=-j Sd et =) 


c décroît donc linéairement dans le temps. 

2° Le champ magnétique (M) en un point M doit respecter la symétrie du système physique étudié : 

a) symétrie de rotation autour de la normale MH menée de M au plan P. B(M) doit donc être colinéaire 
à cette normale 3 

b) symétrie par rapport à tout plan H passant par MH. Le pseudo-vecteur B(M) doit se confondre avec 
opposé de son symétrique par rapport à H : il est normal à H, et parallèle au plan P. 

Ces deux exigences sont a priori contradictoires : on ne peut les réconcilier qu’en admettant que B(M) 
est nul, quel que soit le point M. 

3° Si nous appliquons le théorème d’ Ampère sous sa forme primitive, valable en magnelosiaudug (ce qui 
semble a priori justifié puisque le courant j est stationnaire) : 


rot B = Moj 
nous trouvons que la circulation de B sur une courbe C quelconque, plane, et normale aux lignes de courant, 
est égale au flux de j à travers C multiplié par Ho, SOitHo JS. Cette valeur non nulle contredit le résultat du 2°. 
En fait, il convient d'introduire le courant de déplacement, et de prendre : 


— 8 > 0Ë 
rot B = n + o) 


Ē est ici le champ électrique créé par le plan chargé avec la densité superficielle e(t). On démontre 


classiquement, en Électrostatique, et par application du théorème de Gauss, que ce champ est en 2, 
0 
Comme øg, il varie au cours du temps. D’après le calcul du 1°, on a : 
dË| 1 do j dE à 
= | = — et Ego = — j. 
dt &o dé Eo dt 
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Équations de Maxwell. (SUITE) 


Le courant de déplacement contrebalance ainsi exactement le courant de conduction et rot 8 peut être nul, 
de même que B lui-même. 


(Remarquons que l'expression de Ë («théorème de Coulomb ») est tirée du théorème de Gauss qui reste 
valable dans ce régime qui n’est plus électrostatique, ce qui justifie le raisonnement présenté ici.) 
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PASSAGE DES CHAMPS ÉLECTRIQUE ET MAGNÉTIQUE AUX POTENTIELS. 


Nous avons déjà démontré que les deux champs Ë et B peuvent s'exprimer en fonction des deux 
potentiels scalaire et vecteur V et À par les relations : 


B = rot À (Électromagn. 5-6) 
Ē = — grad V — t (Électromagn. 6-10) 


(5-6) traduit l'annulation de div B. (6-10) traduit alors la loi de Faraday (6-9), compte tenu de (5-6). Ces 
deux équations sont donc équivalentes à deux des équations de Maxwell. Nous verrons dans ce chapitre 
comment les deux autres équations de Maxwell peuvent être transformées en des « équations aux poten- 


tiels », autrement dit en des équations définissant directement À et V sans passer par les champs Ë et B. 


INVARIANCE DE JAUGE. JAUGE ET CONDITION DE LORENTZ. 


On a déjà remarqué que les potentiels À et V dont dérivent deux champs Ë et B donnés ne sont pas 
entièrement définis par ceux-ci. La «jauge» (c’est-à-dire le couple des deux potentiels) peut se voir 
imposer une condition supplémentaire qui achève sa détermination : nous décrirons ci-dessous la condition 
qui fixe la jauge la plus utilisée en Électromagnétisme: la jauge de Lorentz. Auparavant, deux remarques 
préliminaires : 


1° Si À, V) est une première jauge associée à un couple de champs (Ë, B) donné, un autre potentiel 
vecteur possible est : 


À' = À + grad À (Électromagn. 8-1) 
(où À est une fonction de point à valeur scalaire quelconque). On a en effet : 
rot À’ = rot (A + grad À) = rot À = Ë 


et, le gradient disparaissant sous le rotationnel, (5-6) est encore satisfaite. V doit alors devenir V” tel que : 


à — GN 
= _— yV ee anea 
E grad z 
soit aussi égal à : 
, ôf 
= grad y — ETS 


D'où la condition : 
aan n À 
— grad V -%5 — grad V’ — 2 - À (grad À) 
et 


grad V’ = ma (v-Ẹ) 


Une solution possible est : 
V'= V - Z (Électromagn. 8-2) 


(8-1) et (8-2) explicitent ainsi le passage d’une jauge à une autre. 


LA SOLUTION GÉNÉRALE 


DES ÉQUATIONS DE MAXWELL: 
LES POTENTIELS RETARDES. (SUITE) 


IL 


2° Même si l’on admet que différentes jauges peuvent être associées au même champ électromagnétique, 
encore faut-il s’assurer que, au niveau des phénomènes observables (actions électromagnétiques sur les 
particules chargées par exemple), les prédictions de la théorie sont insensibles à ces changements de jauge. 
Autrement dit, il doit y avoir « invariance de jauge ». En particuiier, la jauge de Lorentz se révèle particu- 
lièrement pratique, mais toute autre jauge devrait également permettre de décrire une situation électro- 
magnétique donnée. Nous ne ferons ici que signaler cette exigence. 


La jauge de Lorentz est définie par la condition supplémentaire 


| 
. > oV à 
| div À + &olto "ar = 0 (Électromagn. 8-3) 


qui généralise, au cas des régimes variables la relation (5-8) qui fixe la jauge de Coulomb. (Autrement dit, 
en régime stationnaire, la jauge de Lorentz se confond avec cette dernière.) Il sera intéressant de vérifier 


l’homogénéité de la condition de Lorentz (8-3) (cf. Exercice 8-1}. 
Si, comme c’est le cas pratiquement, les charges et courants du système étudié sont, tout au long de 


son histoire, maintenus à l’intérieur d’un domaine borné, on imposera de plus à V et À de s'annuler à l'infini, 
ce qui généralise les conditions traditionnellement imposées à ces potentiels en Électrostatique et en 
Magnétostatique. Il n’est pas toujours possible de satisfaire à cette condition supplémentaire (ce sera, par 
exemple, le cas de «londe plane » qu’on étudiera en détail plus loin) : on prendra soin alors de ne pas 
utiliser indûment des formules ou propriétés qui ne sont justifiées que par la condition d’annulation à 
l'infini. 


ÉQUATIONS AUX POTENTIELS. 


De (7-6) avec (1-33), (5-6) et (6-10) nous tirons : 
© rot B = rot (rot À) = grad (div À) — A4 


> êl — ôA 
= Hoj + ag ( grad V — 5) 


Soit en réarrangeant, et en permutant les dérivées par rapport au temps et aux coordonnées d'espace : 


+ 2Â d — 3 ed OV 5 
Ad — eono + = — Hoj + grad ai À + soto | (Electromagn. 8-4) 
De la dernière équation de Maxwell (3-11) il vient de même : 
z f — ô ; 
div E = div |- grad V — 4] = (Electromagn. 8-5) 
0 


que nous réécrirons, pour faire apparaître une analogie avec l’équation (8-4) : 


ôy p ôf, > ôv z 
AV — sgt = — — | div À + Eoo — Électromagn. 8-6 
V9 2 Eo A g q | EL 


He OV 
aux deux membres est ajouté : £o#o DE ; 
t 
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Nous voyons immédiatement que la condition de Lorentz (8-3) va nous apporter une simplification 
considérable de (8-4) et (8-6) qui deviennent : 


je 

-> A > z 

AA — solto D = — Hoj | (Electromagn. 8-7) 
Ca 4 p | 


| AV — eh = -> | (Électromagn. 8-8) 


Le passage des régimes statiques ou stationnaires au cas général des régimes variables se traduit donc, 
au niveau des équations définissant les potentiels, par l'apparition des dérivées secondes par rapport au 
temps, qui viennent compléter l’équation de Laplace (3-14) et l'équation (5-12) de la Magnétostatique. 
Ces dernières équations ne sont d’ailleurs que des cas particuliers de (8-7) et (8-8), obtenus quand la 
dépendance en temps disparaît. 

On notera au passage la simplicité de (8-7) et (8-8), et leur profonde analogie : c’est ce qui justifie qu’on 
parle d’une « esthétique » des théories physiques. Tous les développements ultérieurs montreront d’ailleurs 
que cette harmonie n’est en fait que la marque de correspondances très intimes : entre Aet V d’une part, entre 
les dépendances spatiales et temporelles de l’autre. 

On peut donner une forme plus condensée au système (8-7)-(8-8) en introduisant l'opérateur 
d’Alembertien, noté : 


2 
ie 
c* ôt : 
i (Electromagn. 8-9) 
a = £obo 
Il vient alors : 
OA=-pj OV=- E (Électromagn. 8-10) 


Remarque. Le système des équations aux potentiels, de la condition de Lorentz, et des relations de 
passage des potentiels aux champs est équivalent aux équations de Maxwell. Dans la dérivation qui précède 
en effet, nous n’avons fait que transcrire celles-ci, en utilisant, à chaque pas des identités. Donné les 
«sources » du champ électromagnétique (c’est-à-dire les fonctions p et j), on peut déterminer celui-ci en 
résolvant d’abord le système (8-3)-(8-7)-(8-8), puis trouver Ë par (6-10), B par (5-6). Ce faisant, aucune 
solution étrangère ne sera introduite. 


IV. ÉQUATIONS RÉSOLUES EN Ë ET À. 


Si l’on souhaite cependant effectuer une étude directe de È et B, on peut transformer les équations de 
Maxwell, qui mélangent ces deux champs, en des équations où ne figure plus soit B, soit E, autrement dit 


des équations résolues en Æ ou È. 
Partons de (6-9) : 


rot E = — ETS 
En prenant le rotationnel des deux membres, on obtient : 
ooer A2 a(S) = -2G 
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(7-6) permet alors de transformer le deuxième membre en : 


-ifu F aÈ) 
— ĝt Loj ĉoKo “5r ôt 


En regroupant, et avec (3-11), il vient finalement : 
Ci 
ôt 


2 


AE — colo DE Le grad p + Lo (Électromagn. 8-11) 
t Eo 


De la même manière, nous écrirons : 


2. 
wam 


rot (rot B) = grad (div B) — AB = — AB = rot nor + Eolto A 


sy CS à M ‘uen ĉ aB 
Bo TOt j + Eoo; (TOt E) = Ho TOt j + otto z| — Gr |: 


soit : 
> 028 — > s 
AB — oho z7 = — Ho fOtj (Electromagn. 8-12) 
t 


Dans (8-11) comme dans (8-12), cest encore le d’Alembertien du champ considéré qui intervient. Les 
seconds membres, par contre, ne sont guère susceptibles de se transformer de manière intéressante. 
Autre problème : (8-11) et (8-12) sont-elles équivalentes aux Equations de Maxwell? (cf. Exercice 8-2). 
Cette parenthèse faite, nous pouvons revenir aux équations fondamentales (8-7) et (8-8), et expliciter 
leurs solutions. 


V. POTENTIELS RETARDÉS. 


Soit une distribution de charge définie par la fonction p(P, t), une distribution de courant définie par 
le champ de vecteurs JP, t), l'une et l’autre variant avec le temps # et avec les coordonnées Xp» Yp Z 
du « point source » P, en restant confinées dans le volume #” limité par une surface fermée S 

En un « point d'observation » M, de coordonnées xy, Yy, Zy et à l'instant t apparaissent les potentiels 
V(M, à et AM, t) définis par les trois équations (8-7)-(8-8)-(8-3). Nous montrerons que l’on a : 


j (Électromagn. 8-13) 
t 
meee A 


Fic. 8.1 


Les potentiels définis par (8-13) sont les « potentiels retardés » : on remarquera en effet que ce qui 


détermine le potentiel à l'instant r, ce n’est pas la densité à cet instant, mais à un instant retardé | £ — E 
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fonction de la distance r entre le point source considéré et le point d'observation. Ces densités 
- PM ; i r 
pl P, t — =h ji P,t — ER sont dites densités retardées. 


Tout se passè comme si’les potentiels, au point et à l'instant d'observation, se construisaient à partir 
d'informations sur la présence de charges ou de courants aux points 
sources qui se propageraient à la vitesse finie c (et non de manière instan- 
tanée) : les points sources P}, Pa, Pae. «rayonnent » en permanence 
cette information autour d’eux, à l'attention de tous les observateurs M,, 
M, etc. ; mais chaque observateur ne reçoit ces informations qu'avec le 
retard dû à la propagation à la vitesse c. Il « se construit » alors un potentiel, 
scalaire ou vecteur, à partir de la densité, de charge ou de courant, qu'il 
connaît à chaque instant, c’est-à-dire à partir des densités retardées. 
A ceci près, la « construction » des potentiels à partir des densités se fait 
sur le modèle de l’Electrostatique ou de la Magnétostatique : (8-13) se 
distingue de (3-9) ou (5-7) par la seule intervention des retards. 

Remarquons d’ailleurs que (3-9) ou (5-7) sont cohérents avec (8-13) : 
dans les régimes, statique ou permanent, auxquels ces équations s’ap- 
pliquent, p et j ne dépendent plus du temps, et l’opération « introduction d’un retard » est sans effet. 

La simplicité de (8-13) est frappante; la possibilité d’une description très intuitive des retards (pro- 
pagation à vitesse finie de « l'information » sur la situation électromagnétique) est également prometteuse. 
Néanmoins les calculs qui suivent, et qui nous assureront de la validité de (8-13), nous convaincront 
des sérieuses complications mathématiques qui apparaissent quand on abandonne les régimes stationnaires : 
peut-être est-ce une incitation à l'introduction d’une description mieux adaptée à l’Électromagnétisme; la 
suite suggérera que c'est un des mérites de la Relativité Restreinte. 


FIG. 8.2 


DÉMONSTRATIONS RELATIVES AUX POTENTIELS RETARDÉS. 


Nous allons montrer que le potentiel scalaire donné par (8-13) satisfait à (8-8). (L'identité des formules nous assure 
alors que le potentiel vecteur est effectivement solution de (8-7).) Ensuite, nous vérifierons que les deux potentiels retardés 
obéissent à la condition de Lorentz (8-3). 


1° Remarques préliminaires : 


La difficulté des calculs qui suivent tient au fait que les densités retardées doivent être calculées, pour un point source 
donné et un instant d’observation donné, en tenant compte de leur dépendance vis-à-vis du point d'observation. Celui-ci 
intervient pour définir le retard qu’on impose aux dites densités. Il est essentiel ici de préciser par rapport à quelles variables 
on dérive : coordonnées du point d’observation, du point source, temps d'observation etc... 

D'autre part, on se souviendra que les densités retardées, de courant ou de charge, sont, en quelque sorte, des inter- 
médiaires de calcul, dépourvus de véritable réalité physique. (Des exercices (Exerc. 8-3 p. ex.) aideront à percevoir l'écart 
entre ces densités retardées et les « vraies » densités.) Il en résulte qu'il faudra s'entourer de précautions pour appliquer aux 


. PM 5 PM ; ; ; js ; 3 oies 
deux champs scalaire p| P, t — —- | et vecteur j | P,t — —- |des relations désormais familières, l’équation de continuité 
(4-7) par exemple. e c 

Les « vraies » densités p (P, t), j(P, À sont des fonctions de quatre variables x,, X3, X3, x, avec : 


p(P, t) = pri, Xz X3; X4) 


3 P Xi = Xp X3 = Yp X3 = Zp X4 = t l (Électromagn. 8-14) 
JP, t) = HET Xas X3s Xa) 


Les densités retardées sont alors représentées par les mêmes fonctions des quatre variables, mais les valeurs des x, ont 
changé; dans (8-13) interviennent : 


PX Xa X3 X4) PM 
avec Xi = Xp X3 = Yp X3 = Zp X4 = t — ——— 


r (Électromagn. 8-15) 


ir Xz, X39 Xa) 
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: ĉp 
Par suite quand nous calculerons -— par exemple, il conviendra de prendre garde à ce que, quand xp varie, toutes choses 


Xp 
égales d’ailleurs (c’est-à-dire pour Yp et Zp, M et t fixes), x, et x, varient. La règle de dérivation des fonctions de fonctions 


se traduit par : 
op _ Op ôx z ôp 0x4 


ôxp ôx, xp Ôx, xp 
ôx; 


Onax, = Xp et > 1 
De 
PM f ôx, 1 Ô(PM) 
x4 = t- — nous tirons : — = —-: 
c Xp €  Ôxp 

soit 2 

0x, RR 1 A (MP), 

xp c MP 


(avec (1-6), et en prenant M comme origine des coordonnées). 
Au total, il vient : 


M z 
Oox caeh = MP (Électromagn. 8-16) 
ôxp ôx, c ĉx, MP 


De la même manière, on voit que la « vraie » densité définie par (8-14) ne dépend pas de M, alors que c’est le cas de 
la densité retardée (8-15), comme nous l’avions déjà remarqué. Quantitativement, on peut par exemple écrire : 


ô êp ô ô 1 (MP ô 1 (PM i 
E A E 2(- ae = (- 1.0) (Électromagn. 8-17) 


(toujours d’après (1-6) mais en prenant ici P comme origine). 
2° Étude du potentiel scalaire : 


Dans (8-8), le potentiel est calculé au point M à l'instant ż; le laplacien est lui-même calculé par rapport au point M : 
c’est Ay. D'où à transformer l'expression : 


PM PM 
dP, p= 2) dP, t- z 
E ENE se ex i 
Mans, III? PM 4e, JJJ 7 ML PM | 
des) 
pl P, t- — 
2 c 


Il nous faut donc évaluer les dérivées du type : o ee Commençons par le calcul des dérivées premières. 
On a ainsi : xir PM 


(7 
pl P, t =- 
ô CP EM E EAE PE 
x PM A bx IPM) PM At tu) 


re 

A P, t — — 

ê? € PM\ 6? 1 ô 1 ô PM 
ol P,t go PS sus 

ôx? PM c / ôx? \PM xu \PM/ Oxy c 


1 ð? [e zn] 
T l 
PM ôx? g 


et par suite : 
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Avec les notations de (8-15), et en mettant à profit (8-17), on a : 


à = 
ô Arc 1) | _ _ 12 (PM) 
ÊXy c c ôx, PM 


Ale ele) 


è fè , ð 3 ; : ; 
On calculera rS (2) comme on a calculé l'origine de la variation avec x,, de la quantité qu’on dérive est à trouver 
x X A 


| see : . Op ý 
dans le retard et le résultat s'obtient simplement en substituant 2 à p dans (8-17). Il vient : 
X4 


et 


ð (2) : 1 (PM), p 


Axy \ôx4 c PM ôx? 


Rassemblant tous les termes nous pouvons écrire : 


a a 
1 (PME òp is (PM) à (1) Doa [E 


c? PM? ôx? c àx, 


PM  Ôôxy \PM. 
a f/1 1 (PM @p 16 ð ê` 1 ©? 
ox \PM/ oc? PM? ôx? côx, | ôxy Ôxy \PM/ PM ôx} 


ô 2 
(en réutilisant (8-16). (o2 et =E sont prises en P et, naturellement, « retardées ». } Le dernier crochet suggère 


ÔX4 4 i 
l'intervention d’une dérivée seconde d’un produit. (on sait que 


ð? ô? ôf à ð? 
~— (fg) vaut : 2 +2 f , 09 f: 1) 
ôx ôx ôx ôx ôx 


Nous pouvons appliquer ici à f = PM, g = 1/PM et prendre : 


ô ô f1 1 ô ô? 1 2/1 
2— (PM) Z) +—:— (PM) = Z (ru x Z) - PM— (>) 
xx ôxy \PM/ PM ôx} ôx? PM x2, \PM 


2/1 
-EG 
x3, \PM 


Un calcul analogue peut être mené pour les dérivées secondes par rapport à Ya et z,. Au total, nous aurons : 


1 
(ex x PM vaut 1 et sa dérivée est nulle |. 


M n 
(r - 2) 7 
c 1 1 ô 1 1ô 1 
PM PM/ œ 6x? PM c ôx, PM 


2 
: a P P i 
(avec, ici encore, des fonctions p, — et — « retardées » |. 
ôx, ôx? 
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Si nous effectuons l'intégration sur P, if vient donc : 


PM 1 
drp” dr, t- x) alza) 
4 c PM 
1 ií 1 ©? PM 
rs Île (e.-# 
ce? Ans, JJJ, PM ôx? c 


t 1 ôp 1 
—: d PM |: A À — 
i C AREg ii z Tp (2 ) a 


De (2-23) et compte tenu des propriétés de la « fonction » de Dirac, nous déduisons : 


PM 
a) que ła première intégrale vaut (— 4r) fois la fonction (2, t — e quand P est en M. Le retard s’annule alors; 
il reste simplement £ 
p(M, t) 

eo 


l pue ui tn x ô SE ; : 
b) que la dernière intégrale vaut de même (— 4x). (Zu) pour P en M, soit zéro, du fait de la présence du 
facteur PM dans l’intégrande. ôx, 

Si nous nous souvenons que la «vraie» densité p est définie par (8-14), nous voyons enfin que la dérivée 


seconde Ÿe qui figure dans la deuxième intégrale n’est autre que a LA dérivée seconde par rapport au temps, mais qu’il 


x4 
convient de calculer à Pinstant retardé. Cette intégrale s'écrit donc : 


JL &? P ð? 
: ZJ [i aeiee] 1 t). 
a or 4TEo PM c e at? 


En fin de compte, nous voyons que le laplacien du potentiel scalaire retardé (8-13) vaut 


(toutes valeurs calculées en M à Finstant f) et que ce potentiel, comme annoncé, satisfait effectivement à l'équation (8-8). 


3° Potentiels retardés et conditions de Lorentz : 


Il nous reste à montrer que les potentiels retärdés (8-13) satisfont à la condition de Lorentz (8-3). Dans ce cas seulement 
sera-t-on assuré de l’équivalence aux équations de Maxwell des équations aux potentiels, 
En ignorant un facteur multiplicatif, nous devons donc montrer l’égalité à zéro de : 


J P,t se êp| P, t 2 
c 1 A Cfpo (Électromagn. 8-18) 
[fa A PM ôt 
. pa 
Le calcul du premier terme de l'intégrande fait intervenir des dérivées telles que LE el 2), soit, en 
transposant (8-17) : Xy PM 
PM\ à 1 1 1 ð (PM 
(pr) (ee Hess . +0 À PM) 
c } Oxy \PM PM c ôx, PM 


En sommant les dérivées analogues obtenues par rapport aux trois coordonnées de M, on trouve : 


. ÑP, £ — PMjc) _P ÿ PM 
divy n- {P =). grady F 1 )-2 2. PM 


K PM/ c ôx, PM? 
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Le premier terme (cf. (1-34) et (1-20)) s'écrit encore : 


PM 
a 
a I c 1 PM 
— j- grad, | —— | = — di di i| P,t — — 
Pr (5) e| -PM TPM El Ne ] 
(Notre idée directrice, comme pour les manipulations analogues d’Électrostatique et de Magnétostatique, est de faire 
apparaître des opérateurs divergence calculés par rapport au point d’intégration, soit P.) 


S 


i disparaît de (8-18) car il se transformera en le flux à travers 4 du vecteur a = 
retardée, mais peu Wak ici) : les charges restent dans le volume Ÿ”, et j j (donc JPM) n’a, jamais, et nulle part sur #, 
de composante normale à cette surface, (le raisonnement est repris de la Magnétostatique, cf. chap. 5-II) et le flux 
considéré est nul. 

Reste donc à calculer : 


[lee a (ie e- -t as À 8x, T T + [ffe n aele PME) 


(Électromagn. 8-19) 


Le terme en div, (i (avec une densité 


\PM) 


Dans la première intégrale intervient la divergence de la densité de courant retardée, somme de trois termes du type : 
ô 1 TR 2 dj 1 êj (MP), 
6x} [7 c = ôx, c ôx} MP 


(il suffit de transposer 8-16). Le premier terme correspond à une dérivation par rapport à Xp mais à x, constant, done 
à'retard fixe (tandis que l'expression totale tient compte de la variation du retard quand P se déplace). On a donc 


wla PM m 1 a MP 
div EC meni z] = [divs J]e.:- 24) m z D MP 


Après division par (PM), le deuxième terme annulera la deuxième intégrale de (8-19) ce qui nous laissera 


1 ooa òp PM ; 
Mos Eri ai, VE, EM) + À (r. t — z) (Electromagn. 1-20) 


y 


expression dans laquelle intervient la divergence retardée de la densité de courant (quantité qu’il importe de distinguer, 
comme on le voit, de la divergence de la densité retardée). 

Bien qu’elle soit retardée, cette divergence de la «vraie» densité récupère sa signification physique habituelle; 
en particulier, elle satisfait à équation de continuité si on l’associe à la dérivée (par rapport au temps), retardée, de la densité 
de charge, quantité qui, elle, se confond avec la dérivée (par rapport au temps) de la densité de charge retardée. (En effet, 
les difficultés que nous rencontrons sont liées à des variations du retard avec la position de P ou de M. Ici P et M sont fixes, 


et le retard aussi : seul le temps # varie.) Or (1-20) associe précisément cette dérivée = à la divergence de j. L’équation 


de continuité (4-7) assure donc l’annulation de l’intégrande dans (1-20) et nous savons désormais que les potentiels 
retardés satisfont à la condition de Lorentz. 


EXERCICES La solution générale des équations 
de axwell. 
Les potentiels retardés. 


rar 


[8-1] Vérifier l’homogénéité de la condition de Lorentz (8-3). 


On sait que le potentiel scalaire V est homogène au produit par une vitesse du potentiel vecteur. (En régime 


stationnaire, nous avons ainsi 
1 dr > Ho dr 
ÂTE [E . aet pS ? 


j (cf. (4-3)) est homogène au produit (p x vitesse.) D'où|V/A X (okto) t: (elÿ) = 


Væ 


c? -= où c? est le carré 
(vitesse) 


d’une vitesse}. div A fait intervenir des dérivées spatiales premières de 4° E , .) : cette grandeur est, en 
Í i (A) ôV (v) ôV 1 (V) (V) . f 
ilM, eSt t RSS a H 
dimensions, ETET est en he et Eolo = zr en REET reproduit (A), le produit c x (temps) 
(4) 


constitue une longueur. On retrouve bien une expression du type ==, 
(longueur) 


[8-2] Les équations résolues en È et B, (8-11) et (8-12), sont-elles équivalentes aux équations de Maxwell ? 


(8-11) et (8-12) s’obtiennent, à partir des équations de Maxwell traduisant la loi de Faraday ou le 
théorème d'Ampère, en prenant le rotationnel des deux membres puis en effectuant des substitutions. 


Or il est évident que, B et Č étant deux champs de vecteurs, l’équation 


+ 


rot B = rot Č 


n’est pas équivalente à B = Č. Elle n’en est qu’une conséquence : elle implique en fait B=C+ grad f (où f est 
une fonction quelconque), ce qui est en particulier vrai si B = Č, mais constitue, en général, une relation 
différente. (De la même manière, deux fonctions F(x) et G(x) d’une variable dont les dérivées sont égales ne sont 
pas nécessairement égales : en général, elles diffèrent d’une constante. Un processus de dérivation, comme 
le passage d’un champ de vecteurs à son rotationnel, fait perdre une certaine information sur l'être mathématique 
de départ, fonction, champs de vecteurs, etc. Il conduit à des relations qui concernent des ensembles de 
solutions plus vastes que les solutions des relations initiales.) 

Une fois résolues (8-11) et (8-12), il faut soumettre directement les champs Eet B auxquels elles ont conduit 
aux équations de Maxwell pour s'assurer que des solutions étrangères ne sont pas apparues, ou les éliminer 
dans le cas contraire. 


{8-3} Un plan M contient une densité superficielle e de charge variant en 6 = 9, cos Ô où 6 est langle polaire (Ox, OM) 
du point M où e est mesurée. Le plan est mis en rotation à la vitesse angulaire constante œ autour de l’axe Oz qui lui est 
perpendiculaire en ©. 


1° calculer la densité de courant j qui apparaît alors dans J, en précisant ses coordonnées cylindriques. Calculer la 
divergence de cette densité, et expliquer le résultat obtenu. 
2° Si l’on cherchait à calculer les potentiels retardés V(O, £t) et A(O, t) qui apparaissent en © du fait de la présence 


107 


EXERCICES La solution générale des équations 
de Maxwell. 
Les potentiels retardés. (SUITE) 


A OP OP 
de e et de j, il conviendrait de faire intervenir les densités retardées dr. — (ee — 22), Donner l'expression, 
€ c 


résolue par rapport au temps £ et aux coordonnées polaires de P, de ces deux fonctions. 


° L'axe Ox du plan J a, par rapport à une direction OX fixe, un angle polaire qui varie linéairement avec le temps, soit 
@ = wt. Un point P de IT a donc, si 8 est son angle polaire dans Ñ par rapport à Ox, un angle polaire par rapport à la direction 
fixe OX : @ = 0 + p = 0 + ot. i 
. La densité de courant (normalement pọ où p est la densité volumique de charge, et D la vitesse des porteurs) devient ici 
j = cv (ü sagit, comme pour o, d'une densité superficielle). En un point P distant de p de l’origine O, la vitesse des charges 
est purement tangentielle et vaut œw : p. La seule composante de j est donc : 


j =0: w: p = (gw): p cos 0 


soit, en utilisant les coordonnées du repère fixe : 
Ja = (co@)p cos (© — ot). 


Nous calculerons ici div j en restreignant à deux dimensions l’expression générale (1-18); d’où : 


JS 1 ðh 
dvi == (pj) + =% 
u=, ap io) PA 
qui se réduit à 
di RC anO ct) 
Ness À . 


Rapportée elle aussi au repère fixe, la densité superficielle de charge est en 
o = go cos 0 = 6, cos (© — œt) 
D'où | 
ôo in (@ i) 
— = wo, Sin (© — ot). 
ôt £ 
i f £ sa fase p p A0 
Nous vérifions bien la relation de continuité | div j} + P 0, qui se réduit à div j + P 0}. 
2° Soit o*(P, t) = o*(p, ©, £) la densité de charge retardée au point p de coordonnées p et © (dans le système fixe), pour 
un instant { d'observation du potentiel V en O. Soit de même ;*(p, ©, t) la densité de courant retardée dans les mêmes conditions. 


ne > . se | à OP À ne Aon 
o* et j* se distinguent de ø et j par la simple substitution de l'instant retardé t — —— = t — à t. On a ainsi immédiatement : 
c c 


ot = oo'cos| 0 — oft -£)| = o:cosfo + En = ar) (#0) 
Z 


RSI 
et la seule composante non nulle de j* est de même p, À NE N 
© ` 
B= -oop sinlo+2p- ot] 7 NS 
| 7 NZ 


On remarquera que, tandis que ø garde même valeur, à un instant donné, en tous les 


do x ; Q 
points d’une droite issue de O (© = C"), o* prend mêmes valeurs en tous les points LL 
d’une courbe LÉ 


me pO” 19" maxi 
c . ; NS 
c'est-à-dire d'une spirale d’Archimède; ce sont en particulier deux telles courbes qui V (TC 
séparent les régions du plan où o* est, respectivement > 0 et < 0. On voit que la 
fonction o* n’a plus du tout l'allure de la fonction © : or c'est o* qui «crée» le 


potentiel V retardé. Fi. 8.3 
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EXERCICES 


{8-4] Montrer que les potentiels « avancés » 


PM 
de. e+ IM) 
se 
c dr» 
REg y PM 
e PM 
je, t + = 


+ Ho 
ġ =É ||| — à 
AS I$ PM Te 


V(M, t) = 


sont, au même titre que les potentiels retardés (8-13), solutions des équations aux potentiels (8-7) et (8-8). Voyez-vous 
une raison de ne pas les prendre en considération? 


Dans les équations aux potentiels, c intervient par son carré. On peut donc, dans l'expression des potentiels retardés, 
remplacer c par — c, puis dérouler la suite des démonstrations du chapitre et vérifier au bout du compte les équations aux 
potentiels (de même d’ailleurs que la condition de Lorentz). Dans cette transformation, les potentiels retardés se sont trans- 
formés en potentiels avancés, et nous voyons que ceux-ci satisfont à toutes les conditions imposées aux potentiels. 

(On peut donner une démonstration plus sophistiquée : les équations aux potentiels sont invariantes par « renversement 


du temps », car elles ne font intervenir que la dérivée seconde = si bien que le changement de ¢ en — # en respecte la forme. 


t , ; 
On peut donc changer t en — t dans la formule définissant, par exemple, le potentiel retardé : 


1 je 


Y(M, — t) = —— dt 
(Maa) ATE PM 2 
, : , + ; PM Lg à ; A 
puis « renverser » tous les temps : — ¢, qui redevient 4, et l’instant « retardé » — | t + —- j, qui devient un instant « avancé » 
c 


t + 2). Nous obtenons ainsi le potentiel avancé, à partir du potentiel retardé, par un processus qui ne compromet pas 
l'obéissance à l'équation (8-8). 
ai Ne E ; Gali s ôV 
En fait, les choses sont un peu plus compliquées : la condition de Lorentz fait intervenir la dérivée première -gr POY qu’elle 


reste satisfaite dans le renversement du temps, il faudra changer le sens (en même temps qu’il passe de « retardé » à « avancé ») 


du potentiel vecteur À. C’est somme toute normal : renverser le temps, c’est aussi changer le sens des vitesses (qu’on songe 
aux films comiques, projetés « à l’envers »); or, toutes choses égales d’ailleurs, le potentiel vecteur varie linéairement avec 
les densités de courant, elles-mêmes proportionnelles aux vitesses des particules, et changeant de sens avec elles, entraînant 
ainsi le retournement de À) 

Si les potentiels avancés ne sont pas, en général, pris en considération sur le même pied que les potentiels retardés, 
c’est que ceux-ci conduisent à une interprétation beaucoup plus satisfaisante de la relation (charge, courant)/(potentiels). Les 
charges et courants apparaissent en effet comme la «cause» des potentiels, et des champs gui en dérivent : il semble 
nécessaire qu’ils apparaissent avant les potentiels qu’ils créent, autrement dit à un instant retardé, et non avancé, par 
rapport à l’instant d'observation des potentiels. (Cette exigence est connue sous le nom de principe de causalité.) 

(Signalons que certains auteurs ont néanmoins proposé soit de tenir compte de potentiels avancés, soit de justifier par 
d’autres hypothèses leur élimination.) 


{8-51 Une plaque conductrice F7, d'épaisseur ô négligeable, plane, infinie et électriquement neutre, est parcourue par un courant 
de densité j(y, z, £) qui n’a de composante non nulle que j (y, z, f) = j, e71 
Calculer les potentiels retardés en un point M (aw’on prendra par exemple d’abscisse x, et situé sur laxe Ox perpen- 
dicuiaire à la plaque), et montrer qu’ils ne dépendent que de x et £. 
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(Ces potentiels ne peuvent être annulés à l'infini (car ici le système de courants n’est 
pas confiné à un domaine fermé de l’espace) : À, en particulier, n’est défini qu’à une constante 
près; c’est donc à sa variation avec x qu’on s’intéressera au premier chef.) 

Donner les composantes des deux champs ÊetË. 


Il n'y a pas de densité p de charge, donc pas de potentiel scalaire V. Seul reste à 


Aei) 
Le í 


évaluer le potentiel vecteur ĀM, t) == —— dr, (r = PM) Fia. 8.4 


J n’a qu’une composante à tout D] y 3 mi même pour À : nous calculerons A, L'élément de volume dtp, dans l’inté- 
grale (8-7) sera le volume associé à un élément d’aire dS de la plaque; celle-ci est d’ épaisseur ô très mince, et nous négligerons la 
variation du retard, ou du dénominateur r, quand P décrit l'épaisseur de la plaque. Avec dr, = ô : dS, il vient 


-iof ~$) 
AJM, i) = A fx, t) = F2 að f i P tas 
H 


Tous les points P d'une même couronne circulaire de centre O, comprise entre les cercles de rayons p et p + dp, sont 


à une même distance r (r? = x? + p?)de M, et affectés d'un même retard. La contribution d’une telle couronne à l'intégrale 
—io(t—r/c) 


est en (2xp dp). = et, nous trouvons finalement : 
ETE — 
A fx, t) = Ho "jo" à = dp-L. er 
en P r 


Prenant r comme variable (à x fixé, on a p dp = r dr), l'intégrale devient : 


La fonction intégrée n’est pas définie pour la borne supérieure d'intégration; ignorons cette difficulté, comme le suggère 
l'énoncé, et ne gardons que le terme fonction de la borne inférieure ; il vient : 


Hojo ` ô S 
Ax, t) = i 
da 2w 


~iolt = x/c) 


qui ne dépend que de x, comme l'exige d’ailleurs la symétrie de translation du problème parallèlement au plan /7 de la plaque. 
On en déduit un champ électromagnétique en : 


Hojoôe iot- re 
E, = se E lot ~ rc) 
B = — Holod | o- iot=rjo 
7 2 


(Nous verrons au chapitre suivant (deuxième fascicule) que la situation électromagnétique ainsi décrite est celle d'une onde 
plane.) 
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